Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,199,872 T→C Y365C (TAC→TGC)  ygiF ← inorganic triphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,199,8720TC97.9% 344.1 / ‑1.1 97Y365C (TAC→TGC) ygiFinorganic triphosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (43/52);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TCAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGTAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTACCAA  >  NC_000913/3199725‑3200017
                                                                                                                                                   |                                                                                                                                                 
tCAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGc                                                                                                                                                >  1:10931/1‑151 (MQ=255)
tCAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGc                                                                                                                                                >  1:706374/1‑151 (MQ=255)
 cAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCg                                                                                                                                               >  1:570638/1‑151 (MQ=255)
  aaTTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCg                                                                                                                                               <  2:1142955/150‑1 (MQ=255)
  aaTTTCGAGGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGa                                                                                                                                              <  1:446404/151‑1 (MQ=255)
     ttCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCg                                                                                                                                           >  1:426019/1‑151 (MQ=255)
       cGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCg                                                                                                                                           <  2:1286540/149‑1 (MQ=255)
        gATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGTAGCGATCGCCt                                                                                                                                        >  2:1611298/1‑151 (MQ=255)
        gATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCt                                                                                                                                        >  2:573391/1‑151 (MQ=255)
        gATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCt                                                                                                                                        >  2:1273520/1‑151 (MQ=255)
         aTGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCt                                                                                                                                        >  2:1123111/1‑150 (MQ=255)
             gTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACgg                                                                                                                                   >  1:70277/1‑151 (MQ=255)
               tGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGc                                                                                                                                  >  1:1412665/1‑150 (MQ=255)
                 cccGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTg                                                                                                                                <  2:930416/150‑1 (MQ=255)
                 cccGGCCGCAATAGCGGGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTgg                                                                                                                               <  1:678197/151‑1 (MQ=255)
                    ggTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGc                                                                                                                              >  2:630175/1‑149 (MQ=255)
                      tCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGa                                                                                                                           >  2:1011217/1‑150 (MQ=255)
                         cAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTACCCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGaaaa                                                                                                                        <  2:854303/150‑1 (MQ=255)
                         cAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAAc                                                                                                                       >  1:978004/1‑151 (MQ=255)
                         cAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGGATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGaaaa                                                                                                                        >  1:851883/1‑150 (MQ=255)
                             aGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCt                                                                                                                    <  2:110011/150‑1 (MQ=255)
                              gCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTACCCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCttt                                                                                                                  <  2:1526981/151‑1 (MQ=255)
                              gCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCttt                                                                                                                  >  1:1042500/1‑151 (MQ=255)
                                gTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttc                                                                                                                >  1:1037778/1‑151 (MQ=255)
                                  gATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttc                                                                                                                >  1:233796/1‑149 (MQ=255)
                                    tGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttcagtt                                                                                                            >  2:1654756/1‑151 (MQ=255)
                                    tGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttcagtt                                                                                                            >  2:1278120/1‑151 (MQ=255)
                                        gCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttcagttcag                                                                                                         >  2:1147585/1‑150 (MQ=255)
                                        gCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGc                                                                                                        >  2:1176389/1‑151 (MQ=255)
                                         ccccTGCCGGTTCCCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGc                                                                                                        <  1:172006/150‑1 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTttcagttca                                                                                                          <  1:430416/148‑1 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGc                                                                                                        <  2:1205136/150‑1 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGc                                                                                                        <  2:1438964/150‑1 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCg                                                                                                       >  1:527092/1‑151 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCg                                                                                                       <  2:1321118/151‑1 (MQ=255)
                                         ccccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCg                                                                                                       <  2:1689032/151‑1 (MQ=255)
                                           ccTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGc                                                                                                     <  1:902190/151‑1 (MQ=255)
                                            cTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGc                                                                                                     >  2:934881/1‑150 (MQ=255)
                                               ccAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTACCCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATg                                                                                                  <  2:173048/150‑1 (MQ=255)
                                                    tctcCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCggg                                                                                             >  1:847993/1‑150 (MQ=255)
                                                     ccacagcaagGCGTGCGCGCCGACAGGACCATATTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGTCGTGGCAACTGGGCACGGCAGCGATCGCCTAACGACTGGCAGAAAACTCTTTTCAGTCCAGCGGCATGCCggg                                                                                             <  1:70208/146‑1 (MQ=255)
                                                       ccagccagGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGaaa                                                                                          <  1:1353346/150‑1 (MQ=255)
                                                         agccagGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGGAAGa                                                                                        >  2:1064359/1‑150 (MQ=255)
                                                          gccagGCTTGCGCGACGCAAGGATAATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGaaa                                                                                          <  2:570689/147‑1 (MQ=255)
                                                           ccagGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGAt                                                                                       <  2:593394/149‑1 (MQ=255)
                                                           ccagGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATgg                                                                                     <  1:686455/151‑1 (MQ=255)
                                                                  tGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGGTATCgg                                                                              >  2:684322/1‑151 (MQ=255)
                                                                       gacgacAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGGAAGATGGAAATCGGCaa                                                                           >  1:704121/1‑149 (MQ=255)
                                                                          gacAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGGTAACCGCAAAGc                                                                        >  1:189911/1‑149 (MQ=255)
                                                                              ggATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGa                                                                  <  1:1326745/151‑1 (MQ=255)
                                                                               gATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGttt                                                                    >  2:888542/1‑148 (MQ=255)
                                                                                 tCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGGTATCGGCCAAGCGGTTTaa                                                                 >  1:933448/1‑149 (MQ=255)
                                                                                 tCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGaa                                                                 <  1:1358358/149‑1 (MQ=255)
                                                                                 tCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAAATGGATATCGGCAAAGCGGTTGaa                                                                 >  1:126540/1‑149 (MQ=255)
                                                                                  cATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAgg                                                               <  2:198075/150‑1 (MQ=255)
                                                                                   aTAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGa                                                              <  1:417781/150‑1 (MQ=255)
                                                                                     aGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGt                                                            <  1:500097/150‑1 (MQ=255)
                                                                                          cccGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTg                                                       <  2:1432783/150‑1 (MQ=255)
                                                                                           ccGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGAt                                                     <  1:340840/151‑1 (MQ=255)
                                                                                               cAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCTATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGAtt                                                    <  2:1495295/148‑1 (MQ=255)
                                                                                                aGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGAAGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTg                                                 >  2:1465098/1‑150 (MQ=255)
                                                                                                    aCAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGGCATCGGCAAAGCGGTTGGAGGAGTCGCTGATTTTGccc                                              >  1:488241/1‑149 (MQ=255)
                                                                                                      aGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTgcg                                          <  2:1256323/151‑1 (MQ=255)
                                                                                                       gTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTgcgc                                         <  2:1579120/151‑1 (MQ=255)
                                                                                                       gTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAACGTTTGAAGGAGTCGCTGATTTTGCCCTgcg                                          >  1:1363748/1‑150 (MQ=255)
                                                                                                         aTTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCtt                                       <  2:1576476/151‑1 (MQ=255)
                                                                                                           tGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCttt                                      <  2:967166/150‑1 (MQ=255)
                                                                                                            gAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTgg                                    <  2:735521/151‑1 (MQ=255)
                                                                                                            gAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTg                                     <  2:817251/150‑1 (MQ=255)
                                                                                                            gAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTg                                     <  1:463424/150‑1 (MQ=255)
                                                                                                              gTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCa                                  <  1:97448/151‑1 (MQ=255)
                                                                                                              gTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCa                                  <  2:890677/151‑1 (MQ=255)
                                                                                                                cAATATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATc                                <  1:355268/151‑1 (MQ=255)
                                                                                                                 aaTATCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCt                               <  2:933516/151‑1 (MQ=255)
                                                                                                                  atatCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCt                               >  2:1485929/1‑150 (MQ=255)
                                                                                                                  atatCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTa                              <  1:47848/151‑1 (MQ=255)
                                                                                                                   tatCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTa                              >  2:1197416/1‑150 (MQ=255)
                                                                                                                    atCACGCGTCAGGCGTGGCAACTGGTCACGGTAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                            <  1:1611200/151‑1 (MQ=255)
                                                                                                                    atCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaa                             <  2:1042579/150‑1 (MQ=255)
                                                                                                                     tCACGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                            >  2:877286/1‑150 (MQ=255)
                                                                                                                     tCACGCGTCAGGCGTGGCAACTGGTCACAGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAAGAGGCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                            >  1:242653/1‑150 (MQ=255)
                                                                                                                       aCGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGGGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaaaa                          >  2:34761/1‑150 (MQ=255)
                                                                                                                       aCGCGTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAAt                         <  2:1037856/151‑1 (MQ=255)
                                                                                                                         gcgTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATgg                       <  1:795681/151‑1 (MQ=255)
                                                                                                                          cgTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGc                      <  1:979861/151‑1 (MQ=255)
                                                                                                                          cgTCAGGCGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGGTTGAAGGAGTCGCTGATTTTTCCCTGCGCTTTGGCATCTaaaaa                          >  1:890763/1‑147 (MQ=255)
                                                                                                                                 cGTGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtg               <  1:1251540/151‑1 (MQ=255)
                                                                                                                                   tGGCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGAAGGGTATCGGCGAAGCGGTTGAAAGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAAAAGCTGcc                  >  1:1205237/1‑146 (MQ=255)
                                                                                                                                    ggCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgctt            <  1:1123025/151‑1 (MQ=255)
                                                                                                                                     gCAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgcttt           <  2:1228292/151‑1 (MQ=255)
                                                                                                                                      cAACTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgctttg          <  2:1595719/151‑1 (MQ=255)
                                                                                                                                       aaCTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTTCCCTGCGCTTTGGCATCTAAAAAAGGCTGCCAtgctttg          >  2:6121/1‑150 (MQ=255)
                                                                                                                                       aaCTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAAATGGATATCGGCAAAGCGTTTGAAAGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAAATGCTGCCAtgcttt           >  1:284964/1‑149 (MQ=255)
                                                                                                                                        aCTGGTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgctttgct        <  2:1047492/151‑1 (MQ=255)
                                                                                                                                           ggTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgctttgctt       <  2:1167498/149‑1 (MQ=255)
                                                                                                                                           ggTCACGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTac     <  1:630117/151‑1 (MQ=255)
                                                                                                                                               aCGGCAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTTCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTaccaa  >  2:1416327/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGACAGGATCATAGTAACCCGCCAGCAACAGTATTGAGTCAATATCACGCGTCAGGCGTGGCAACTGGTCACGGTAGCGATCGCCTAACGGCTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTACCAA  >  NC_000913/3199725‑3200017

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