Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 900,202 T→A 6.2% K125I (AAA→ATA) ‡ artJ ← arginine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913900,2020TA 6.2% 269.7 / 5.4 95K125I (AAA→ATA) ‡artJarginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (49/40);  new base A (4/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCA  >  NC_000913/900054‑900349
                                                                                                                                                    |                                                                                                                                                   
gtCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCttt                                                                                                                                                    <  2:710652/150‑1 (MQ=255)
gtCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCttt                                                                                                                                                    >  2:601896/1‑150 (MQ=255)
 tCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTc                                                                                                                                                   <  1:1178117/150‑1 (MQ=255)
  cACCAAGTACCCCATCAATAAGACCTTTTTTCAGATCGTTAAAGGCATTCTGATAACTGTCATAAGTTACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCa                                                                                                                                                  <  1:1597831/150‑1 (MQ=255)
   aCCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAg                                                                                                                                                 <  2:336970/150‑1 (MQ=255)
   aCCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAg                                                                                                                                                 >  2:309641/1‑150 (MQ=255)
    ccAAATACCCCATCAAAACGACCATTTTTCAGATCGAAAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGCGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTATACATACCAATACGTTTGACTTTCaaa                                                                                                                                                >  2:189752/1‑148 (MQ=255)
     cAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACATCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGAt                                                                                                                                               >  2:173231/1‑150 (MQ=255)
      aaaTTCACCATCAATAAGACCATTTTTCATATCGATAAAGGCATTCTTATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTTCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATc                                                                                                                                              <  1:732768/150‑1 (MQ=255)
       aaTACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCg                                                                                                                                             <  2:1534993/150‑1 (MQ=255)
         tagaccATCAATACGACCTTTATTGAGATCGTTAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCAGTGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGc                                                                                                                                           <  1:25809/146‑1 (MQ=255)
         tACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGc                                                                                                                                           <  1:671924/150‑1 (MQ=255)
         tACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGc                                                                                                                                           <  1:1295854/150‑1 (MQ=255)
              cATCAATACGACCATTTTTCATATCGATTAAGGCATTCTGATAACTGTCTTAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACg                                                                                                                                      <  2:1048573/150‑1 (MQ=255)
                 cAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGtat                                                                                                                                   <  1:1628984/150‑3 (MQ=255)
                 aaaTACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGtat                                                                                                                                   <  2:1296063/149‑3 (MQ=255)
                  aaTACGACCATTTTTCAGATCGATAAAGGCATTATGATAACTGTCATAAGAGACAGTTTTCACTTCCGTGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGtttt                                                                                                                                  <  1:88294/150‑1 (MQ=255)
                   aTACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTg                                                                                                                                 >  2:1044594/1‑150 (MQ=255)
                      cGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAg                                                                                                                              >  1:1608487/1‑150 (MQ=255)
                      cGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAg                                                                                                                              >  1:816949/1‑150 (MQ=255)
                             ttttCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtt                                                                                                                       <  1:601321/150‑1 (MQ=255)
                             ttttCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtt                                                                                                                       <  1:1140544/150‑1 (MQ=255)
                             ttttCAGATCGAAAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGACTTTCAGATCTTCAAACTTTTTGTAGGTATCtt                                                                                                                       >  1:268787/1‑150 (MQ=255)
                                 cAGATCGATAAAGGCATTCTGATAACTGTCATAAGATACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtttttt                                                                                                                   <  2:149984/150‑1 (MQ=255)
                                  aGATCGATAACGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTAAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTg                                                                                                                  <  2:382220/150‑1 (MQ=255)
                                  aGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTg                                                                                                                  <  1:901185/150‑1 (MQ=255)
                                  aGATCGATAAAGGCATTCTGATAACTGTCATAAGAGAAAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGCTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTg                                                                                                                  <  2:941904/150‑1 (MQ=255)
                                   gATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTgg                                                                                                                 >  2:1731060/1‑150 (MQ=255)
                                     tCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCa                                                                                                               >  1:220426/1‑150 (MQ=255)
                                       gATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTTGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTATTTGGCAAt                                                                                                             >  1:332763/1‑150 (MQ=255)
                                          aaaGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGTCAATCac                                                                                                          >  2:1517950/1‑150 (MQ=255)
                                           aaGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCacg                                                                                                         >  1:236264/1‑150 (MQ=255)
                                               cATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACgg                                                                                                     >  2:446747/1‑150 (MQ=255)
                                               cATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACgg                                                                                                     >  1:1428470/1‑150 (MQ=255)
                                                aTTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGc                                                                                                    >  2:964824/1‑150 (MQ=255)
                                                   cTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACCTTTGCCTTACAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGa                                                                                                 >  2:978710/1‑150 (MQ=255)
                                                     gATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGt                                                                                               >  2:1311687/1‑150 (MQ=255)
                                                       tAACTGTTATATGAGTCTGTTTTCACTTCCGGGTGCTGATCCTTAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGttt                                                                                             <  1:570752/150‑1 (MQ=255)
                                                        aaCTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGtttt                                                                                            <  1:309641/150‑1 (MQ=255)
                                                          cTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCa                                                                                          >  1:595655/1‑150 (MQ=255)
                                                           tGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAt                                                                                         <  2:352220/150‑1 (MQ=255)
                                                               aTAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGAAGGCTGAGTTTTCAtagta                                                                                     >  2:809842/1‑150 (MQ=255)
                                                                tAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtagtag                                                                                    >  2:763897/1‑150 (MQ=255)
                                                                 aaGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAgg                                                                                   >  2:406554/1‑150 (MQ=255)
                                                                  agagACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                  >  2:420524/1‑150 (MQ=255)
                                                                  agagACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                  <  1:886596/150‑1 (MQ=255)
                                                                  agagACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                  >  1:610757/1‑150 (MQ=255)
                                                                       cAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggcggtgt                                                                             >  1:693139/1‑144 (MQ=255)
                                                                         gTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTgg                                                                             >  1:206293/1‑148 (MQ=255)
                                                                         gTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTgg                                                                             <  2:206293/148‑1 (MQ=255)
                                                                         gTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCATTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTa                                                                           >  1:1239594/1‑150 (MQ=255)
                                                                            ttCACTTCCGGGTGCTGATCCTGAATATATATCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAAc                                                                        <  1:1565863/150‑1 (MQ=255)
                                                                               aCTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCATACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTATTCATAGTAGGGCGTGGTAAACGAt                                                                     >  1:1175404/1‑150 (MQ=255)
                                                                                  tCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATAcc                                                                  <  1:580244/150‑1 (MQ=255)
                                                                                  tCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGAACCGTTTTCCATCCCAATACGTTTTCCTTTCAGATCGGCACACGATTTGTAGGTATCTGTTTTGGCAATCACGACGGCTCAGTTTACATAGTAGAGCGTGTTAAACGATAcc                                                                  >  2:1089422/1‑150 (MQ=255)
                                                                                     gggTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATTCCTGt                                                               >  2:787939/1‑150 (MQ=255)
                                                                                      ggTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGtt                                                              >  2:548652/1‑150 (MQ=255)
                                                                                       gTGCTGATCCTGAATATATTTCTTGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGttt                                                             <  2:332763/150‑1 (MQ=255)
                                                                                        tGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTg                                                            <  2:220426/150‑1 (MQ=255)
                                                                                        tGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTg                                                            >  1:446482/1‑150 (MQ=255)
                                                                                              tCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGc                                                      <  1:326548/150‑1 (MQ=255)
                                                                                                 tGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTATTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGAGTGGTAAACGATACCTGTTTGCTACGCTca                                                   >  1:1422007/1‑149 (MQ=255)
                                                                                                     tatatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGTCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGTGCGTGGTAAACGATACCTGTTTGCTACGCTCCgggg                                               <  2:986258/150‑1 (MQ=255)
                                                                                                       tatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTg                                             <  2:1239271/150‑1 (MQ=255)
                                                                                                       tatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTg                                             <  2:700868/150‑1 (MQ=255)
                                                                                                       aatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGAGCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTg                                             <  2:186803/149‑1 (MQ=255)
                                                                                                          ttCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGata                                          <  2:1024020/150‑1 (MQ=255)
                                                                                                            gtGGTGCGTCGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGAAGGTATCTTTTTTGCCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATc                                        <  2:1117776/149‑1 (MQ=255)
                                                                                                            cTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATc                                        >  2:1468508/1‑150 (MQ=255)
                                                                                                               gTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCAt                                     >  2:384354/1‑150 (MQ=255)
                                                                                                                 gCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATAc                                   >  2:591957/1‑150 (MQ=255)
                                                                                                                  cGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTTATATCCATAcc                                  <  1:173231/150‑1 (MQ=255)
                                                                                                                   gtagtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCg                                 >  1:1134907/1‑150 (MQ=255)
                                                                                                                    tagtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCgg                                >  1:1474441/1‑150 (MQ=255)
                                                                                                                      gtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGgag                              <  2:1264790/150‑1 (MQ=255)
                                                                                                                      gtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTCCGCTCCGGGGTGATATCCATACCGgaga                             >  1:875745/1‑150 (MQ=255)
                                                                                                                       taCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACTATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGgaga                             >  2:561/1‑150 (MQ=255)
                                                                                                                       taCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGgaga                             >  2:1415369/1‑150 (MQ=255)
                                                                                                                       taCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCAGGGGTGATATCCATACCGgaga                             >  1:576136/1‑150 (MQ=255)
                                                                                                                        aCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGCTACCGGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGAt                            >  2:393055/1‑150 (MQ=255)
                                                                                                                           gTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCGGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTAc                         <  1:971691/150‑1 (MQ=255)
                                                                                                                            ttttCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACg                        <  2:268787/150‑1 (MQ=255)
                                                                                                                              ttCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGc                      <  2:392843/150‑1 (MQ=255)
                                                                                                                                 cATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGGCGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTc                   >  1:246294/1‑150 (MQ=255)
                                                                                                                                 cATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTc                   <  1:446747/150‑1 (MQ=255)
                                                                                                                                  cctacAAAACGTTTGCCGATCTGAAAGGCAAGCGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCa                  >  1:845599/5‑150 (MQ=255)
                                                                                                                                     ccAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTAGCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGCTTACGGCGTCagat               >  2:500528/1‑150 (MQ=255)
                                                                                                                                     ccAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATAtt              >  1:1684648/1‑150 (MQ=255)
                                                                                                                                      cAAAACGTTTGCCGATCTGAAAGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATAtt              >  2:1006807/1‑150 (MQ=255)
                                                                                                                                          aCGTTTGCCGATCTGAAAGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCCCGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCt          >  1:1585690/1‑150 (MQ=255)
                                                                                                                                          aCGTTTGCCGATCTGAAAGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCt          >  2:1385707/1‑150 (MQ=255)
                                                                                                                                            gTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGa        >  1:766231/1‑150 (MQ=255)
                                                                                                                                               tGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAAtt     >  1:318563/1‑150 (MQ=255)
                                                                                                                                                gCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAAttt    <  2:1647436/150‑1 (MQ=255)
                                                                                                                                                 ccTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTc   >  1:1734152/1‑150 (MQ=255)
                                                                                                                                                 ccTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTc   >  2:1070421/1‑150 (MQ=255)
                                                                                                                                                  cTTTCATATCGGCAAACGTTTTGTAGGCATATTTTTTGGCAATCACGACGGCGGAGTTATCATAGTAGGGCGAGGTAATCGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAAATTCa  <  1:341447/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCA  >  NC_000913/900054‑900349

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: