Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 78,831 T→A 10.7% intergenic (+32/+17) setA → / ← leuD broad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,8310TA10.7% 538.0 / 45.3 168intergenic (+32/+17)setA/leuDbroad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit
Reads supporting (aligned to +/- strand):  ref base T (63/87);  new base A (8/10);  total (71/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

Only 100 of 172 total aligned reads displayed.

GCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGG  >  NC_000913/78683‑78978
                                                                                                                                                    |                                                                                                                                                   
gctggcGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                     <  1:1230658/149‑1 (MQ=255)
       ttATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCgg                                                                                                                                              <  4:1255070/149‑1 (MQ=255)
          ttCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGct                                                                                                                                           <  4:624586/149‑2 (MQ=255)
            cAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGcttt                                                                                                                                         <  4:809057/149‑4 (MQ=255)
             aGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGctt                                                                                                                                          >  1:1171488/1‑145 (MQ=255)
              gggAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAggg                                                                                                                                       <  3:224952/149‑1 (MQ=255)
                  gCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAtt                                                                                                                                   >  4:142455/1‑149 (MQ=255)
                     aTTGCACAAAGTTGGGGGCATTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGTTAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                                <  2:553801/149‑1 (MQ=255)
                     aTTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                                <  1:1310854/149‑1 (MQ=255)
                     aTTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                                <  4:103153/149‑1 (MQ=255)
                         cacaAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCAt                                                                                                                            >  3:864978/1‑149 (MQ=255)
                         cacaAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCAt                                                                                                                            <  3:908189/149‑1 (MQ=255)
                          acaAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATa                                                                                                                           >  3:617076/1‑149 (MQ=255)
                          acaAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATa                                                                                                                           >  3:683439/1‑149 (MQ=255)
                            aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGa                                                                                                                                                  >  1:1452712/1‑124 (MQ=255)
                            aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGa                                                                                                                                                  <  2:1452712/124‑1 (MQ=255)
                            aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaaa                                                                                                                         >  2:685231/1‑149 (MQ=255)
                            aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAACTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaaa                                                                                                                         <  3:1088889/149‑1 (MQ=255)
                               gTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGc                                                                                                                      >  1:1563560/1‑149 (MQ=255)
                                ttGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTATTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGATGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAGTCATAACCGCa                                                                                                                     <  4:910828/149‑1 (MQ=255)
                                   ggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgt                                                                                                                  <  4:804878/149‑1 (MQ=255)
                                   ggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgt                                                                                                                  <  2:1584182/149‑1 (MQ=255)
                                    gggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgtt                                                                                                                 <  1:475873/149‑1 (MQ=255)
                                      gCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTTTTTTTAACCGCATAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgttgt                                                                                                               >  4:1495032/1‑149 (MQ=255)
                                        aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                             >  4:571794/1‑149 (MQ=255)
                                        aCTTTGCTGTCTACTGGGTAATCGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                             >  1:1147914/1‑149 (MQ=255)
                                          tttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCGAGGGGATTAATTCATAAACGCAGGTTGTTTTg                                                                                                           >  4:1422793/1‑149 (MQ=255)
                                          tttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTg                                                                                                           <  3:1434900/149‑1 (MQ=255)
                                           ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGc                                                                                                          <  3:571794/149‑1 (MQ=255)
                                           ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGc                                                                                                          <  1:703116/149‑1 (MQ=255)
                                            tGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                         <  1:1057907/149‑1 (MQ=255)
                                             gCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCtt                                                                                                        <  4:1241905/149‑1 (MQ=255)
                                              cTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTc                                                                                                       >  2:1160231/1‑149 (MQ=255)
                                               tGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCa                                                                                                      >  2:1316086/1‑149 (MQ=255)
                                                gTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCAt                                                                                                     <  2:1563560/149‑1 (MQ=255)
                                                 tCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATa                                                                                                    >  3:33248/1‑149 (MQ=255)
                                                   tACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                             <  4:594932/138‑1 (MQ=255)
                                                    aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGc                                                                                                 <  2:873445/149‑1 (MQ=255)
                                                    aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGc                                                                                                 <  3:167383/149‑1 (MQ=255)
                                                        ggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCa                                                                                             >  1:254678/1‑149 (MQ=255)
                                                        ggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCa                                                                                             >  1:769777/1‑149 (MQ=255)
                                                         gTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCaa                                                                                            <  2:307235/149‑1 (MQ=255)
                                                          tAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAAt                                                                                           >  4:266463/1‑149 (MQ=255)
                                                          tAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                               >  4:1104628/1‑97 (MQ=255)
                                                           aaTTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATg                                                                                          <  3:1082233/149‑1 (MQ=255)
                                                             ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGc                                                                                        <  2:400370/149‑1 (MQ=255)
                                                                cGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtc                                                                                     <  4:683439/149‑1 (MQ=255)
                                                                cGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtc                                                                                     <  3:1509371/149‑1 (MQ=255)
                                                                   ttATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCCGCAATGGcgtcgtc                                                                                  <  2:546590/149‑1 (MQ=255)
                                                                     aTTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtcgt                                                                                <  4:36654/149‑1 (MQ=255)
                                                                      tttCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtg                                                                               <  3:267512/149‑1 (MQ=255)
                                                                        tCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtgct                                                                             <  3:1202192/149‑1 (MQ=255)
                                                                        tCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtgct                                                                             <  2:490621/149‑1 (MQ=255)
                                                                            ttgtCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtgc                                                                              >  2:1121958/1‑144 (MQ=255)
                                                                            ttgtCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCaa                                                                         >  4:863234/1‑149 (MQ=255)
                                                                               tCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAt                                                                                                                                    <  3:760148/87‑1 (MQ=255)
                                                                                 gCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAACGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTaa                                                                    >  4:1322547/1‑149 (MQ=255)
                                                                                  cATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAg                                                                   <  2:254678/149‑1 (MQ=255)
                                                                                   attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                  <  2:1229219/149‑1 (MQ=255)
                                                                                   attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                  >  1:99007/1‑149 (MQ=255)
                                                                                    ttattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcg                                                                                    <  1:1354260/130‑1 (MQ=255)
                                                                                    ttattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGcc                                                                 <  1:1316086/149‑1 (MQ=255)
                                                                                    ttattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGcttt                                                                                                                                         >  1:1106994/1‑74 (MQ=39)
                                                                                    ttattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGcttt                                                                                                                                         <  2:1106994/77‑4 (MQ=38)
                                                                                              ccGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTc                                                       <  1:1489417/149‑1 (MQ=255)
                                                                                                 cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAg                                                    >  3:1238985/1‑149 (MQ=255)
                                                                                                 cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAACCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAg                                                    >  3:707917/1‑149 (MQ=255)
                                                                                                  aaaGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGa                                                   >  3:664116/1‑149 (MQ=255)
                                                                                                   aaGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAc                                                  <  3:320563/149‑1 (MQ=255)
                                                                                                       ttAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGtt                                              <  2:1392131/149‑1 (MQ=255)
                                                                                                       ttAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGtt                                              >  3:113450/1‑149 (MQ=255)
                                                                                                         aaaGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtca                                            <  3:509222/149‑1 (MQ=255)
                                                                                                         aaaGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtca                                            >  1:267424/1‑149 (MQ=255)
                                                                                                          aaGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                     >  1:800655/1‑43 (MQ=255)
                                                                                                          aaGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                     <  2:800655/43‑1 (MQ=255)
                                                                                                          aaGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtcat                                           >  1:632965/1‑149 (MQ=255)
                                                                                                           aGACGTTTGATGACGTGGACGATAGCGGAAATCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtcatc                                          <  1:1483010/149‑1 (MQ=255)
                                                                                                           aGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtc                                             <  1:903577/146‑1 (MQ=255)
                                                                                                           aGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtc                                             >  2:903577/1‑146 (MQ=255)
                                                                                                                tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCa                                     <  3:424748/149‑1 (MQ=255)
                                                                                                                   gatgaCGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTg                                  >  3:426029/1‑149 (MQ=255)
                                                                                                                     tgaCGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTggc                                <  1:1172301/149‑1 (MQ=255)
                                                                                                                             aCGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAg                        <  4:793455/149‑1 (MQ=255)
                                                                                                                             aCGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAg                        <  2:1539520/149‑1 (MQ=255)
                                                                                                                                   gCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCg                  >  3:569289/1‑149 (MQ=255)
                                                                                                                                    cGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGa                 <  4:617076/149‑1 (MQ=255)
                                                                                                                                       cttgCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCa                                     >  4:485844/4‑126 (MQ=255)
                                                                                                                                       cttgCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCa                                     <  3:485844/123‑1 (MQ=255)
                                                                                                                                       cttgCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATg               <  3:1275265/145‑1 (MQ=255)
                                                                                                                                       cttgCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATg               >  4:1275265/4‑148 (MQ=255)
                                                                                                                                        ttgCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg              >  1:531372/3‑148 (MQ=255)
                                                                                                                                          gCCCGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGCCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTaa           >  3:44761/1‑149 (MQ=255)
                                                                                                                                           cccGGTCATTTTACCGGGCAAGGGGATTAATTCATAATCGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTaaa          <  3:1020435/149‑1 (MQ=255)
                                                                                                                                            ccGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAg         >  4:1045603/1‑149 (MQ=255)
                                                                                                                                            ccGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAg         >  4:24048/1‑149 (MQ=255)
                                                                                                                                             cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGc        >  1:1396977/1‑149 (MQ=255)
                                                                                                                                              ggTCATTTGACCGGGCAAGGGGATTAATTCATAAACGc                                                                                                                      <  1:311553/38‑1 (MQ=255)
                                                                                                                                               gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGa      >  3:997919/1‑149 (MQ=255)
                                                                                                                                                tCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGAt     <  4:707917/149‑1 (MQ=255)
                                                                                                                                                  aTTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAg   >  3:862322/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGG  >  NC_000913/78683‑78978

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