Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,954,273 G→T 5.8% E266* (GAA→TAA) ‡ ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,954,2730GT 5.8% 594.9 / 19.5 173E266* (GAA→TAA) ‡ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (72/91);  new base T (6/4);  total (78/95)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 176 total aligned reads displayed.

GGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954125‑3954421
                                                                                                                                                     |                                                                                                                                                    
gggCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTg                                                                                                                                                      >  3:847187/1‑149 (MQ=255)
     aCGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAAc                                                                                                                                                 <  3:1299130/149‑1 (MQ=255)
       ggCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGc                                                                                                                                               <  4:845874/149‑1 (MQ=255)
             ctgctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGAc                                                                                                                                         <  2:678070/149‑1 (MQ=255)
                ctgGGAACCCACGGCCGACCTAAAGCACTGTTCCTTAAAGCTGGTAAACCGATTTTTGGATTGGACCAACCTTAATACGAGGAAAACGACGAAAGTGGACTGGCGCCGAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCt                                                                                                                                      <  4:257476/149‑1 (MQ=255)
                ctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCt                                                                                                                                      <  1:547656/149‑1 (MQ=255)
                 tgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTg                                                                                                                                     <  2:177223/149‑1 (MQ=255)
                     aaCCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGata                                                                                                                                 >  1:862898/1‑149 (MQ=255)
                      aCCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGatat                                                                                                                                >  4:1432555/1‑149 (MQ=255)
                        ccACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATcg                                                                                                                              >  2:1174016/1‑149 (MQ=255)
                           cGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGa                                                                                                                           <  2:593169/149‑1 (MQ=255)
                             gCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGAt                                                                                                                          <  2:535671/149‑1 (MQ=255)
                               cGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATgg                                                                                                                        <  2:1268855/149‑1 (MQ=255)
                                  ccGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTg                                                                                                                     <  1:1013943/149‑1 (MQ=255)
                                         cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  1:308040/149‑1 (MQ=255)
                                          agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             <  2:643470/149‑1 (MQ=255)
                                            cTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTaa                                                                                                           >  1:1138628/1‑149 (MQ=255)
                                              gTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAaca                                                                                                         <  4:282173/149‑1 (MQ=255)
                                                tCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACAcc                                                                                                       <  3:639106/149‑1 (MQ=255)
                                                  cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCATTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  2:1264215/1‑149 (MQ=255)
                                                  cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  1:805599/1‑149 (MQ=255)
                                                  cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  2:375107/1‑149 (MQ=255)
                                                  cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  3:1122562/1‑149 (MQ=255)
                                                   ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    >  4:1249372/1‑149 (MQ=255)
                                                   ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    <  2:1461202/149‑1 (MQ=255)
                                                   ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    <  3:601008/149‑1 (MQ=255)
                                                    tAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTAc                                                                                                   <  4:737080/149‑1 (MQ=255)
                                                         cTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAc                                                                                              <  3:437620/149‑1 (MQ=255)
                                                          tGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                                                                             <  4:621809/149‑1 (MQ=255)
                                                            gTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCtcg                                                                                                                                    <  4:1001511/108‑3 (MQ=255)
                                                            gTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTgg                                                                                                                                    >  3:1001511/1‑108 (MQ=255)
                                                                aCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  4:806521/126‑1 (MQ=255)
                                                                 cGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggc                                                                                      >  1:846741/1‑149 (MQ=255)
                                                                       gttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcggc                                                                                >  1:17316/1‑149 (MQ=255)
                                                                        ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              >  2:523155/1‑118 (MQ=255)
                                                                             ttGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    <  1:377464/123‑1 (MQ=255)
                                                                               gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAg                                                                        >  3:138776/1‑149 (MQ=255)
                                                                                 ccAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGct                                                                      >  4:520809/1‑148 (MQ=255)
                                                                                 ccAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      <  3:583578/149‑1 (MQ=255)
                                                                                  cAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCgg                                                                     >  2:162447/1‑149 (MQ=255)
                                                                                  cAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCgg                                                                     >  2:887658/1‑149 (MQ=255)
                                                                                  cAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCgg                                                                     >  4:441610/1‑149 (MQ=255)
                                                                                   aaaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGa                                                                    <  2:467159/149‑1 (MQ=255)
                                                                                    aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaa                                                                   >  2:1044857/1‑149 (MQ=255)
                                                                                    aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaa                                                                   >  4:354791/1‑149 (MQ=255)
                                                                                     aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaaa                                                                  <  1:984298/149‑1 (MQ=255)
                                                                                        ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCg                                                               >  2:1395803/1‑149 (MQ=255)
                                                                                        ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCg                                                               <  1:1044857/149‑1 (MQ=255)
                                                                                         tattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTTCTGGCGGCGGCGCAGGAAGCGGAAATCGa                                                              <  3:1249372/149‑1 (MQ=255)
                                                                                          attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCTAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:596207/149‑1 (MQ=255)
                                                                                          attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:536720/149‑1 (MQ=255)
                                                                                          attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:165044/149‑1 (MQ=255)
                                                                                          attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:702340/149‑1 (MQ=255)
                                                                                           ttaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATgg                                                                                                                        <  3:548470/89‑1 (MQ=255)
                                                                                           ttaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACt                                                            <  4:1447598/149‑1 (MQ=255)
                                                                                              cGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCa                                                         <  2:89106/149‑1 (MQ=255)
                                                                                               gAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGTAAGCGGAAATCGACTTCAc                                                        >  2:1060319/1‑149 (MQ=255)
                                                                                                    aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   >  3:729932/1‑149 (MQ=255)
                                                                                                    aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   <  2:846741/149‑1 (MQ=255)
                                                                                                    aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   <  1:1394178/149‑1 (MQ=255)
                                                                                                    aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   <  3:822035/149‑1 (MQ=255)
                                                                                                     aaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAg                                                  >  2:1362938/1‑149 (MQ=255)
                                                                                                      acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGa                                                                                                                   <  1:101830/83‑1 (MQ=255)
                                                                                                      acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGa                                                                                                                   >  2:101830/1‑83 (MQ=255)
                                                                                                          cgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGata                                             <  3:49282/149‑1 (MQ=255)
                                                                                                           gaAAGTGCAGTGCGGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGatat                                            <  2:1322402/149‑1 (MQ=255)
                                                                                                           gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGatat                                            >  1:1146324/1‑149 (MQ=255)
                                                                                                               gTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGTAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGAt                                        <  3:1150298/149‑1 (MQ=255)
                                                                                                                tGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaaa                                                                  <  2:1110741/122‑1 (MQ=255)
                                                                                                                tGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATa                                       >  1:1181158/1‑149 (MQ=255)
                                                                                                                 gCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATaa                                      <  1:1104514/149‑1 (MQ=255)
                                                                                                                  cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGa                                                                                                                                   >  2:444657/1‑55 (MQ=255)
                                                                                                                       ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTc                                <  3:354791/149‑1 (MQ=255)
                                                                                                                        cgcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTcc                               <  1:420630/149‑1 (MQ=255)
                                                                                                                          cgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCg                             >  1:492679/1‑149 (MQ=255)
                                                                                                                            tAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCa                           <  4:729932/149‑1 (MQ=255)
                                                                                                                             aaTATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCaa                          >  2:709119/1‑149 (MQ=255)
                                                                                                                               tatCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGata                                             <  1:385566/128‑1 (MQ=255)
                                                                                                                               tatCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGata                                             >  2:385566/1‑128 (MQ=255)
                                                                                                                                    ccAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTgg                                                                                       <  3:461097/81‑1 (MQ=255)
                                                                                                                                    ccAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCca                   >  2:991682/1‑149 (MQ=255)
                                                                                                                                       gTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAg                >  3:1100656/1‑149 (MQ=255)
                                                                                                                                       gTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAg                <  2:862898/149‑1 (MQ=255)
                                                                                                                                       atAAGGCGGCGTTCGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATGGATAAGCTTTCCCGCAAGGTTACGCAg                >  4:1262400/2‑149 (MQ=255)
                                                                                                                                         aaGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCt              <  4:911468/149‑1 (MQ=255)
                                                                                                                                         aaGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCt              <  3:1432555/149‑1 (MQ=255)
                                                                                                                                         aaGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCt              <  2:406449/149‑1 (MQ=255)
                                                                                                                                          tcatggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtg             >  2:737220/5‑149 (MQ=255)
                                                                                                                                           ggcggcGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgt            <  2:240129/149‑1 (MQ=255)
                                                                                                                                            atggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      >  4:1208958/3‑90 (MQ=255)
                                                                                                                                             tggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATaa                                      <  3:1090078/120‑1 (MQ=255)
                                                                                                                                             tggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAg                         >  3:1409909/2‑134 (MQ=255)
                                                                                                                                             tggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTgtgt          >  3:1426264/2‑149 (MQ=255)
                                                                                                                                             tggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTgtgt          >  3:763568/2‑149 (MQ=255)
                                                                                                                                              ggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   <  4:1380784/107‑1 (MQ=255)
                                                                                                                                              ggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   >  3:1380784/1‑107 (MQ=255)
                                                                                                                                               gcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTaa        <  4:1133716/149‑1 (MQ=255)
                                                                                                                                                 gTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      <  2:325432/85‑1 (MQ=255)
                                                                                                                                                 gTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      >  1:325432/1‑85 (MQ=255)
                                                                                                                                                   ttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGtt    <  2:544104/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954125‑3954421

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