Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,957,484 T→C T217A (ACC→GCC)  torY ← cytochrome c quinol dehydrogenase TorY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,957,4840TC100.0% 396.4 / NA 109T217A (ACC→GCC) torYcytochrome c quinol dehydrogenase TorY
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (55/54);  total (55/54)

Only 100 of 111 total aligned reads displayed.

TCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGTTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGC  >  NC_000913/1957351‑1957617
                                                                                                                                     |                                                                                                                                     
tcttcACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAgg                                                                                                                                 >  1:32691/1‑140 (MQ=255)
tcttcACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAgg                                                                                                                                 >  1:1223837/1‑140 (MQ=255)
  ttcACTTAATGTGGCTAGCACCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGAt                                                                                                                               >  1:1177170/1‑140 (MQ=255)
    cACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTc                                                                                                                             >  2:460435/1‑140 (MQ=255)
     aCTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCTCTGGCATGAACGACGGTGAGAGCGGCTCCAGGATTCa                                                                                                                            >  1:607128/1‑140 (MQ=255)
       ttAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACg                                                                                                                          <  2:1512368/140‑1 (MQ=255)
         aaTGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGTCATCAACGACGCTGAGATCGGCTTCAGGATTCACGGt                                                                                                                        >  1:821978/1‑140 (MQ=255)
           tgtgGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAg                                                                                                                      >  1:1401558/1‑140 (MQ=255)
            gtgGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGc                                                                                                                     <  2:96999/140‑1 (MQ=255)
            gtgGCTAGCGCCAGCCGTTGACATGCCGCCACGTAGAGCCTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGc                                                                                                                     <  2:586749/140‑1 (MQ=255)
               gCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGAGCGGCTCCAGTATTCACGGTAGCCAg                                                                                                                  >  1:485865/1‑140 (MQ=255)
                 tAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGTTGAGATGGGCTCCAGGATTCATTGTAgtctgcc                                                                                                                >  2:1766472/1‑134 (MQ=255)
                   gcgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTa                                                                                                                       <  1:1823936/131‑1 (MQ=255)
                   gcgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTa                                                                                                                       >  2:1823936/1‑131 (MQ=255)
                   gcgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCg                                                                                                              <  1:82350/140‑1 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACTGTGAGATAGGCACCAGGATTCCCGGTAGCCAGCTCGc                                                                                                             >  1:548756/1‑140 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCa                                                                                                                            >  1:1220419/1‑125 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCa                                                                                                                            <  2:1220419/125‑1 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGc                                                                                                             >  1:1700537/1‑140 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGc                                                                                                             >  1:1436994/1‑140 (MQ=255)
                    cgcCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGc                                                                                                             >  1:981030/1‑140 (MQ=255)
                      ccAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCt                                                                                                           <  1:460435/140‑1 (MQ=255)
                         gCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTAtg                                                                                                        >  1:1558134/1‑140 (MQ=255)
                         gCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTAtg                                                                                                        >  1:1419645/1‑140 (MQ=255)
                         gCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTAtg                                                                                                        >  2:2427/1‑140 (MQ=255)
                          ccGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTAtgt                                                                                                       >  2:311536/1‑140 (MQ=255)
                             ttGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTgag                                                                                                    >  2:1396484/1‑140 (MQ=255)
                              tGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTgaga                                                                                                   <  2:1030166/140‑1 (MQ=255)
                              tGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTgaga                                                                                                   <  2:1126122/140‑1 (MQ=255)
                              tGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTgaga                                                                                                   >  1:1329264/1‑140 (MQ=255)
                               gACCTGCCGCCAGGTAGTGCGTTTTTTCACTCCCCCGCAGTTGATAACACTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGaa                                                                                                  <  2:548756/140‑1 (MQ=255)
                                     ccgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTg                                                                                            >  1:1414716/1‑140 (MQ=255)
                                     ccgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:1690997/140‑1 (MQ=255)
                                      cgccAGGTGGAGCGTGGCTTCATTGCCCTGCATGTGATACCCCTGCAAGAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:894056/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:1436994/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:1401558/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:33724/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:1579519/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:84779/140‑1 (MQ=255)
                                      cgccAGGTAGAGCGTGTATTCAATGCCCTGCTTTTTATAGCCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGa                                                                                           <  2:485865/140‑1 (MQ=255)
                                        ccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGaaa                                                                                         <  2:218876/140‑1 (MQ=255)
                                        ccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGaaa                                                                                         <  1:1022711/140‑1 (MQ=255)
                                        ccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGaaa                                                                                         <  2:1403060/140‑1 (MQ=255)
                                        ccAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGaaa                                                                                         <  1:105704/140‑1 (MQ=255)
                                          aGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAgg                                                                                       <  2:981030/140‑1 (MQ=255)
                                               gagCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATaaa                                                                                  >  1:1822893/1‑140 (MQ=255)
                                                agCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCTCCTATGAGAGAAGTTTTGAAAGGATAAat                                                                                 >  1:634452/1‑140 (MQ=255)
                                                   gtgtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGtttt                                                                                             >  2:1259962/1‑125 (MQ=255)
                                                   gtgtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGtttt                                                                                             <  1:1259962/125‑1 (MQ=255)
                                                    tgtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGAt                                                                                     <  2:1577775/132‑1 (MQ=255)
                                                    tgtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGAt                                                                                     >  1:1577775/1‑132 (MQ=255)
                                                    tgtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATAtgt                                                                             <  1:2427/140‑1 (MQ=255)
                                                     gtTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAATGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCTAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATAtgtg                                                                            <  2:1197869/140‑1 (MQ=255)
                                                      ttttCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGaaa                                                                                         <  1:1596491/126‑1 (MQ=255)
                                                      ttttCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTgc                                                                           >  1:1826466/1‑140 (MQ=255)
                                                       tttCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTgcg                                                                          <  2:409019/140‑1 (MQ=255)
                                                       tttCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTgcg                                                                          <  2:1177170/140‑1 (MQ=255)
                                                        ttCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTgcgc                                                                         <  1:261193/140‑1 (MQ=255)
                                                          cACTGCCCTGCATTTGAGAACCCTGCAACAGAACGCTCGTCTGTTTGCTACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCa                                                                       <  2:607128/140‑1 (MQ=255)
                                                           aCTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAt                                                                      >  2:1744223/1‑140 (MQ=255)
                                                            cTGCCCTGCATTTGATAACCCTGCATCAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAtt                                                                     <  1:742574/140‑1 (MQ=255)
                                                            cTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAtt                                                                     <  2:1700537/140‑1 (MQ=255)
                                                               cccTGCATTTGATAACCCTGCAACAGAACGATCGGTTGTCTGCCACTGACATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGt                                                                  <  2:634452/140‑1 (MQ=255)
                                                                ccTGCATTTGATAACCCCGCAACAGAACGATCGGCCGTTCGCCACTGGAATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTa                                                                 <  2:531332/140‑1 (MQ=255)
                                                                  tGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGTCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAtt                                                                     >  1:461575/1‑134 (MQ=255)
                                                                  tGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGTCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAtt                                                                     <  2:461575/134‑1 (MQ=255)
                                                                   gCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCAtt                                                                     >  2:1686535/1‑133 (MQ=255)
                                                                    cATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTg                                                             >  2:148529/1‑140 (MQ=255)
                                                                    cATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTg                                                             <  2:1223837/140‑1 (MQ=255)
                                                                         gATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTa                                                        >  2:1802044/1‑140 (MQ=255)
                                                                         gATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTa                                                        >  1:421598/1‑140 (MQ=255)
                                                                         gATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTa                                                        >  2:319016/1‑140 (MQ=255)
                                                                             aCCCTGCATCAGAACGTTCGGCTGTTTGCCTTTCGCATCCACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCgg                                                    <  2:901362/140‑1 (MQ=255)
                                                                              cccTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGc                                                   >  2:1649500/1‑140 (MQ=255)
                                                                              cccTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGc                                                   >  1:1088651/1‑140 (MQ=255)
                                                                                cTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCt                                                 <  2:1414716/140‑1 (MQ=255)
                                                                                 tGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTg                                                >  1:1034626/1‑140 (MQ=255)
                                                                                 tGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTg                                                >  1:1688008/1‑140 (MQ=255)
                                                                                 tGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTg                                                >  2:960987/1‑140 (MQ=255)
                                                                                 tGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTg                                                >  1:958514/1‑140 (MQ=255)
                                                                                    aaCAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCCCCATTACTAACTGAAGTAGCGGCCTGACt                                             >  1:638631/1‑140 (MQ=255)
                                                                                        gAACGATCGGCTGTTTGTCACTCGCATCAACAACGCTGACATCGGCTCAAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTcc                                         <  2:603768/140‑1 (MQ=255)
                                                                                        gAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTcc                                         <  1:319016/140‑1 (MQ=255)
                                                                                         aaCGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAACCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCa                                        >  2:1264810/1‑140 (MQ=255)
                                                                                             aTCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTc                                    >  2:100086/1‑140 (MQ=255)
                                                                                             aTCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTc                                    >  1:1625884/1‑140 (MQ=255)
                                                                                              tCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCa                                   >  1:441643/1‑140 (MQ=255)
                                                                                                     tttGCCACTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGcgccg                            <  2:32691/140‑1 (MQ=255)
                                                                                                           aCTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGt                      <  2:1822893/140‑1 (MQ=255)
                                                                                                           aCTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGt                      <  1:1826298/140‑1 (MQ=255)
                                                                                                            cTGGCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTc                     <  1:1397521/140‑1 (MQ=255)
                                                                                                              ggCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCAt                   <  2:958514/140‑1 (MQ=255)
                                                                                                               gCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGTCTTTCCAGCTCATGCGCCGCGTTGTCATc                  >  1:599807/1‑140 (MQ=255)
                                                                                                               gCATCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTACCGGCCTGACTTGCCAGCTCATACGCCGCGTTGTCATc                  >  1:782496/1‑140 (MQ=255)
                                                                                                                 aTCAACGACGGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCa                >  1:1045940/1‑140 (MQ=255)
                                                                                                                         cgGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  2:301982/1‑140 (MQ=255)
                                                                                                                         cgGTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  1:1797129/1‑140 (MQ=255)
                                                                                                                          ggTGAGATCGGCTCCAGGATTCCCGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAttt       <  2:817539/140‑1 (MQ=255)
                                                                                                                          ggTGAGATCGGCTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAttt       <  2:421598/140‑1 (MQ=255)
                                                                                                                               gaTCGGCTTCAGGATTCACGGTGGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGc  <  2:821978/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTGTTTTCACTGCCCTGCATTTGATAACCCTGCAACAGAACGATCGGCTGTTTGCCACTGGCATCAACGACGGTGAGATCGGTTCCAGGATTCACGGTAGCCAGCTCGCCTATGTGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTAACTGAAGTAGCGGCCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGC  >  NC_000913/1957351‑1957617

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