Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,095,033 C→A 8.8% S16S (TCG→TCT yggR ← putative PilT family AAA+ ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,095,0330CA 8.8% 237.9 / 2.6 92S16S (TCG→TCTyggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (27/55);  new base A (8/2);  total (35/57)
Fisher's exact test for biased strand distribution p-value = 5.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAA  >  NC_000913/3094884‑3095181
                                                                                                                                                     |                                                                                                                                                    
aGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATcc                                                                                                                                                      <  2:463506/150‑1 (MQ=255)
   tGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCgag                                                                                                                                                   <  1:1166082/150‑1 (MQ=255)
         ttCTCCAGCAATATTCCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACTTCCGCCGAATCAAACGCCGCAGCTTCCAATCTGCCGCGAATGCCCCATCGTGCCGGCCAGGCGCTCCCCAGGTGTAGATCCGAGACGtaa                                                                                                                                             >  2:537050/1‑148 (MQ=255)
           ctcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGTCTTTATg                                                                                                                                           >  2:169770/1‑150 (MQ=255)
              cAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                                                                                                                        >  1:396280/1‑150 (MQ=255)
              cAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                                                                                                                        <  1:244064/150‑1 (MQ=255)
              cAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                                                                                                                        <  1:688564/150‑1 (MQ=255)
              cAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCACAGAAGTAAAGCtt                                                                                                                                        >  2:890864/1‑150 (MQ=255)
                gCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCATTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa                                                                                                                                      <  1:233602/150‑1 (MQ=255)
                gCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa                                                                                                                                      >  1:1746961/1‑150 (MQ=255)
                 cAATATTGCGCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGAATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTac                                                                                                                                     <  2:1434596/150‑1 (MQ=255)
                 cAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTCGCTCTTCGACGGCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCTCTGCACAGGTGTAGATCCGAGACGTTATGCTTTac                                                                                                                                     <  2:1702094/150‑1 (MQ=255)
                 cAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTac                                                                                                                                     >  2:1440283/1‑150 (MQ=255)
                 cAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTac                                                                                                                                     <  1:1642558/150‑1 (MQ=255)
                   atatTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTacac                                                                                                                                   <  2:658604/150‑1 (MQ=255)
                        gCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTCCACTAAgg                                                                                                                              >  2:1439410/1‑150 (MQ=255)
                        gCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAgg                                                                                                                              <  1:666040/150‑1 (MQ=255)
                        gCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGAGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAgg                                                                                                                              <  1:1758767/150‑1 (MQ=255)
                         cccGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAggg                                                                                                                             <  1:1830784/150‑1 (MQ=255)
                         cccGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAggg                                                                                                                             <  1:1499658/150‑1 (MQ=255)
                         cccGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGAGGGCCAGGCGCTGCACCGGTGTAGCTCCGAGGCGTTATGCTTTACACTAAggg                                                                                                                             >  1:61972/1‑150 (MQ=255)
                           cGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGcc                                                                                                                           <  2:1454509/150‑1 (MQ=255)
                           cGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGcc                                                                                                                           >  1:1401178/1‑150 (MQ=255)
                           cGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGcc                                                                                                                           >  1:219070/1‑150 (MQ=255)
                           cGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGACGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGcc                                                                                                                           <  1:1317304/150‑1 (MQ=255)
                              tGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCcaca                                                                                                                        <  2:1371427/150‑1 (MQ=255)
                                   gTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttc                                                                                                                   >  1:496468/1‑150 (MQ=255)
                                   gTCATCCAGCCACTCCCGCAGAAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttc                                                                                                                   <  2:883910/150‑1 (MQ=255)
                                     cATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttctt                                                                                                                 <  2:792221/150‑1 (MQ=255)
                                     cATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttctt                                                                                                                 <  1:1574075/150‑1 (MQ=255)
                                      aTCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttcttc                                                                                                                <  2:439709/150‑1 (MQ=255)
                                       tCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCGATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTcc                                                                                                               <  2:483192/150‑1 (MQ=255)
                                       tCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTcc                                                                                                               <  2:399007/150‑1 (MQ=255)
                                       tCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTcc                                                                                                               <  2:1380093/150‑1 (MQ=255)
                                         cAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCat                                                                                                             >  1:399319/1‑150 (MQ=255)
                                         cAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATGTCTTCCat                                                                                                             >  1:1270741/1‑150 (MQ=255)
                                             cACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCACAGACGTAATCCTTTACAAAAAAAAAAACAATTTCTTCCATATTc                                                                                                         >  2:763165/1‑150 (MQ=255)
                                               cTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTACATCCGAGACGTAATGCTATACACTAAGGGCCACAATTTCTTCCATATTCAt                                                                                                       >  1:52850/1‑150 (MQ=255)
                                                tCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATa                                                                                                      <  1:1482990/150‑1 (MQ=255)
                                                tCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATa                                                                                                      <  1:1490525/150‑1 (MQ=255)
                                                  ccGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACt                                                                                                    <  1:300081/150‑1 (MQ=255)
                                                   cGCAGGAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAAAGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTa                                                                                                   <  1:781844/150‑1 (MQ=255)
                                                            tctTCGACGTCCGGCGTATCAAACGGCGCCGCTTCCATTCTGAGGCGAATGCGGCATCGTTCGGGCCAGGCGCAGCCCAGGTGTAGATCCCAGCAGTGATGCGTTACACTAAGGGCCACAAGTTCTGGCATATTCTTACTATGATCCTCg                                                                                          >  1:927763/1‑150 (MQ=255)
                                                            tctTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATGCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCg                                                                                          <  1:81534/150‑1 (MQ=255)
                                                            tctTCGACGTCCGGCGTATCAAACGGCGCAGCTGCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGGTATGCTTTACACTAAGGTCGACAATTTGTTCCATATTCATACTAAGAACCTCg                                                                                          >  1:931648/1‑150 (MQ=255)
                                                             ctTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCCAGACGTAAAGCTATACACAAATCGACACAATTTCGTCCATATTCATACTAAGATCCGCgg                                                                                         >  2:1768156/1‑150 (MQ=255)
                                                                    gTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACCCTAAGGGCCACATTTTCTTCCATATTCATACTAAGATCCTCGGAAAATGa                                                                                  >  1:1386432/1‑150 (MQ=255)
                                                                      ccGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGaat                                                                                <  1:317148/150‑2 (MQ=255)
                                                                      ccGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAAc                                                                                <  1:169479/150‑1 (MQ=255)
                                                                          ggTATCTAACGGCGCCGCGTCGATTCTGCCGCGTTTGCGCCATCGGGCGGGCCAGGCGCTGAACTGGTGTAGATCCGATAAGTTATGCTTTACAGTAAGGGCCAGAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGata                                                                            <  1:1157700/149‑1 (MQ=255)
                                                                           gTATCACACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCGAGGCGCTGCACAGGTGTAGATCCGTGACGTTATGCTTAACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGatat                                                                           <  1:1415332/150‑1 (MQ=255)
                                                                             aTCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCTTGCGGGCCAGGCGCTGCACAGGTGTAGATCCTAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCTGAAAATGAACGATATTg                                                                         <  1:1103554/150‑1 (MQ=255)
                                                                               cAAACGGCGCAGCTTCCTTTCTGCCGCGAATGAGCCATCGTGAGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGTGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTgcg                                                                       <  2:1231952/150‑1 (MQ=255)
                                                                               cAAACGGCGCAGCTTCCATTCTGCCGCGACGGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGGTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTgcg                                                                       <  2:1627185/150‑1 (MQ=255)
                                                                                 aaCGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCa                                                                     <  2:1790585/150‑1 (MQ=255)
                                                                                        gaGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCGAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTg                                                              <  2:972393/149‑1 (MQ=255)
                                                                                         aGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTgg                                                             <  1:698372/150‑1 (MQ=255)
                                                                                         aGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTgg                                                             <  1:525130/150‑1 (MQ=255)
                                                                                          gCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGc                                                            >  1:1762489/1‑150 (MQ=255)
                                                                                             tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCAGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGcaca                                                         >  1:107189/1‑150 (MQ=255)
                                                                                               cATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGCCGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCCTACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAgg                                                       <  1:1768156/150‑1 (MQ=255)
                                                                                                aTTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGt                                                      <  1:1564892/150‑1 (MQ=255)
                                                                                                aTTCTGCCGCGAATACGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCACAGACGTAAACCTATAAAAAAAAAAAAACAATTTCTTCCATATTCATACTAAGCTCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGt                                                      >  2:354866/1‑150 (MQ=255)
                                                                                                   cTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCg                                                   <  1:656393/150‑1 (MQ=255)
                                                                                                      ccGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGa                                                <  2:219070/150‑1 (MQ=255)
                                                                                                          gAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGATCGATATTGCGAATAACCTGGCACAGGTCCGGGACaaa                                            >  1:1380183/1‑150 (MQ=255)
                                                                                                             tGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAAtc                                         <  1:193060/150‑1 (MQ=255)
                                                                                                             tGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAAtc                                         <  2:1401178/150‑1 (MQ=255)
                                                                                                                 caaTCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCAGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAg                                     <  1:271322/148‑1 (MQ=255)
                                                                                                                    tCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGACCGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCgg                                  >  1:338306/1‑150 (MQ=255)
                                                                                                                      gTGCGGGCCAGGCGCTGCACAGGTGTAGATCACAGAAGTCAACCTATCAAAAAAAAAAAAAAAAATCTTCCATATTCATACTAAGATGCTCGGAAAAGGACCGATATTGCGCATTACCTGGCACAGGGCCGGGACAAAAGCTCAGCGGCt                                >  2:1174162/1‑150 (MQ=255)
                                                                                                                       tGCGGGCCAGGCGCTGCACAGGTGTAGATCACAGACGTAAACCTATAAAAAAAAAAAAAAAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTg                               >  2:1515573/1‑150 (MQ=255)
                                                                                                                          aggCCAGGCGCTGCACAGGTGTAGATCACAGACGTAAACCTATCAAAAAAAAAAAAAAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCaa                            >  2:1523742/2‑150 (MQ=255)
                                                                                                                          aggCCAGGCGCTGCACAGGTGTAGATCACAGACGTAAACCTATCAAAAAAAAAAAAAAATTTCTTCCATATTCATAATAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGaaa                            >  2:957151/2‑150 (MQ=255)
                                                                                                                            gCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCCAcg                          >  1:1851898/1‑150 (MQ=255)
                                                                                                                            gCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAAcg                          >  1:1228482/1‑150 (MQ=255)
                                                                                                                             ccAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATACCCTGGCACATGTCCGGGACAAAATCTCAGCGGCTGCAAcgc                         >  1:976028/1‑150 (MQ=255)
                                                                                                                             ccAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAAcgc                         <  2:396280/150‑1 (MQ=255)
                                                                                                                              cAGGCGCTGCACAGGTGTAGCTCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACTAAATCTCCGCGGCTGCAAcgcg                        >  1:490726/1‑150 (MQ=255)
                                                                                                                               aGGCGCTGCACTGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTTAGATCCTCGGAATATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGt                       <  2:650047/150‑1 (MQ=255)
                                                                                                                                  cgcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGc                    <  2:399319/150‑1 (MQ=255)
                                                                                                                                   gcgtcacaGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGAGCCTCGGAAACTGACCGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCg                   <  1:82944/146‑1 (MQ=255)
                                                                                                                                   gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCg                   <  2:496468/150‑1 (MQ=255)
                                                                                                                                   gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCg                   >  1:300490/1‑150 (MQ=255)
                                                                                                                                      gCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCTGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGcc                <  2:1232344/150‑1 (MQ=255)
                                                                                                                                         caGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGtt             <  1:1575931/150‑1 (MQ=255)
                                                                                                                                            gtgtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTtctc          <  1:514739/150‑1 (MQ=255)
                                                                                                                                            gtgtAGATCAGATACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCCGCAACGCGTTGCAGCCgctgag          <  1:242394/150‑6 (MQ=255)
                                                                                                                                             tgtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTcc         <  2:1281755/150‑1 (MQ=255)
                                                                                                                                             tgtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTcc         >  1:932357/1‑150 (MQ=255)
                                                                                                                                              gtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATAATGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCa        <  1:169770/150‑1 (MQ=255)
                                                                                                                                                   tCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCagaaga   <  2:1520802/150‑1 (MQ=255)
                                                                                                                                                   tCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCagaaga   >  1:1424787/1‑150 (MQ=255)
                                                                                                                                                    ccGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAgaagaa  >  1:108080/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACTAAGATCCTCGGAAAATGAACGATATTGCGCATAACCTGGCACAGGTCCGGGACAAAATCTCAGCGGCTGCAACGCGTTGCGGCCGTTCTCCAGAAGAA  >  NC_000913/3094884‑3095181

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: