Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 443,022 G→A 9.5% intergenic (+25/+29) thiI → / ← yajL tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,0220GA 9.5% 210.8 / 8.6 74intergenic (+25/+29)thiI/yajLtRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (3/4);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGG  >  NC_000913/442877‑443170
                                                                                                                                                 |                                                                                                                                                    
tGGCGATCACGACCAGAACAAAACCTCGCTGCTGTGGTGTGAGCGCGCGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCGTTAACAATGTGAAGGTATCTCGCCCGTAATTTGTGGTTTTTAAGTCGCATCTGGTCac                                                                                                                                                  >  1:163397/1‑149 (MQ=255)
   cGATCTCGACCAGAACAAAACCTGGCTGCGGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATg                                                                                                                                               >  1:229431/1‑150 (MQ=255)
    gATCTCGACCAGAACAAAACCTGGCAGCTGTGGTGTGAGCGCGGGGTGATGACCCGTCTGCAGCCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATAACGCCCGTAATTTGTGGTTTTTACTACGCATCTAGTCAGATGc                                                                                                                                              >  1:101169/1‑150 (MQ=255)
       ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                           <  2:947700/150‑1 (MQ=255)
       ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGAc                                                                                                                                           <  1:186078/150‑1 (MQ=255)
        tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCTAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTCTACGTCGCATCTGGTCAGATGCGACg                                                                                                                                          <  2:1317343/150‑1 (MQ=255)
        tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                                                                          <  1:746251/150‑1 (MQ=255)
        tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                                                                          <  2:1261096/150‑1 (MQ=255)
        tcGACCAGAACAAAACCTGGCTACTGTGGTGTGAGCGCGGGGTGATGAGCCGGCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                                                                          <  1:142683/150‑1 (MQ=255)
            ccAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTg                                                                                                                                      >  2:1033708/1‑150 (MQ=255)
                       ccTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGac                                                                                                                           <  2:696390/150‑1 (MQ=255)
                        cTGGCAGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacc                                                                                                                          >  1:355180/1‑149 (MQ=255)
                         tGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacac                                                                                                                         >  2:1384771/1‑150 (MQ=255)
                         gggCTGCTGGGGTGTGAGCGCGGGGGGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacac                                                                                                                         <  1:1507373/149‑1 (MQ=255)
                             tgctgTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGAAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactact                                                                                                                     <  2:1657685/150‑1 (MQ=255)
                                 gtggtgTGAGCGCGGGGTGATGAGCAGTCTGCAGGCGATCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCAGACACTACTCGta                                                                                                                 >  2:377409/1‑150 (MQ=255)
                                     tgtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATAACGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGATTGCCGCATACGACACAACTCGtaataa                                                                                                             >  2:1592790/1‑150 (MQ=255)
                                     tgtgAGCGCGGGGTGATGAGCCGACTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCAGACACTACTCCtaataa                                                                                                             >  2:482909/1‑150 (MQ=255)
                                      gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                                                                            >  1:18686/1‑150 (MQ=255)
                                      gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                                                                            >  1:960466/1‑150 (MQ=255)
                                        gAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCTTCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATta                                                                                                          <  1:1270890/150‑1 (MQ=255)
                                          gcgcGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                                                                        <  2:1340376/150‑1 (MQ=255)
                                          gcgcGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                                                                        <  1:9169/150‑1 (MQ=255)
                                                   gatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgc                                                                                               >  1:512610/1‑150 (MQ=255)
                                                     tgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccg                                                                                             <  1:19909/150‑1 (MQ=255)
                                                      gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccgc                                                                                            <  2:1601827/150‑1 (MQ=255)
                                                        gCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGACa                                                                                          >  1:483482/1‑150 (MQ=255)
                                                           gTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                       >  1:22413/1‑150 (MQ=255)
                                                               gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAg                                                                                   >  1:940953/1‑150 (MQ=255)
                                                                  ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGTTGCTACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                                                                                <  1:1304219/150‑1 (MQ=255)
                                                                     gcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAACGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCTTCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                             <  2:88447/150‑1 (MQ=255)
                                                                       tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATg                                                                           <  1:107854/150‑1 (MQ=255)
                                                                         tATCTGCGCGTGCAGGGCTTTAACATTGTGAAGTTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGTTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                         <  1:607643/150‑1 (MQ=255)
                                                                         tATCTGCGCGAGCAGGGCTTTCACAAGGTGACGGGATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCCGATGCGACGTTTGCCGCATCCGACACTACTCGTAATACTTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                         <  1:847862/150‑1 (MQ=255)
                                                                         tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                         <  1:25947/150‑1 (MQ=255)
                                                                            cTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACt                                                                      <  1:441583/150‑1 (MQ=255)
                                                                            cTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACt                                                                      <  2:116930/150‑1 (MQ=255)
                                                                            cTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACt                                                                      >  1:208964/1‑150 (MQ=255)
                                                                                   aGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg                                                               <  2:940953/150‑1 (MQ=255)
                                                                                   aGCAGGGCTTTAACAATGTGAAGGAATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTTGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg                                                               <  2:1191976/150‑1 (MQ=255)
                                                                                     cAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt                                                             <  1:807049/150‑1 (MQ=255)
                                                                                      aGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtt                                                            <  1:23081/150‑1 (MQ=255)
                                                                                       gggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACAATACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttt                                                           >  1:255152/1‑150 (MQ=255)
                                                                                       gggCTTTAACAATGTGAAGGTATATCGCCAGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTATTAATTTTAATTCCCTGCCGCCATCACCAGTTGTCATGCCACTTAATGGGCttt                                                           >  1:702124/1‑150 (MQ=255)
                                                                                         gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                         <  1:353068/150‑1 (MQ=255)
                                                                                         gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                         >  2:1454334/1‑150 (MQ=255)
                                                                                          cTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                        <  2:1089073/150‑1 (MQ=255)
                                                                                           tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCa                                                       >  2:971104/1‑150 (MQ=255)
                                                                                             tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACg                                                     >  1:1217702/1‑150 (MQ=255)
                                                                                                cAATGTGAAGGTATATCGCCCGTAATTTGTGGGTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGccc                                                  <  2:163397/150‑1 (MQ=255)
                                                                                                    gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaacc                                              >  1:8530/1‑150 (MQ=255)
                                                                                                      gAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaaccaa                                            >  1:1703719/1‑150 (MQ=255)
                                                                                                       aaGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCcaaccaac                                           <  1:1033708/150‑1 (MQ=255)
                                                                                                          gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAgg                                        <  1:377409/150‑1 (MQ=255)
                                                                                                           tataTCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                       >  1:792590/1‑150 (MQ=255)
                                                                                                            atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCACTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACATGTc                                      >  1:1132825/1‑150 (MQ=255)
                                                                                                             tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCg                                     <  2:22413/150‑1 (MQ=255)
                                                                                                              atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                    >  2:1663383/1‑150 (MQ=255)
                                                                                                                 gCCCGTTATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATaa                                 <  1:1067454/150‑1 (MQ=255)
                                                                                                                  cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAt                                >  2:756822/1‑150 (MQ=255)
                                                                                                                   ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGCCGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAtt                               >  2:1137007/1‑150 (MQ=255)
                                                                                                                      tAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTc                            <  2:521740/150‑1 (MQ=255)
                                                                                                                       attttGTGGTTTTTACGTCGCATCTTGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATATTTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                           <  1:1592790/148‑1 (MQ=255)
                                                                                                                       aaTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                           <  2:18686/150‑1 (MQ=255)
                                                                                                                          ttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCACTGCCGCCATCCCCAGTTGTGATGCCACTTCATGGGCTTTTTAACGCCCAACCAACAGGTCGATAATTTTCAGAc                        >  1:165238/1‑150 (MQ=255)
                                                                                                                             tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaa                     >  1:1707881/1‑150 (MQ=255)
                                                                                                                              ggTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                    >  1:799463/1‑150 (MQ=255)
                                                                                                                                 ttttACGTCGCATCTGGTCAGATGCGACGGTTGCCGCATCCGACACGACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGGGATGCCACTTCATGGGCTTTTTCACGCCCAACCAGCAGGGCGATAATTTTCAGGCCAAAGTc                 >  1:607282/1‑150 (MQ=255)
                                                                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg            >  2:1026614/1‑150 (MQ=255)
                                                                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg            >  1:933737/1‑150 (MQ=255)
                                                                                                                                       gTCGCATCTGGTCCGTTGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc           <  2:868261/150‑1 (MQ=255)
                                                                                                                                       gTCGCATCTGACCAGATGCGACGTTTGCCGCATACGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGGTTTATCACGCACAACCAACAGGTCGATAATTTTCAGACTAAAGTCGATAGc           >  1:1017823/1‑150 (MQ=255)
                                                                                                                                            aTCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTAc      >  2:746505/1‑150 (MQ=255)
                                                                                                                                             tCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATAACCAGTTGTGATGCCACTTCATGGGCATTTTCACCCCCAACCAACAGGTAGATAATTTTCAGACCAAAGTCGATAGCTGTAcc     >  2:147446/1‑150 (MQ=255)
                                                                                                                                                ggTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCATAGTCGATAGCTGTACCCgg  >  1:19912/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGG  >  NC_000913/442877‑443170

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