Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 900,202 T→A 6.6% 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.6% 308.9 / 8.6 106K125I (AAA→ATA) artJarginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (55/44);  new base A (2/5);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 109 total aligned reads displayed.

GTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGG  >  NC_000913/900054‑900351
                                                                                                                                                    |                                                                                                                                                     
gtCACCAAATACCCCATCAATACGACCATATTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGGGACAGTTTTCACGTCCGGGTGCTGATCCTGATTATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGACttt                                                                                                                                                      >  2:79241/1‑150 (MQ=255)
   aCCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAg                                                                                                                                                   <  2:1875730/150‑1 (MQ=255)
    ccAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGa                                                                                                                                                  >  2:624199/1‑150 (MQ=255)
    ccAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCATTTCAGa                                                                                                                                                  >  2:693736/1‑150 (MQ=255)
      aaaTACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATc                                                                                                                                                <  2:1609452/150‑1 (MQ=255)
      aaaTACCCAATCAATACGAACATTTTTCATATCTATAAAGGCATTCTGATAACTGTCATAAGAGACTGTTTTCACTTCCGGGTGCTTATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCGATCCCAATACGTTTGCCTTTCAGATc                                                                                                                                                <  1:388611/150‑1 (MQ=255)
       aaTACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATACGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCg                                                                                                                                               <  2:1652843/150‑1 (MQ=255)
        aTACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCgg                                                                                                                                              <  1:553542/150‑1 (MQ=255)
           ccccATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCaa                                                                                                                                           >  2:539823/1‑150 (MQ=255)
              cATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACg                                                                                                                                        <  2:681827/150‑1 (MQ=255)
              cATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATAATTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACg                                                                                                                                        <  2:1714800/150‑1 (MQ=255)
                 cAATACGACCATTTTTCAGATCGATAAAGGCATTCTTATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGtat                                                                                                                                     <  1:1061238/150‑3 (MQ=255)
                 aaaTACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTCCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGtat                                                                                                                                     <  1:115999/149‑3 (MQ=255)
                  aaTACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGtttt                                                                                                                                    <  2:558016/150‑1 (MQ=255)
                  aaTACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGtatt                                                                                                                                    <  2:458658/150‑1 (MQ=255)
                     aCGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTTAATATATTTCTGGTGCGTAGTAACGTTTTCCATCCCAATACGTTTGACTTTCAGATCGGCAAACGTTTTGTa                                                                                                                                 >  1:256158/1‑150 (MQ=255)
                      cGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAg                                                                                                                                <  1:1312579/150‑1 (MQ=255)
                       tatcATTTTTCAGATCTATAAAGTCATTCTGATAACTGTCATAAGAGTCAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAgg                                                                                                                               <  1:1603411/147‑1 (MQ=255)
                       gACCATTTTTCAGATCGATAAAGGCATTCTGATAAATGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAgg                                                                                                                               <  1:918311/150‑1 (MQ=255)
                       gACCATTTTTCAGATCGAAAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCGATCTGAAAGGCAAACGTATTgggat                                                                                                                               <  2:350919/150‑5 (MQ=38)
                             ttttCAGATCGATAAAGGCTTTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtt                                                                                                                         <  1:449081/150‑1 (MQ=255)
                             ttttCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtt                                                                                                                         >  1:998905/1‑150 (MQ=255)
                               ttCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtttt                                                                                                                       >  1:958265/1‑150 (MQ=255)
                                 cAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCtttttt                                                                                                                     >  2:194227/1‑150 (MQ=255)
                                  aGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTg                                                                                                                    <  2:1119665/150‑1 (MQ=255)
                                   gATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTgg                                                                                                                   >  1:1777941/1‑150 (MQ=255)
                                   gATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAATACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTgg                                                                                                                   <  2:386298/150‑1 (MQ=255)
                                    aTCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTgga                                                                                                                  >  2:653706/1‑149 (MQ=255)
                                    aTCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGc                                                                                                                  >  2:799118/1‑150 (MQ=255)
                                    aTCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGc                                                                                                                  <  1:1489065/150‑1 (MQ=255)
                                     tCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCa                                                                                                                 <  2:462611/150‑1 (MQ=255)
                                       gATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAAt                                                                                                               >  2:1351427/1‑150 (MQ=255)
                                         tAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCa                                                                                                             >  1:1045545/1‑150 (MQ=255)
                                           aaGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATAGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCacg                                                                                                           >  1:1498242/1‑150 (MQ=255)
                                            aGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTAACGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCacga                                                                                                          <  2:1056136/150‑1 (MQ=255)
                                                aTTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCATGACGGc                                                                                                      >  1:1843639/1‑150 (MQ=255)
                                                    tGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAg                                                                                                  >  2:383678/1‑150 (MQ=255)
                                                        aaCTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGtttt                                                                                              <  2:1060558/150‑1 (MQ=255)
                                                         aCTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTc                                                                                             <  1:195635/150‑1 (MQ=255)
                                                          cTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTGTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCa                                                                                            >  2:1477547/1‑150 (MQ=255)
                                                           tGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAt                                                                                           >  2:898534/1‑150 (MQ=255)
                                                            gTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAta                                                                                          >  2:193952/1‑150 (MQ=255)
                                                             tCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtag                                                                                         >  1:1571836/1‑150 (MQ=255)
                                                             tCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtag                                                                                         >  1:1564519/1‑150 (MQ=255)
                                                               aTAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtagta                                                                                       >  1:566723/1‑150 (MQ=255)
                                                                tAAGAGACAGTTTTCACTTCCGTGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtagtag                                                                                      >  2:466299/1‑150 (MQ=255)
                                                                tAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCAtagtag                                                                                      <  1:1440909/150‑1 (MQ=255)
                                                                 aaGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAgg                                                                                     >  2:1820861/1‑150 (MQ=255)
                                                                  agagACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                    >  2:1766210/1‑150 (MQ=255)
                                                                  agagACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                    <  2:1401330/150‑1 (MQ=255)
                                                                  agagACAGTTTTCACTTACGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAggg                                                                                    <  2:1393711/150‑1 (MQ=255)
                                                                     gaCAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGt                                                                                 >  2:572586/1‑150 (MQ=255)
                                                                     gaCAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGt                                                                                 >  2:458385/1‑150 (MQ=255)
                                                                     gaCAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGt                                                                                 >  2:1498461/1‑150 (MQ=255)
                                                                       aaGTTTTCAATTCCGGGTGCTGCTCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATTCGTTTGCCTTTCAGATCGTCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTgg                                                                               <  2:177432/149‑1 (MQ=255)
                                                                         gTTTTCATTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGCCGGCTGAGTTTTCATAGTAGGGCGTGGTa                                                                             >  2:1386561/1‑150 (MQ=255)
                                                                           tttCACTTCCGGGTGCTGCTCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGGTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTaaa                                                                           >  2:512231/1‑150 (MQ=255)
                                                                                cTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATa                                                                      <  1:1901945/150‑1 (MQ=255)
                                                                                  tCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATAcc                                                                    >  1:1229945/1‑150 (MQ=255)
                                                                                  tCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATAcc                                                                    <  2:1859979/150‑1 (MQ=255)
                                                                                      ggTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGGGTTTTCATAGTAGGGCGTGGTAAACGATACCTGtt                                                                <  2:1886149/150‑1 (MQ=255)
                                                                                      ggTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGATAAACGATACCTGtt                                                                <  1:1002104/150‑1 (MQ=255)
                                                                                        tGCTGATCCTGAATATATTTATGGTGCGTAGTCCCGGTGTCCCTCCCAATACGTTTGCCTTTCTGATCTGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGCCTGAGTTTTCATAGTAGGGCGGGGTAAACGATACCTGTTTg                                                              <  2:307444/150‑1 (MQ=255)
                                                                                           tGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTa                                                           >  1:676597/1‑150 (MQ=255)
                                                                                            gATCCTTAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTAc                                                          <  1:1790268/150‑1 (MQ=255)
                                                                                                 tGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTGAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCAAGACGGCTGAGTTTTCATAGTAGGGCGTTGTAAACGATACCTGTTTGCTACGCTcc                                                     >  1:357108/1‑150 (MQ=255)
                                                                                                  gaattaatTTCTAGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGGGGTAAACGATACCTGTTTGCTACGCTCCg                                                    <  1:101393/144‑1 (MQ=255)
                                                                                                  gAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCg                                                    >  1:135339/1‑150 (MQ=255)
                                                                                                  gAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCg                                                    <  1:146646/150‑1 (MQ=255)
                                                                                                  gAATAGATTTCTGGAGCGTAGTACAGTTTTCCATCCCTATACGTGTGCCTATCAGCGCGGCAAACGTTTTGTATGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGGGTGGTCAACGATACCTGTTTGCTACGCTCCg                                                    <  1:105909/150‑1 (MQ=255)
                                                                                                    atatatTTCTGGTGCGTAGTACGGTTTGCCTTCCCAATACGGTTGCCTTTCAGATCGGAAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCggg                                                  <  2:1683571/150‑1 (MQ=255)
                                                                                                      atatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGt                                                <  1:213340/150‑1 (MQ=255)
                                                                                                      atatTTCTGGTGAGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGt                                                <  1:1168560/150‑1 (MQ=255)
                                                                                                            cTGGTGCGTAGTACCGGTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGTTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATc                                          <  1:828806/150‑1 (MQ=255)
                                                                                                              ggTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGTGGTGATATCCa                                        >  1:1403667/1‑150 (MQ=255)
                                                                                                              ggTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCa                                        >  2:946301/1‑150 (MQ=255)
                                                                                                              ggTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCa                                        >  2:401682/1‑150 (MQ=255)
                                                                                                                 gCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATAc                                     >  1:1694508/1‑150 (MQ=255)
                                                                                                                   gtagtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTATACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCg                                   >  2:44719/1‑150 (MQ=255)
                                                                                                                    tagtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCgg                                  <  1:1094645/150‑1 (MQ=255)
                                                                                                                     agtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGga                                 <  1:539823/150‑1 (MQ=255)
                                                                                                                       taCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGgaga                               >  1:803830/1‑150 (MQ=255)
                                                                                                                         ccGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGAtt                             >  1:1831575/1‑150 (MQ=255)
                                                                                                                         ccGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGAtt                             >  1:951540/1‑150 (MQ=255)
                                                                                                                         ccGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGAtt                             >  1:1156310/1‑150 (MQ=255)
                                                                                                                           gTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTAc                           >  1:885472/1‑150 (MQ=255)
                                                                                                                            ttttCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACg                          <  1:1338762/150‑1 (MQ=255)
                                                                                                                             tttCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACgg                         <  2:1777941/150‑1 (MQ=255)
                                                                                                                               tCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCg                       <  1:458385/150‑1 (MQ=255)
                                                                                                                               tCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCg                       >  1:1261104/1‑150 (MQ=255)
                                                                                                                                ccATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGTGATTACGGCGt                      >  2:1133817/1‑150 (MQ=255)
                                                                                                                                ccATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGt                      <  1:712363/150‑1 (MQ=255)
                                                                                                                                 cATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTc                     <  2:1214783/150‑1 (MQ=255)
                                                                                                                                     ccAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCatat                 >  2:1809884/1‑150 (MQ=255)
                                                                                                                                     acAAAACGTTTGCCGATCTGAAAGGCAAACGTATTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAATCGATACCTGTTTGCTACGCTCCGGGGTGATATCCATTCCGGAGATTTCGGCGTCatat                 >  1:1303216/2‑150 (MQ=38)
                                                                                                                                      cAAAACGTTTGCCGATCTGAAAGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTATTCATAGTACGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATAtt                >  1:1256573/1‑150 (MQ=255)
                                                                                                                                           cGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTg           >  1:1665711/1‑150 (MQ=255)
                                                                                                                                             tttGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGaa         >  2:1855465/1‑150 (MQ=255)
                                                                                                                                                 cgtttCTGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTc     <  2:1498242/148‑1 (MQ=255)
                                                                                                                                                    ttCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAgg  <  1:676242/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGG  >  NC_000913/900054‑900351

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: ‑