Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 390,251 Δ14,565 bp [yaiT][dgcC] 16 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 390251 404815 14565 49 [0] [0] 49 [yaiT]–[dgcC] [yaiT],ytiB,insF2,insE2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,yaiZ,ddlA,iraP,phoA,psiF,[dgcC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3902500 (0.000)49 (0.820) 42/494 0.3 100% intergenic (‑523/‑1) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 404816 = 0 (0.000)coding (1114/1116 nt) dgcC diguanylate cyclase DgcC

GATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGC  >  NC_000913/404573‑405060
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
atatttattatgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataatggatggatcaacttattacgtccctgaggagggatgacaaTGACGTcc                                                                                                                                                                                                                                               <  1:138718‑M2/8‑1 (MQ=255)
  atttattatgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCgg                                                                                                                                                                                                                                             >  1:354910‑M2/242‑251 (MQ=255)
               atgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttctctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGACAGtcatc                                                                                                                                                                                                                                  >  2:443735‑M2/229‑249 (MQ=255)
                     caaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTg                                                                                                                                                                                                                          >  2:156419‑M2/223‑251 (MQ=255)
                      aaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGc                                                                                                                                                                                                                         >  1:154806‑M2/222‑251 (MQ=255)
                            aaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGtt                                                                                                                                                                                                                   <  1:207978‑M2/36‑1 (MQ=255)
                                ttatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTc                                                                                                                                                                                                               >  2:186246‑M2/212‑251 (MQ=255)
                                       caacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAtttt                                                                                                                                                                                                        >  2:275656‑M2/205‑251 (MQ=255)
                                                          gcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCg                                                                                                                                                                                     <  1:402082‑M2/66‑1 (MQ=255)
                                                                atgaatttgatgcgtgtttttatccataaagtttatatgtccacattcggacttagggtaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttgtctatgcattgataattgatggatcaacttattacggccctgaggagtgatgacaaTGACGTTCGGCGAAAGTCATCATGATTTGCTGAGTGTTTTTGATTTTGCCAtactcttcgtctgcggttcta                                                                                                                                                                               <  2:417033‑M2/72‑11 (MQ=255)
                                                                   aatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATc                                                                                                                                                                             >  1:218669‑M2/177‑250 (MQ=255)
                                                                   aatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATc                                                                                                                                                                             >  1:153614‑M2/177‑250 (MQ=255)
                                                                                     ctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCtt                                                                                                                                                          >  2:39362‑M2/159‑251 (MQ=255)
                                                                                                      gtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGt                                                                                                                                         >  1:40620‑M2/142‑251 (MQ=255)
                                                                                                         cacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAcc                                                                                                                                      <  1:186246‑M2/113‑1 (MQ=255)
                                                                                                         cacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAcc                                                                                                                                      <  2:212451‑M2/113‑1 (MQ=255)
                                                                                                          acattcggccttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCg                                                                                                                                     >  2:431421‑M2/138‑251 (MQ=255)
                                                                                                            attcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaactcattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCt                                                                                                                                   >  1:465221‑M2/136‑251 (MQ=255)
                                                                                                              tcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTc                                                                                                                                 >  2:53572‑M2/134‑251 (MQ=255)
                                                                                                                  acttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAAtg                                                                                                                             >  1:40249‑M2/130‑251 (MQ=255)
                                                                                                                  acttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAAtg                                                                                                                             <  1:294925‑M2/122‑1 (MQ=255)
                                                                                                                            aagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTc                                                                                                                   >  1:324381‑M2/120‑251 (MQ=255)
                                                                                                                                   attgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTg                                                                                                             >  2:332947‑M2/113‑250 (MQ=255)
                                                                                                                                         ccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGAc                                                                                                      >  2:450554‑M2/107‑251 (MQ=255)
                                                                                                                                           attcgtctgtaacgcagcataatcgttagcgctaaacataatatgttttctatgcattgctaattgatggttcaacttattacgtcactgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTGTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCACGGTGAGCAGAcc                                                                                                     >  2:445521‑M2/105‑250 (MQ=255)
                                                                                                                                               gtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAgg                                                                                                                                                                                                                               >  2:71405‑M2/101‑124 (MQ=255)
                                                                                                                                               gtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAgg                                                                                                                                                                                                                               <  1:71405‑M2/24‑1 (MQ=255)
                                                                                                                                                         tagcatactcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCATTGcc                                                                                      >  2:445871‑M2/91‑251 (MQ=255)
                                                                                                                                                          agcataatcgttagcgcgcaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGccc                                                                                     >  2:137199‑M2/90‑251 (MQ=255)
                                                                                                                                                            cataatcgttagagcgaaacataatatgttttctatgcattgatacttgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTTATCAGGATTTGCTGAGTTTTTCTGATTTTTCTATAAACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCAGGTGAGCAGACCATATCTTTCAGTGccccc                                                                                   >  2:436218‑M2/88‑251 (MQ=255)
                                                                                                                                                                  cgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGa                                                                                                                                                                >  2:42046‑M2/82‑168 (MQ=255)
                                                                                                                                                                  cgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGa                                                                                                                                                                <  1:42046‑M2/87‑1 (MQ=255)
                                                                                                                                                                   gttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGtt                                                                            >  1:137445‑M2/81‑251 (MQ=255)
                                                                                                                                                                        cgcgaaacataatatgttttctatgcattgataagtgacggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTGCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTccc                                                                       <  2:151883‑M2/176‑1 (MQ=255)
                                                                                                                                                                           gaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGGGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGtt                                                                    >  1:299523‑M2/73‑251 (MQ=255)
                                                                                                                                                                              acataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTcc                                                                 >  1:115183‑M2/70‑251 (MQ=255)
                                                                                                                                                                                         tttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAtttt                                                      >  2:117286‑M2/59‑251 (MQ=255)
                                                                                                                                                                                          ttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAtttt                                                      <  2:40620‑M2/193‑1 (MQ=255)
                                                                                                                                                                                              atgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCTTACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                  <  2:324381‑M2/197‑1 (MQ=255)
                                                                                                                                                                                              atgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCg                                                  <  2:218669‑M2/197‑1 (MQ=255)
                                                                                                                                                                                                gcaatgataattgatggatcaccttattacgtacctgaggtgggatgacaaTGACGTCCGGCGAATGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGAACCGCTTCAATGGTGGTGCTTCCCGGTGTGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTGTTGCGGa                                                >  2:289550‑M2/52‑250 (MQ=255)
                                                                                                                                                                                                 cattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                              >  1:83633‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                         attgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACgg                                                                                                                                                                 <  2:444691‑M2/86‑1 (MQ=255)
                                                                                                                                                                                                         attgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACgg                                                                                                                                                                 >  1:444691‑M2/43‑128 (MQ=255)
                                                                                                                                                                                                                   caacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGCTCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCAGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCgg                            >  2:397687‑M2/33‑251 (MQ=255)
                                                                                                                                                                                                                     ccttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCat                                                                                                   >  2:182199‑M2/31‑178 (MQ=255)
                                                                                                                                                                                                                     acttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCat                                                                                                   <  1:182199‑M2/148‑1 (MQ=255)
                                                                                                                                                                                                                               tccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGACCCCATTACCGCCTGAGCGGCAAATTTATacc                >  2:61587‑M2/21‑250 (MQ=255)
                                                                                                                                                                                                                                              gacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGc  <  1:111023‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGC  >  NC_000913/404573‑405060

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: