Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,912,832 6827 bp→TGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT [atpC][atpI] [atpC], atpD, atpG, atpA, atpH, atpF, atpE, atpB, [atpI]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3912832 3919658 6827 70 [0] [0] 63 [atpC]–[atpI] [atpC], atpD, atpG, atpA, atpH, atpF, atpE, atpB, [atpI]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39128310 (0.000)48 (0.860)
+TGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT
35/376 0.5 100% coding (390/420 nt) atpC F1 sector of membrane‑bound ATP synthase, epsilon subunit
?NC_000913 3919659 = 0 (0.000)coding (30/381 nt) atpI ATP synthase, membrane‑bound accessory factor

GAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAATGTGCTCTTCAGCCTTACGTTTCGCTTCCATGGCTCGCGCTTCGTCGAGATCCTGGCCGCGAATTGCGGTGTCGGCCAGAACGGTCACGTTGCCAGGCTGCACTTCAAGAATGCCGCCAGACAGATAGATAAACTCTTCGTGACCGTGCTGTTTCACGATGCGAATCATACCAGGCTTAATGGCGGTGAGCAGCGGTGCGTGGCCAGGGTAGATCCCCAGTTC  >  NC_000913/3912531‑3913127
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        
gagGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTTCATCGCTTTTTTGGTCAACTCGATAAcg                                                                                                                                                                                                                                                                                                          >  2:366524/1‑301 (MQ=255)
              cAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagc                                                                                                                                                                                                                                                                                             <  2:645250‑M1/300‑14 (MQ=255)
              cAATAACGGATTAACCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagc                                                                                                                                                                                                                                                                                             <  1:460412‑M1/300‑14 (MQ=255)
                aTAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagc                                                                                                                                                                                                                                                                                             <  2:37933‑M1/298‑14 (MQ=255)
                                 cTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                          <  1:147346‑M1/300‑33 (MQ=255)
                                   gACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatacttt                                                                                                                                                                                                                                                                       >  2:203660‑M1/1‑266 (MQ=255)
                                       aaaaCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctag                                                                                                                                                                                                                                                                   <  1:620329‑M1/301‑40 (MQ=255)
                                       aaaGCGCCATGTTCGACCATCGTCACGACTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGGTGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAATCAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactt                                                                                                                                                                                                                                                                        <  2:558829‑M1/296‑35 (MQ=255)
                                          gcgcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggtgt                                                                                                                                                                                                                                                                                             >  1:57056‑M1/1‑259 (MQ=255)
                                          gcgcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagc                                                                                                                                                                                                                                                                                             <  2:57059‑M1/272‑14 (MQ=255)
                                           cgcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggatggagctgcttcgaagttcctatactttctagagaa                                                                                                                                                                                                                                                               >  2:177655‑M1/1‑258 (MQ=255)
                                               aTTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATCTCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggatctgcttcgaagttccttttctttctagagaatag                                                                                                                                                                                                                                                            >  1:605104‑M1/1‑254 (MQ=255)
                                                                cGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggt                                                                                                                                                                                                                                          >  2:571114‑M1/1‑237 (MQ=255)
                                                                  ccTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcg                                                                                                                                                                                                                                        >  2:499928‑M1/1‑235 (MQ=255)
                                                                  ccTTATGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtc                                                                                                                                                                                                                                         <  2:613346‑M1/300‑66 (MQ=255)
                                                                            gCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaacttcaggtcgacggatcc                                                                                                                                                                                                                                >  1:492369‑M1/1‑225 (MQ=255)
                                                                                     tAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATAAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTAAACGCTATTCGGGACAAGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAGGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtagacggatccccggaataa                                                                                                                                                                                                                       <  1:448746‑M1/299‑84 (MQ=255)
                                                                                         gCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgcagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtt                                                                                                                                                                                                                   >  1:107599‑M1/1‑212 (MQ=255)
                                                                                              aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtc                                                                                                                                                                                                                                         <  2:242707‑M1/272‑66 (MQ=255)
                                                                                              aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggt                                                                                                                                                                                                                                          >  1:242697‑M1/1‑207 (MQ=255)
                                                                                                        ataagagacagggatgtttttgctCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataac                                                                                                                                                                                                                      <  1:589163‑M1/270‑85 (MQ=255)
                                                                                                                ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaaattcgaactgaaggtcgacggatccccggaataacgtttcgactcccgagcgaaacag                                                                                                                                                                                               >  2:285212‑M1/1‑189 (MQ=255)
                                                                                                                 gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTTTTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatca                                                                                                                                                                                         >  1:428300‑M1/1‑188 (MQ=255)
                                                                                                                   ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGGGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcaacgcatccccggaataacctgtctcttatacaactctcctagcccacg                                                                                                                                                                                        >  2:589177‑M1/1‑186 (MQ=255)
                                                                                                                       gTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtttaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcccgtttt                                                                                                                                                                                   >  1:385712‑M1/1‑182 (MQ=255)
                                                                                                                           ttGCTCTGATTTCGATCAATGAATAAACAGGAAGACAGGGGGGTTTTATATCAGAAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTTGTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcac                                                                                                                                                                               <  1:366512‑M1/301‑124 (MQ=255)
                                                                                                                                                 aTAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACTtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatcacgttttcactcaagaacgaacacgacatcacgtttttctctcctgctcccttcctggt                                                                                                                                                             >  1:489499‑M1/1‑156 (MQ=255)
                                                                                                                                                                   gggTTTTATATCAGCAGGATCTATGTGAACGCTAAGCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaac                                                                                                                                        <  2:107608‑M1/300‑163 (MQ=255)
                                                                                                                                                                                   ggATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacttcacgtctttcaatcctgctcccttcgaggtatgccgcgttgcgtataaaactttcttttatgatta                                                                                                                         >  1:619477‑M1/1‑122 (MQ=255)
                                                                                                                                                                                   ggATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcaacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttag                                                                                                                        >  2:104167‑M1/1‑122 (MQ=255)
                                                                                                                                                                                    gATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgc                                                                                                                                                         <  2:497561‑M1/266‑146 (MQ=255)
                                                                                                                                                                                     aTCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagag                                                                                                                      <  2:481627‑M1/300‑181 (MQ=255)
                                                                                                                                                                                      tctagtttaacgctattcAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTTTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagt                                                                                                                     <  1:203651‑M1/292‑182 (MQ=255)
                                                                                                                                                                                      tCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagag                                                                                                                      <  1:177646‑M1/299‑181 (MQ=255)
                                                                                                                                                                                       cTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagt                                                                                                                     <  2:385724‑M1/299‑182 (MQ=255)
                                                                                                                                                                                        taaaagagacagattcaggacgtgtCACACTCGCAAAAAAAAGCCAGCCGGTGTCCAGACTGGCTTTTGTTCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGggtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggc                                                                                                                            <  1:655408‑M1/279‑175 (MQ=255)
                                                                                                                                                                                         aTGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaa                                                                                                                  >  2:649264‑M1/1‑116 (MQ=255)
                                                                                                                                                                                         aTGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                          >  2:454601‑M1/1‑116 (MQ=255)
                                                                                                                                                                                         aTGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                          <  1:454588‑M1/148‑33 (MQ=255)
                                                                                                                                                                                         aTGTGAAAGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagag                                                                                                                      <  2:375520‑M1/296‑181 (MQ=255)
                                                                                                                                                                                             gAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagc                                                                                                                <  2:619491‑M1/298‑187 (MQ=255)
                                                                                                                                                                                                  cTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaa                                                                                                        >  2:126680‑M1/1‑107 (MQ=255)
                                                                                                                                                                                                   tATTCAGGACGGGTCACACGCCAAAAAAAAAGCCAGCCTGGTTCCAGACTGGCTTGTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgatgtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaag                                                                                                       <  2:380129‑M1/301‑196 (MQ=255)
                                                                                                                                                                                                    aTTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTTCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggatctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttctactcacgagcgacacagacatcacgtttttcactcttgctcccttccaggtatgccgcgggccgtataaaactttctttaaggcttagaggcaagcatcaaaaa                                                                                                        >  2:56984‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                    aTTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcacaaa                                                                                                        >  1:7223‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                    aTTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgacccacgagcgacccagacatcaagtttttaactcctgctcccttcgaggtatgccgcgcgtcgtctaaaactttctttaaggcttagcgtcaagcatcaaaaag                                                                                                       >  1:341455‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                    aTTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacctcgacctgcaggtcgacggatccccggaataacgtttcgacttacgagcgacacagacatcacgtttttcactcctgttcccttcgaggtctgccgcgtgtcgtatacaactttcttta                                                                                                                                >  2:655426‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                            cGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaa                                                                                               <  2:605118‑M1/300‑204 (MQ=255)
                                                                                                                                                                                                              ggTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAAGGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcggctcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggt                                                                                                                                                             <  2:564605‑M1/236‑142 (MQ=255)
                                                                                                                                                                                                              ggTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttcgggtcaagcatcaaaaagcggtcaaatt                                                                                             >  2:625252‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                              ggTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcactcatcaaaaagcggtcaaatt                                                                                             >  2:205170‑M1/1‑95 (MQ=255)
                                                                                                                                                                                                                   cacGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacgg                                                                                       <  1:649246‑M1/301‑212 (MQ=255)
                                                                                                                                                                                                                    acGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacgg                                                                                       <  1:126670‑M1/300‑212 (MQ=255)
                                                                                                                                                                                                                     cgcgCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatac                                                                                         <  1:571101‑M1/297‑210 (MQ=255)
                                                                                                                                                                                                                           aaaaaaaGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccc                                                                               <  2:428313‑M1/301‑220 (MQ=255)
                                                                                                                                                                                                                           aaaaaaaGACAGCTTGTTTCCAGACTGCATTTTGTGTTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacgg                                                                                       <  1:285202‑M1/293‑212 (MQ=255)
                                                                                                                                                                                                                            aaaaaaGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcc                                                                                 >  2:255086‑M1/1‑81 (MQ=255)
                                                                                                                                                                                                                                          gTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgttttgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtatacaactttcttcaaggcttagagtcaagcatcaaaaaccggtcaaattatacggtgcgcccccgcgatttaaaa                                                                   >  1:636052‑M1/1‑67 (MQ=255)
                                                                                                                                                                                                                                            ttCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaag                                                                                                                                                                                                                                                                                   >  2:267741‑M1/1‑65 (MQ=255)
                                                                                                                                                                                                                                            ttCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaag                                                                                                                                                                                                                                                                                   <  1:267731‑M1/88‑24 (MQ=255)
                                                                                                                                                                                                                                                         tttGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtgatttcaaacaataagtagccaaaagg                                                 <  2:492382‑M1/301‑250 (MQ=255)
                                                                                                                                                                                                                                                               cTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtgatttcaaacaataagtagccaaaaggtgaat                                            <  2:7223‑M1/300‑255 (MQ=255)
                                                                                                                                                                                                                                                                           gtgtTACATCGCTTTTTTGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccc                                                                               >  1:318745‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                                                           gtgtTACATCGCTTTTTTCGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccc                                                                               <  2:318756‑M1/253‑220 (MQ=255)
                                                                                                                                                                                                                                                                                       ttttttGGTCAACTCGATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtgatttcaaacaataagtagccaaaaggtgaataaatgtttaaatatttttccagtgc                   <  1:625238‑M1/301‑280 (MQ=255)
                                                                                                                                                                                                                                                                                                   acagagaacgtttaggctggagatgcttcgaagttcctatactttctagcgaatatgaacttcggactgcaggtcgacggatcccctgaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtaggccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtgatttcaaacaataagtagccaaaaggtgaataaatgtttaaatatttttccagtg                    <  1:231213‑M1/282‑279 (MQ=255)
                                                                                                                                                                                                                                                                                                      gATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtga                                                                          <  1:71633‑M1/231‑225 (MQ=255)
                                                                                                                                                                                                                                                                                                      gATAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgtga                                                                          >  2:71636‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                                                                                       aTAACGtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctcccttcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatccaaaagcggtcaaattatacggtgcgcccccgtgatttcaaacaataagtagccaaaaggtgaataaatgtttaaatatttttcca                       >  2:231223‑M1/1‑6 (MQ=255)
                                                                                                                                                                                                                                                                                                           ngtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaataacgtttcgactcacgagcgacacagacatcacgtttttcactcctgctccattcgaggtatgccgcgtgtcgtataaaactttctttaaggcttagagtcaagcatcaaaaagcggtcaaattatacggtgcgcccccgcgatttcaaacaataagaagccaaaaggtgaatacatgtttaaatatttttacagtgaataaaattgcgactttt  >  1:486151‑M1/2‑2 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        
GAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTTTTGTGCTTTTCAAGCCGGTGTTACATCGCTTTTTTGGTCAACTCGATAACGCGCAGCTGCGCGATCGCTTTGGCCAGTTCCGCAGACGCCTGAGCGTAATCTACGTCGCCGTGAGAGCTGCTAATGTGCTCTTCAGCCTTACGTTTCGCTTCCATGGCTCGCGCTTCGTCGAGATCCTGGCCGCGAATTGCGGTGTCGGCCAGAACGGTCACGTTGCCAGGCTGCACTTCAAGAATGCCGCCAGACAGATAGATAAACTCTTCGTGACCGTGCTGTTTCACGATGCGAATCATACCAGGCTTAATGGCGGTGAGCAGCGGTGCGTGGCCAGGGTAGATCCCCAGTTC  >  NC_000913/3912531‑3913127

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: