Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,969 T→A 5.8% Y168N (TAC→AAC)  hisB → histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,9690TA 5.8% 357.1 / 9.6 103Y168N (TAC→AAC) hisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (47/50);  new base A (4/2);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 105 total aligned reads displayed.

AGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATG  >  NC_000913/2093821‑2094115
                                                                                                                                                    |                                                                                                                                                  
aGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                                                                                                                    >  3:1555/1‑149 (MQ=255)
aGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                                                                                                                    >  3:20298/1‑149 (MQ=255)
   gATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACg                                                                                                                                                 >  4:764383/1‑149 (MQ=255)
   gATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACg                                                                                                                                                 >  4:535926/1‑149 (MQ=255)
   gATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACg                                                                                                                                                 >  4:441382/1‑149 (MQ=255)
       gATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCa                                                                                                                                             <  2:67364/149‑1 (MQ=255)
        aTCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCAc                                                                                                                                            <  3:631613/149‑1 (MQ=255)
             gcTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGt                                                                                                                                       <  3:860535/149‑1 (MQ=255)
               tAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGc                                                                                                                                     >  2:297869/1‑149 (MQ=255)
               tAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGc                                                                                                                                     >  4:203466/1‑149 (MQ=255)
                 aCAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                                                                                                                   <  4:671420/149‑1 (MQ=255)
                     ttATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATa                                                                                                                               <  4:73195/149‑1 (MQ=255)
                     ttATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATa                                                                                                                               >  2:247925/1‑149 (MQ=255)
                      tATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATAc                                                                                                                              <  2:19572/149‑1 (MQ=255)
                        tgtgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                                                                                                                            >  1:110874/1‑149 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTTACGTAGTGCGTAATACCaaa                                                                                                                          <  2:565420/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  1:450714/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  3:441382/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  4:579513/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  3:382093/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  1:887120/149‑1 (MQ=255)
                          tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  2:110874/149‑1 (MQ=255)
                           gATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAag                                                                                                                         >  3:377970/1‑149 (MQ=255)
                            aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                        <  2:211237/149‑1 (MQ=255)
                            aTTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga                                                                                                                        >  3:130509/1‑149 (MQ=255)
                                  gATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa                                                                                                                  >  3:491909/1‑149 (MQ=255)
                                   aTCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAg                                                                                                                 >  4:390339/1‑149 (MQ=255)
                                         cgaccgacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                                                                                           <  4:377970/149‑1 (MQ=255)
                                         cgaccgacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                                                                                           <  2:385688/149‑1 (MQ=255)
                                         cgaccgacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                                                                                           <  2:325950/149‑1 (MQ=255)
                                         cgaccgacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                                                                                           <  1:358357/149‑1 (MQ=255)
                                              gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          <  4:465954/129‑1 (MQ=255)
                                              gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                          >  3:465954/1‑129 (MQ=255)
                                              gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                      >  2:114442/1‑149 (MQ=255)
                                              gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                      >  1:123363/1‑149 (MQ=255)
                                                 aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                   >  3:418460/1‑149 (MQ=255)
                                                 aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                   >  1:18206/1‑149 (MQ=255)
                                                 aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                   <  4:388523/149‑1 (MQ=255)
                                                  ttCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtg                                                                                                  >  4:46509/1‑149 (MQ=255)
                                                  ttCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtg                                                                                                  >  3:130827/1‑149 (MQ=255)
                                                  ttCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACACTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtg                                                                                                  >  2:539013/1‑149 (MQ=255)
                                                      aCTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                                                                                              <  3:133218/149‑1 (MQ=255)
                                                         ggCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTgg                                                                                           <  3:5194/149‑1 (MQ=255)
                                                            ggAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                                                                                        <  4:159580/149‑1 (MQ=255)
                                                             gAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCg                                                                                       >  1:246167/1‑149 (MQ=255)
                                                                 aCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                   <  1:622377/149‑1 (MQ=255)
                                                                  cATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAg                                                                                  >  2:894305/1‑149 (MQ=255)
                                                                    tGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGt                                                                                <  3:493104/149‑1 (MQ=255)
                                                                          ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagc                                                                          >  1:528938/1‑149 (MQ=255)
                                                                           tACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagca                                                                         >  1:868981/1‑149 (MQ=255)
                                                                            aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCaa                                                                        >  1:198098/1‑149 (MQ=255)
                                                                             cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                                                                                                     >  2:452951/1‑119 (MQ=255)
                                                                             cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                                                                                                     <  1:452951/119‑1 (MQ=255)
                                                                             cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAg                                                                       <  4:130509/149‑1 (MQ=255)
                                                                              tGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGa                                                                      >  1:199981/1‑149 (MQ=255)
                                                                                 tttACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTa                                                                   >  4:731781/1‑149 (MQ=255)
                                                                                    aCGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAaca                                                                <  1:174207/149‑1 (MQ=255)
                                                                                        tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCgg                                                            >  3:80416/1‑149 (MQ=255)
                                                                                        tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCgg                                                            >  3:800866/1‑149 (MQ=255)
                                                                                              cgcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTgg                                                      <  1:114442/149‑1 (MQ=255)
                                                                                              cgcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTgg                                                      <  3:390339/149‑1 (MQ=255)
                                                                                              cgcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTgg                                                      <  2:605693/149‑1 (MQ=255)
                                                                                                 gAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGctt                                                   <  3:40258/149‑1 (MQ=255)
                                                                                                  aaaCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttc                                                  <  2:123363/149‑1 (MQ=255)
                                                                                                   aaCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttct                                                 <  1:297869/149‑1 (MQ=255)
                                                                                                     cccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                               >  3:251146/1‑149 (MQ=255)
                                                                                                     cccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                               <  3:784221/149‑1 (MQ=255)
                                                                                                     cccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCttt                                               >  1:239658/1‑149 (MQ=255)
                                                                                                      ccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACTAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCCAGATTAACACCAGGGTTGGCTTCtgtg                                              <  4:434789/149‑1 (MQ=255)
                                                                                                      ccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTg                                              >  1:452498/1‑149 (MQ=255)
                                                                                                       cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                             <  2:868981/149‑1 (MQ=255)
                                                                                                       cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                             <  3:203466/149‑1 (MQ=255)
                                                                                                       cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                             >  1:96036/1‑149 (MQ=255)
                                                                                                        tGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                            <  1:457127/149‑1 (MQ=255)
                                                                                                         gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATc                                           <  3:97516/149‑1 (MQ=255)
                                                                                                         gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATc                                           <  2:18206/149‑1 (MQ=255)
                                                                                                               gCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGc                                     <  2:96036/149‑1 (MQ=255)
                                                                                                                  aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTgg                                  >  1:690348/1‑149 (MQ=255)
                                                                                                                  aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTgg                                  >  1:562278/1‑149 (MQ=255)
                                                                                                                   aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGa                                 <  4:20298/149‑1 (MQ=255)
                                                                                                                      aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCa                              >  4:519926/1‑149 (MQ=255)
                                                                                                                        tGGCGAGCAACTCACCAGACGTGACCGTTTCGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGa                            <  2:616484/149‑1 (MQ=255)
                                                                                                                               cAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAc                     <  3:450380/149‑1 (MQ=255)
                                                                                                                                aaCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAcc                    >  1:445271/1‑149 (MQ=255)
                                                                                                                                     aCCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcgg               >  3:712822/1‑149 (MQ=255)
                                                                                                                                      cactaCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggc              >  3:469049/5‑149 (MQ=255)
                                                                                                                                       cAGACGTGACGGTTACGCTCACGTAGGGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTTGCTGGCTCGTCAAGGTGGCAGCTAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCActgcg             >  4:125968/1‑149 (MQ=255)
                                                                                                                                       cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTTCCCAcggcg             >  3:340003/1‑149 (MQ=255)
                                                                                                                                       cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg             >  1:291270/1‑149 (MQ=255)
                                                                                                                                        aGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcgg            >  3:297546/1‑149 (MQ=255)
                                                                                                                                         taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                   <  1:616988/76‑1 (MQ=39)
                                                                                                                                         taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                   >  2:616988/2‑77 (MQ=39)
                                                                                                                                          aCGTGACCGTTACGCTGACGTAGTGAGTAATACCAAAGATACTGAGATTGACGTTCAGGTGTGGCTGTATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGtg          <  4:597449/149‑2 (MQ=255)
                                                                                                                                             tgagcGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCa                                          <  4:698513/110‑1 (MQ=255)
                                                                                                                                             tGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCa                                          >  3:698513/1‑114 (MQ=255)
                                                                                                                                             tGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTcc       >  4:55151/1‑149 (MQ=255)
                                                                                                                                              gACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGtt                                                        >  1:376240/1‑99 (MQ=255)
                                                                                                                                              gACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGtt                                                        <  2:376240/99‑1 (MQ=255)
                                                                                                                                                 cGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGGTCAGATCGCTACCCACGGCGGTTTCCGCAt   >  1:580250/1‑149 (MQ=255)
                                                                                                                                                  gTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATg  >  1:238418/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATG  >  NC_000913/2093821‑2094115

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