Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,816 T→G 8.6% N190K (AAT→AAG) ‡ yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,8160TG 8.6% 397.8 / 21.2 135N190K (AAT→AAG) ‡yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (63/60);  new base G (5/7);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 135 total aligned reads displayed.

GGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAC  >  NC_000913/4062668‑4062962
                                                                                                                                                    |                                                                                                                                                  
ggCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATg                                                                                                                                                   <  2:1315679/150‑1 (MQ=255)
ggCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATg                                                                                                                                                   <  2:1453492/150‑1 (MQ=255)
  cGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGAc                                                                                                                                                 >  2:1342688/1‑150 (MQ=255)
  cGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTAAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGAc                                                                                                                                                 >  1:526625/1‑150 (MQ=255)
   gTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACa                                                                                                                                                >  2:1717430/1‑150 (MQ=255)
    ttgccttGTCGCTCCCGGGGTGGCATTTTCTGGCTTGGCAGTCTGTCCCTGGTGTTTCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATACCCTTATTGTGATTGCCGTAGGCATTATTAATATCTTTGTCATGAATGACaa                                                                                                                                               <  1:280571/145‑1 (MQ=255)
      gCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaag                                                                                                                                             <  1:1211295/150‑1 (MQ=255)
       cATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaaga                                                                                                                                            >  1:1505722/1‑150 (MQ=255)
         ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAaagaag                                                                                                                                          <  2:647424/150‑1 (MQ=255)
         ttGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaa                                                                                                                                          <  1:797211/150‑2 (MQ=255)
           ggcGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                        <  2:876206/148‑4 (MQ=255)
           gTCGATACCGTGGTGGCATTTTCTGCGTTGACAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCATAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaaaat                                                                                                                                        <  1:592606/150‑4 (MQ=255)
              gATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCa                                                                                                                                     >  1:758403/1‑150 (MQ=255)
                tACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAcc                                                                                                                                   >  2:866769/1‑150 (MQ=255)
                   cGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtc                                                                                                                                >  1:19443/1‑150 (MQ=255)
                       gtggCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGtt                                                                                                                            >  2:1202140/1‑150 (MQ=255)
                         ggCTTTTTCTGAGTTGGCAGTATTTAGCTGGTTTGGGTGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCGCTTATTGTGATTGCCGTAGGCTTTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGttta                                                                                                                          <  2:428655/150‑3 (MQ=255)
                          gCATTTTCTGCGTTGGCAGTATTTACCAGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAg                                                                                                                         <  2:870555/150‑1 (MQ=255)
                           cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAgg                                                                                                                        <  2:1509170/150‑1 (MQ=255)
                           cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaa                                                                                                                                           <  1:1517276/131‑1 (MQ=255)
                           cATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCATTCAGGACAaagaa                                                                                                                                           >  2:1517276/1‑131 (MQ=255)
                              tttCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGa                                                                                                                     <  2:215999/150‑1 (MQ=255)
                               ttCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGaa                                                                                                                    >  2:1370717/1‑150 (MQ=255)
                                 cTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCAccgtc                                                                                                                                <  1:463900/136‑1 (MQ=255)
                                 cTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGaaaa                                                                                                                  >  1:1446461/1‑150 (MQ=255)
                                    cGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTCTGTTTTAAAGCAGGCATCTGGTTCTATTACCTTATTGTGATTGCCGTAGGAATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAaag                                                                                                               <  2:661015/150‑1 (MQ=255)
                                    cGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGTAGAGGCACCGTCCGTTGAGGTGAAAAaag                                                                                                               <  1:428592/150‑1 (MQ=255)
                                     gTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAaaga                                                                                                              <  2:1290228/150‑1 (MQ=255)
                                     gTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAaaga                                                                                                              <  1:1235678/150‑1 (MQ=255)
                                        ggCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAagaaga                                                                                                           >  2:549584/1‑150 (MQ=255)
                                         gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGAc                                                                                                          <  2:1419008/150‑1 (MQ=255)
                                         gCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGAc                                                                                                          <  1:1342688/150‑1 (MQ=255)
                                          cAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACg                                                                                                         >  2:1568050/1‑150 (MQ=255)
                                             tATTTAGCTGTTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCCTCTGGGTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTATTTGTCCTGAATGAAAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAg                                                                                                      <  2:175990/150‑1 (MQ=255)
                                             tATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAg                                                                                                      <  2:969935/150‑1 (MQ=255)
                                              aTTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGc                                                                                                     >  2:1664607/1‑150 (MQ=255)
                                              aTTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGc                                                                                                     >  2:781240/1‑150 (MQ=255)
                                               tttACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTAGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCa                                                                                                    >  2:1069927/1‑150 (MQ=255)
                                                  aCCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCg                                                                                                 >  1:625291/1‑150 (MQ=255)
                                                    cTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaa                                                                                               >  2:944970/1‑150 (MQ=255)
                                                     tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaa                                                                                              >  1:958479/1‑150 (MQ=255)
                                                     tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGAACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaa                                                                                              >  2:824321/1‑150 (MQ=255)
                                                     tGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGACGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAACAAGCCGGAGCTTCGaaa                                                                                              >  2:240358/1‑150 (MQ=255)
                                                      ggTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaa                                                                                             >  2:705067/1‑150 (MQ=255)
                                                      ggTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaa                                                                                             >  2:786538/1‑150 (MQ=255)
                                                       gTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGaaaaa                                                                                            >  2:943536/1‑150 (MQ=255)
                                                        tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAac                                                                                           >  2:339416/1‑150 (MQ=255)
                                                        tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAac                                                                                           <  1:723128/150‑1 (MQ=255)
                                                        tttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACAGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAac                                                                                           >  1:846929/1‑150 (MQ=255)
                                                         ttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAaca                                                                                          <  1:994105/150‑1 (MQ=255)
                                                         ttGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAaca                                                                                          <  1:866769/150‑1 (MQ=255)
                                                           ggCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACAcc                                                                                        >  1:792146/1‑150 (MQ=255)
                                                             cAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAg                                                                                      <  1:1295374/150‑1 (MQ=255)
                                                             cAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAg                                                                                      >  2:1614700/1‑150 (MQ=255)
                                                                tGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTAt                                                                                   >  1:1655881/1‑150 (MQ=255)
                                                                 ggCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCTTTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCCAAAAACACCAGTATg                                                                                  >  2:228049/1‑150 (MQ=255)
                                                                  gCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGa                                                                                 <  2:1383192/150‑1 (MQ=255)
                                                                   cTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTGTGAc                                                                                <  2:711001/150‑1 (MQ=255)
                                                                   cTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGAc                                                                                <  1:54268/150‑1 (MQ=255)
                                                                      attaGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTc                                                                             <  2:846929/150‑1 (MQ=255)
                                                                        tacgTTTTAAAGCAGGCATCTGGTTCTATACCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  2:795594/4‑150 (MQ=255)
                                                                        taGGTTTTTAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  2:1607572/1‑150 (MQ=255)
                                                                        taGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           <  1:1001037/150‑1 (MQ=255)
                                                                        taGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCgg                                                                           >  1:1388464/1‑150 (MQ=255)
                                                                          gggtttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                         <  2:1225633/147‑1 (MQ=255)
                                                                          ggTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTg                                                                         >  1:368112/1‑150 (MQ=255)
                                                                           gTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAATGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGc                                                                        >  1:1097888/1‑150 (MQ=255)
                                                                             tttAAAGCAGGCATCTGGTTCTATTCTCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTg                                                                      >  2:1448180/1‑150 (MQ=255)
                                                                              ttAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGa                                                                     <  2:19443/150‑1 (MQ=255)
                                                                                 aaGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAg                                                                  <  2:973584/150‑1 (MQ=255)
                                                                                      ggCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACaaa                                                             <  2:463959/150‑1 (MQ=255)
                                                                                        cATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAAc                                                           >  1:426498/1‑150 (MQ=255)
                                                                                         aTCTGGTTCTATTACCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACATCGGTGCTGAAAGACAAAACt                                                          >  2:1081234/1‑150 (MQ=255)
                                                                                            tGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTATCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGATAGACAAAACTAta                                                       >  2:705174/1‑149 (MQ=255)
                                                                                             ggTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCt                                                      <  2:776699/150‑1 (MQ=255)
                                                                                             ggTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCt                                                      <  2:645165/150‑1 (MQ=255)
                                                                                              gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     <  1:1369465/150‑1 (MQ=255)
                                                                                              gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     >  1:1588387/1‑150 (MQ=255)
                                                                                              gTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTg                                                     <  1:1016819/150‑1 (MQ=255)
                                                                                                     tCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATc                                              <  1:1615805/150‑1 (MQ=255)
                                                                                                      cccTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCg                                             >  2:1686363/1‑150 (MQ=255)
                                                                                                          tATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCttt                                         <  1:1359719/150‑1 (MQ=255)
                                                                                                            ttGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTa                                       <  2:1265056/150‑1 (MQ=255)
                                                                                                                gATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGTCAAAACTATCTGGCTTATCGCTTTTAACGt                                   <  2:1364806/150‑1 (MQ=255)
                                                                                                                 aTTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGATGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGtc                                  >  1:224182/1‑150 (MQ=255)
                                                                                                                       gTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCttcttc                            <  2:1627563/150‑1 (MQ=255)
                                                                                                                        tAGTCATTTTTTTTTTATTTGTCCTGAATGACATAGAAGTGGCACCGTCGTTTGAGGTGAAAATAGTAGACGGAGCATCGTAAAACACCAGTATGTCCTCGGTGCTGATAGACAAAACTATCTGGCTTATCTCTTTTAACGTCTTCTTCg                           <  2:908732/150‑1 (MQ=255)
                                                                                                                          ggCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGtt                         <  1:795594/150‑1 (MQ=255)
                                                                                                                             attattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTAc                      >  1:1697546/1‑150 (MQ=255)
                                                                                                                                  tattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGt                 <  1:1202140/150‑1 (MQ=255)
                                                                                                                                  cctcttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGt                 >  1:1715750/5‑150 (MQ=255)
                                                                                                                                   ctcttCTTTGTCATTCAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGtt                >  2:1428479/4‑150 (MQ=255)
                                                                                                                                        cTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtg           >  1:669803/1‑150 (MQ=255)
                                                                                                                                        cTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtg           <  1:1568050/150‑1 (MQ=255)
                                                                                                                                         tttGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACtga          >  1:9244/1‑149 (MQ=255)
                                                                                                                                           tGTCCTGAATGACAAAGAAGATGCACCTTCCTTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTgg        <  1:705067/150‑1 (MQ=255)
                                                                                                                                           tGTCCTGAAGGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGCCGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTgg        >  2:1119914/1‑150 (MQ=255)
                                                                                                                                             tgcTGAATGACAAAGAAGAGGATCCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGcc      <  2:1655881/148‑1 (MQ=255)
                                                                                                                                              ccTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCt     <  2:1446461/150‑1 (MQ=255)
                                                                                                                                                 gAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAc  >  2:459332/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGAC  >  NC_000913/4062668‑4062962

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