Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 443,023 T→C 8.7% intergenic (+26/+28) thiI → / ← yajL tRNA uridine 4‑sulfurtransferase/protein/nucleic acid deglycase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,0230TC 8.7% 196.4 / 8.3 69intergenic (+26/+28)thiI/yajLtRNA uridine 4‑sulfurtransferase/protein/nucleic acid deglycase 3
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (4/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

TTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCC  >  NC_000913/442876‑443172
                                                                                                                                                   |                                                                                                                                                     
ttGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                                                                     <  1:417754/150‑1 (MQ=255)
  ggCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTTTGGTTTTTACGTCGCATCTGTTCata                                                                                                                                                   <  2:303228/150‑3 (MQ=255)
   gCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                  >  2:569161/1‑150 (MQ=255)
     gATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                                                                                >  1:1193297/1‑150 (MQ=255)
      aTCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCg                                                                                                                                               >  2:652082/1‑150 (MQ=255)
        ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                             <  1:652082/150‑1 (MQ=255)
         tcGACCAGAACAAACCCTGGCTGCTGTGGTGGGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                                                                            <  1:395218/150‑1 (MQ=255)
          cGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGt                                                                                                                                           <  1:1304945/150‑1 (MQ=255)
               aGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                                                                      <  2:1022761/150‑1 (MQ=255)
               aGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                                                                      <  1:478352/150‑1 (MQ=255)
                  aCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTCGCCGCa                                                                                                                                   >  1:554569/1‑150 (MQ=255)
                          tGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacac                                                                                                                           >  1:903201/1‑150 (MQ=255)
                              tgctgTGGTGTGAGTGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactact                                                                                                                       >  2:916433/1‑150 (MQ=255)
                              tgctgTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactact                                                                                                                       >  2:226892/1‑150 (MQ=255)
                                     gtgtGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactac                                                                                                                        >  1:419780/1‑142 (MQ=255)
                                     gtgtGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactac                                                                                                                        <  2:419780/142‑1 (MQ=255)
                                       gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                                                                              <  1:1358124/150‑1 (MQ=255)
                                           gcgcGGGGTGATGAGCCGTCTGCATGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACTTTTGCCGCATCCGACACTACTCGTAATAATtata                                                                                                          <  2:1193297/150‑1 (MQ=255)
                                           gcgcGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                                                                          <  1:569161/150‑1 (MQ=255)
                                              cgGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAAt                                                                                                       >  1:1233604/1‑150 (MQ=255)
                                                      tgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccg                                                                                               >  1:595523/1‑150 (MQ=255)
                                                         gCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCa                                                                                            <  2:101633/150‑1 (MQ=255)
                                                            gTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                         <  1:876203/150‑1 (MQ=255)
                                                            gTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                         >  2:331040/1‑150 (MQ=255)
                                                              cTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAcc                                                                                       <  1:226892/150‑1 (MQ=255)
                                                                 cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGt                                                                                    >  2:729214/1‑150 (MQ=255)
                                                                 cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGAGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGt                                                                                    >  1:439822/1‑150 (MQ=255)
                                                                   ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGTAGGTATATCGCCCGTAATTTGTGGTTTTTAGGTCGCATCTGGTCAGATGCGACGTTTGCCACATCCGAGACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                                                                                  <  2:668570/150‑1 (MQ=255)
                                                                     cgcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgtg                                                                                <  2:1011302/150‑1 (MQ=255)
                                                                      gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                               >  2:1010409/1‑150 (MQ=255)
                                                                          tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                           <  1:916433/150‑1 (MQ=255)
                                                                           aTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCa                                                                          <  2:793427/150‑1 (MQ=255)
                                                                               gcgcGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTc                                                                      >  2:847095/1‑150 (MQ=255)
                                                                                  cgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                                                                   <  2:502659/150‑1 (MQ=255)
                                                                                      cAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt                                                               >  1:551209/1‑150 (MQ=255)
                                                                                            tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCa                                                         >  2:1148991/1‑150 (MQ=255)
                                                                                                aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGcc                                                     <  1:1010409/150‑1 (MQ=255)
                                                                                                   aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                 <  1:706900/119‑1 (MQ=255)
                                                                                                   aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                 >  2:706900/1‑119 (MQ=255)
                                                                                                       gAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaaccaa                                              >  1:1278886/1‑150 (MQ=255)
                                                                                                         aGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACa                                            <  1:847095/150‑1 (MQ=255)
                                                                                                            tataTCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                         >  1:1347693/1‑150 (MQ=255)
                                                                                                               atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                      >  1:221362/1‑150 (MQ=255)
                                                                                                                tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGAt                                     <  2:1386791/150‑1 (MQ=255)
                                                                                                                 cGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATa                                    >  1:1019306/1‑150 (MQ=255)
                                                                                                                  gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATaa                                   >  1:72956/1‑150 (MQ=255)
                                                                                                                   cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAt                                  >  2:870172/1‑150 (MQ=255)
                                                                                                                   cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCATGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAt                                  >  1:975102/1‑150 (MQ=255)
                                                                                                                       tAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTc                              <  2:1397039/150‑1 (MQ=255)
                                                                                                                       tAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTc                              <  2:595523/150‑1 (MQ=255)
                                                                                                                        aaTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCTCATCCGACACTACTCGTAATAATTATAAATCCCTCCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                             >  2:475209/1‑150 (MQ=255)
                                                                                                                        aaTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                             <  1:729214/150‑1 (MQ=255)
                                                                                                                              tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaa                       >  1:65640/1‑150 (MQ=255)
                                                                                                                                cggcaaaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAg                     >  1:495038/7‑150 (MQ=255)
                                                                                                                                 ggcaaaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGt                    >  1:758533/6‑150 (MQ=255)
                                                                                                                                  ttttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTc                   <  1:146497/150‑1 (MQ=255)
                                                                                                                                     tACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGAt                <  1:870172/150‑1 (MQ=255)
                                                                                                                                     tACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGAt                <  2:598585/150‑1 (MQ=255)
                                                                                                                                      aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATa               >  2:516460/1‑150 (MQ=255)
                                                                                                                                       cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGAcc                         <  1:1347176/139‑1 (MQ=255)
                                                                                                                                       cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGAcc                         >  2:1347176/1‑139 (MQ=255)
                                                                                                                                        gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc             <  2:1233604/150‑1 (MQ=255)
                                                                                                                                            cATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCCCCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTa         >  1:1004398/1‑150 (MQ=255)
                                                                                                                                            cATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTa         >  1:758121/1‑150 (MQ=255)
                                                                                                                                               cTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTAcc       <  2:1190254/149‑1 (MQ=255)
                                                                                                                                               cTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTAcc       >  1:1190254/1‑149 (MQ=255)
                                                                                                                                                tGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCg     >  1:726963/1‑150 (MQ=255)
                                                                                                                                                 ggTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg    >  2:678023/1‑150 (MQ=255)
                                                                                                                                                  gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGc   >  2:726733/1‑150 (MQ=255)
                                                                                                                                                   tCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGcc  <  2:207127/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCC  >  NC_000913/442876‑443172

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