Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,004,029 C→A 8.1% M50I (ATG→ATT)  rarD ← putative chloramphenical resistance permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,0290CA 8.1% 862.6 / 11.2 260M50I (ATG→ATT) rarDputative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base C (169/60);  new base A (2/28);  total (171/89)
Fisher's exact test for biased strand distribution p-value = 6.61e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

AGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCC  >  NC_000913/4003880‑4004176
                                                                                                                                                     |                                                                                                                                                   
aGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCaccacc                                                                                                                                                     >  2:5718875/1‑150 (MQ=255)
  cATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCAt                                                                                                                                                   <  1:1055945/150‑1 (MQ=255)
   atatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATa                                                                                                                                                  <  2:4535638/150‑1 (MQ=255)
   atatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATa                                                                                                                                                  >  1:2004979/1‑150 (MQ=255)
    tatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCTGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATaa                                                                                                                                                 <  2:4470134/150‑1 (MQ=255)
    tatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATaa                                                                                                                                                 <  2:734290/150‑1 (MQ=255)
    tatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATaa                                                                                                                                                 <  1:5010763/150‑1 (MQ=255)
      tGGAGATTGTTCACCGCCCAGATAAACAGTAGACAGTAGACACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTATGTGGCTTCTGAATCAGAGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACAATAAAg                                                                                                                                               >  1:1788341/1‑150 (MQ=255)
        gTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGaa                                                                                                                                             >  2:1008397/1‑150 (MQ=255)
         tGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGaaa                                                                                                                                            <  2:1892313/150‑1 (MQ=255)
           aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGTCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGaaaaa                                                                                                                                          <  2:4885304/150‑1 (MQ=255)
           aTTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGaaaaa                                                                                                                                          <  2:4318225/150‑1 (MQ=255)
              gTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGa                                                                                                                                       <  2:1246636/150‑1 (MQ=255)
              gTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGa                                                                                                                                       >  1:5275439/1‑150 (MQ=255)
                  aCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAg                                                                                                                                   >  2:795947/1‑150 (MQ=255)
                  aCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAg                                                                                                                                   >  1:1170971/1‑150 (MQ=255)
                  aCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAg                                                                                                                                   >  1:4732570/1‑150 (MQ=255)
                      cccctaTAACCAGTTGCATGTTGCCACAAATCTGCAGGGCAGAGACTGCCAGAATAAAATTTTTCTGTGGCGTCTGAATCTGCTTTTTTAAATAGGACCACTGGGGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCa                                                                                                                               <  2:3267504/145‑1 (MQ=255)
                      cccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCa                                                                                                                               >  1:4705609/1‑150 (MQ=255)
                      cccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCa                                                                                                                               >  1:4154678/1‑150 (MQ=255)
                          gATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCAcgcg                                                                                                                           >  1:3391382/1‑150 (MQ=255)
                           aTAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATTCTCATCAGAACAAAAATAAAGAAAAACGACCAGAACAAGCGa                                                                                                                          >  2:5292251/1‑150 (MQ=255)
                             aaactttAGCCTGTTGCCACAAATCAGCAGGGCTGAGACTGCAAGAATAAAATTTTTCTGGGGCGTATGAATCAGCTTTTTTAACTTGCACCCCGGGGGGCAAATGCACATCACCACCCACATAAACAAAAACGACCAGATCACGCGATg                                                                                                                        <  2:30159/144‑1 (MQ=255)
                             aaaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATg                                                                                                                        >  1:925114/1‑150 (MQ=255)
                              aaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGc                                                                                                                       >  2:200094/1‑150 (MQ=255)
                                cAGTAGCCTGTCCCCACCAACCCGCCCGGACGATACGGCCACAAAAAAATTTTTGTGTGGCGTCTGAAGCAGGTTTTTTAAATCGCACCCCGGGGGGAAATTGATCACCAGAACAACAATAAACAAACACGACCAGATCACGCGATGCGt                                                                                                                     <  2:515276/150‑1 (MQ=255)
                                     gCCAGTTCCCACCAAGCCGCACGGCAGAGACTGCCACCATAAAATTTTTCTGTGGCGTCTTACCCCGCTTTTTTAATTTGGACCAGGGGCGGCAATTGCTCAACAGCACCACCATAACCAAACACGACCAGATCACGCGATGCGTCAGGa                                                                                                                <  2:2069411/150‑1 (MQ=255)
                                      ccAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGAt                                                                                                               >  2:1087677/1‑150 (MQ=255)
                                         gTTGCCACCACTCAGCACGGCAGAGACTGCCAGCATAAAAATTTGCTGTGGCGGCTGAATCAGCGTTTTTAAATCGGACCACGGGCGGCAAATGCTCATCCGCACCACCATAAAGCGAAACGACCAGATCACGCGATGCGTCAGGATTTc                                                                                                            <  2:5217398/150‑1 (MQ=255)
                                         gTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTc                                                                                                            >  2:4875838/1‑150 (MQ=255)
                                           tGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCAt                                                                                                          >  1:3858424/1‑150 (MQ=255)
                                           tGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCAt                                                                                                          >  1:1055866/1‑150 (MQ=255)
                                            cccaccaACCAGCACTGCATATACTGACAGAATAAAATTTTTCTGTGGCTTCTTAATCAGGTTTTTAAAATTGCAAACCTGGCGGCAAATGCTCATCAGCACCACAAAAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATc                                                                                                         <  2:167902/149‑1 (MQ=255)
                                              caccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcgg                                                                                                       >  1:4092481/1‑150 (MQ=255)
                                              caccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcgg                                                                                                       >  2:2563336/1‑150 (MQ=255)
                                               accaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcggc                                                                                                      >  1:7010/1‑150 (MQ=255)
                                                ccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcggcg                                                                                                     >  1:1829310/1‑150 (MQ=255)
                                                ccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGAACACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcggcg                                                                                                     >  1:4829896/1‑150 (MQ=255)
                                                  aaTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCggg                                                                                                   >  2:1841616/1‑150 (MQ=255)
                                                   atccccACGGCAGAGCCCCCAACCAAAAAAATTTCCTGTGGCGACTTAAACCGGTTTTTTAAATAGCACCACTGGCGGCAACTGCTCACCCGCACCACCATAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGc                                                                                                  <  2:2202179/145‑1 (MQ=255)
                                                   aTCAGCACGGACGCGCCGGCCAGAAAAAAATTGTTCTGTGGCGTCTGAATCCGCTTTTTTAAATAGCACCACTGGCGGCAAATGCACATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGc                                                                                                  <  2:2760817/150‑1 (MQ=255)
                                                      aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACg                                                                                               >  2:4618973/1‑150 (MQ=255)
                                                        cACGGCAGAGACTGCCAGCCTAAAAATTTTCTGTGGCGTCTGACTCAGCGTTTTTAAATAGGACCACGGGCGGCAAATGCTCATCAGCACCCCCCTAAAGAAAAACGACCAGATCACGCGAGGCGGCAGGAGTTCATCGGCGGGCCCggc                                                                                             >  2:3204861/1‑148 (MQ=255)
                                                        cACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACgta                                                                                             >  2:3516718/1‑150 (MQ=255)
                                                         agggCAGAGACGGCCAGCATAAAATTTTTCTGTGGCGTCTGAATCTGGGTTTTTAACTCGGTCCCGGGGCGGCAATTGCTCACCAGAACCACCATCACGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACgtag                                                                                            <  2:4692230/148‑1 (MQ=255)
                                                         aCGGCAGAGACTGCCAGAATAAAATTTTTCGGGGGCGTCTGAATCCGCTTTTTTAAATAGGACCCGTGGCGGCAACGGCTCATCAGCACCCCCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACgtag                                                                                            <  2:4599753/150‑1 (MQ=255)
                                                                    tGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACtt                                                                                 >  2:3908061/1‑150 (MQ=255)
                                                                    tGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACtt                                                                                 >  2:36230/1‑150 (MQ=255)
                                                                     gCCACCATAAAATTTTTCTGTGGCGTCTCAATCAGGTTTTTAAAATCGCACCCGTGGCGCCAAAGGCCCATCACCACCACCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTg                                                                                <  2:52256/150‑1 (MQ=255)
                                                                      ccAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGa                                                                               >  2:4799870/1‑150 (MQ=255)
                                                                      ccAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGAAGCGTCAGGATTTCATCGGCGGGCCCGTAGTAAATCAACTTGa                                                                               >  1:4396279/1‑150 (MQ=255)
                                                                       cAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGaa                                                                              >  2:4248316/1‑150 (MQ=255)
                                                                       cAGAATAAAATTTTTCTGTGGCGTCTTAATCAGCTTTTTTAAATAGGACCACTGGCGGAAAATGCTAATCACCACCACAATAAAGAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGaa                                                                              <  2:4096752/150‑1 (MQ=255)
                                                                        agaaTAAAATTTTTGTGTGGCTTCTGAATCTGCTTTTTTAAATAGCACCCGGGGCGGAAAATGCTCATCACCACCACAATAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAg                                                                             <  1:364837/147‑1 (MQ=255)
                                                                           aTAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAc                                                                          >  2:4752807/1‑150 (MQ=255)
                                                                            tAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                         >  2:3814522/1‑150 (MQ=255)
                                                                            tAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                         >  2:2124788/1‑150 (MQ=255)
                                                                             aaaatTTTTGTGGGGCGTGTCAATCTGCTTTTTTAAATTGGACCTCTGGGGGCAACTGCTCAACACCACCACAATAAATAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACATAGTAAATCAACTTGAAGTACGc                                                                        <  1:2869696/150‑1 (MQ=255)
                                                                              aaaaTTTTCTGTGGCGTCTTAATCTGCGTTTTTAAATCGCACCAGTGGCGGAAAATCCAAATCCGCACCCCAATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCt                                                                       <  2:5470410/150‑1 (MQ=255)
                                                                              aaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCt                                                                       >  2:1945351/1‑150 (MQ=255)
                                                                              aaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCt                                                                       >  2:1773860/1‑150 (MQ=255)
                                                                                aaTTTTCTGTGGCGTCTGAATCAGCTTTTTTAATTAGGACCACTGGCGGCAAATGCTAATCAGCACCCCCATAATGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTgg                                                                     <  2:1472051/150‑1 (MQ=255)
                                                                                 tttttGTGGGGCTTCTTAATCAGCTTTTTAAAATCGCACCACGGGCGGAAAATGCTAAACAGCACCACAAAAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgc                                                                    <  1:1382443/149‑1 (MQ=255)
                                                                                 aTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgc                                                                    >  1:4937053/1‑150 (MQ=255)
                                                                                   tttCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                  >  1:5513296/1‑150 (MQ=255)
                                                                                   tttCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                  >  1:2886606/1‑150 (MQ=255)
                                                                                          ggCGTCTTAATCAGCTTTTTAAAATTGCCCCACTGGGGGAAACTGCTAATCACCACCACAATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                           <  1:3885999/150‑1 (MQ=255)
                                                                                           gCGTCTTAATCCGCTTTTTTAATTCGGACCCGTGGGGGAAATGGCTCATCCGCCCCACAAAAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcccc                                                          <  1:2202352/150‑1 (MQ=255)
                                                                                              tCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                       >  2:2388596/1‑150 (MQ=255)
                                                                                                    tCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaat                                                 >  1:596409/1‑150 (MQ=255)
                                                                                                       gCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                              >  2:426194/1‑150 (MQ=255)
                                                                                                       gCGTTTTTAAATAGGACCACGGGCGGAAAATGCTCATCAGCACAACAATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                              <  1:5723425/150‑1 (MQ=255)
                                                                                                          tttttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                           >  1:2765725/1‑150 (MQ=255)
                                                                                                            tttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGATCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                         >  2:2985685/1‑150 (MQ=255)
                                                                                                              tAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                       >  1:4481553/1‑150 (MQ=255)
                                                                                                              tAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                       >  1:3020186/1‑150 (MQ=255)
                                                                                                                aaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                     >  2:1609474/1‑150 (MQ=255)
                                                                                                                 aTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                    >  1:1638104/1‑150 (MQ=255)
                                                                                                                 aTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                    >  2:5661541/1‑150 (MQ=255)
                                                                                                                 aTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATGTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                    >  2:4426072/1‑150 (MQ=255)
                                                                                                                   tggACCACGGGGGGAAATTGCTCATCTGCACCACAATAATGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                  <  2:4579530/149‑1 (MQ=255)
                                                                                                                   aGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                  >  2:2713836/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  2:4615572/1‑150 (MQ=255)
                                                                                                                       cctgtGGCGCAAATTGCTAAACAGCACCACAAAAAACAAAAACGACCAGATCCCGCGTTGCGTCAGGATTTCATCGGGGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                              <  1:819042/146‑1 (MQ=255)
                                                                                                                                gcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                     >  1:1331172/1‑150 (MQ=255)
                                                                                                                                  aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                   >  2:478415/1‑150 (MQ=255)
                                                                                                                                      gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               >  1:881921/1‑150 (MQ=255)
                                                                                                                                      gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               >  2:1425602/1‑150 (MQ=255)
                                                                                                                                      gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               >  1:3599382/1‑150 (MQ=255)
                                                                                                                                       cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt              >  1:509081/1‑150 (MQ=255)
                                                                                                                                           tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGTGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt          >  2:1386375/1‑150 (MQ=255)
                                                                                                                                           tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt          >  2:4191436/1‑150 (MQ=255)
                                                                                                                                            ccgCACCACCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt         <  2:674344/148‑1 (MQ=255)
                                                                                                                                            caccacaacAATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACTTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt         <  1:3205297/147‑1 (MQ=255)
                                                                                                                                            caGCACCACCATAACCCACAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCACCTTGAAGTACGCTGGCGCTATACCCCAAAGAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt         <  2:2999131/150‑1 (MQ=255)
                                                                                                                                             aGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt        >  2:5810894/1‑150 (MQ=255)
                                                                                                                                              gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTg       >  1:2529734/1‑150 (MQ=255)
                                                                                                                                               caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTAGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc      >  1:3721241/1‑150 (MQ=255)
                                                                                                                                                   cccaTAACTAACAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATACAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGGGTTTGTTTTGcatcc  <  2:4396279/149‑1 (MQ=255)
                                                                                                                                                   accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc  <  2:2405926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCC  >  NC_000913/4003880‑4004176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑