Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,644 A→C 8.1% intergenic (‑50/+33) yfbU ← / ← yfbV UPF0304 family protein YfbU/UPF0208 membrane protein YfbV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,412,6440AC 8.1% 369.2 / 20.9 124intergenic (‑50/+33)yfbU/yfbVUPF0304 family protein YfbU/UPF0208 membrane protein YfbV
Reads supporting (aligned to +/- strand):  ref base A (53/61);  new base C (4/6);  total (57/67)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 126 total aligned reads displayed.

TTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTG  >  NC_000913/2412495‑2412793
                                                                                                                                                      |                                                                                                                                                     
tttGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTa                                                                                                                                                       <  1:214692/150‑1 (MQ=255)
  tGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATg                                                                                                                                                     <  1:1475969/150‑1 (MQ=255)
    cAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGcc                                                                                                                                                   <  2:437904/150‑1 (MQ=255)
    cAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGcc                                                                                                                                                   <  2:592405/150‑1 (MQ=255)
    cAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGcc                                                                                                                                                   >  2:598753/1‑150 (MQ=255)
      ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                                 >  1:404225/1‑150 (MQ=255)
      ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTAAGGTGAAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                                 >  2:884632/1‑150 (MQ=255)
       gCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttc                                                                                                                                                <  2:1279284/150‑1 (MQ=255)
        cGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttct                                                                                                                                               <  1:527702/150‑1 (MQ=255)
         gACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                                                              <  2:233258/150‑1 (MQ=255)
         gACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                                                              <  2:227311/150‑1 (MQ=255)
           cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCtttt                                                                                                                                            <  2:22830/150‑1 (MQ=255)
             gTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtat                                                                                                                                          >  2:439000/1‑150 (MQ=255)
              tAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtata                                                                                                                                         <  1:1335086/150‑1 (MQ=255)
              tAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttttttt                                                                                                                                         <  1:1007297/150‑4 (MQ=255)
              tAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttttttt                                                                                                                                         <  2:511105/150‑4 (MQ=255)
                  gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                     >  1:936537/1‑150 (MQ=255)
                  gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                     <  1:1380029/150‑1 (MQ=255)
                        cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCa                                                                                                                               >  2:124835/1‑150 (MQ=255)
                          tttGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgg                                                                                                                             >  1:256507/1‑150 (MQ=255)
                          tttGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgg                                                                                                                             <  1:598753/150‑1 (MQ=255)
                             gCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                                                          >  2:1440774/1‑150 (MQ=255)
                             gCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                                                          >  2:1111156/1‑150 (MQ=255)
                               cGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCa                                                                                                                        >  1:827168/1‑150 (MQ=255)
                                ggATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCaa                                                                                                                       >  1:31793/1‑150 (MQ=255)
                                  aTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATATCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAtt                                                                                                                     >  2:1506486/1‑150 (MQ=255)
                                   tCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTa                                                                                                                    >  1:284436/1‑150 (MQ=255)
                                      aGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaa                                                                                                                 <  1:1423191/150‑1 (MQ=255)
                                      aGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaa                                                                                                                 <  1:684345/150‑1 (MQ=255)
                                       gCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                                                                <  2:108083/150‑1 (MQ=255)
                                         aTAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAatc                                                                                                              >  2:460821/1‑150 (MQ=255)
                                          tAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAatca                                                                                                             >  1:1355005/1‑150 (MQ=255)
                                              catcatCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                                                         <  1:596577/150‑1 (MQ=255)
                                               atcatcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAg                                                                                                        <  2:827168/150‑1 (MQ=255)
                                                tcatcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAgg                                                                                                       <  1:1018758/150‑1 (MQ=255)
                                                     ttGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGt                                                                                                  >  1:83992/1‑150 (MQ=255)
                                                        tACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                             >  2:55459/1‑104 (MQ=255)
                                                         aCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                              >  2:1182325/1‑150 (MQ=255)
                                                            gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                           >  2:714155/1‑150 (MQ=255)
                                                            gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                           >  2:1191265/1‑150 (MQ=255)
                                                            gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                           >  2:1594267/1‑150 (MQ=255)
                                                                 gACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttg                                                                                      <  1:460821/150‑1 (MQ=255)
                                                                 gACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttg                                                                                      >  1:1234213/1‑150 (MQ=255)
                                                                   cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                    >  2:1599613/1‑150 (MQ=255)
                                                                   cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                    >  2:810848/1‑150 (MQ=255)
                                                                         aagcccgtGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTCATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAgg                                                                              <  2:349645/143‑1 (MQ=255)
                                                                         aaGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAgg                                                                              <  2:242179/150‑1 (MQ=255)
                                                                             cGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGgcgcg                                                                          <  2:1234213/150‑1 (MQ=255)
                                                                             cGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGgcgcg                                                                          >  1:746182/1‑150 (MQ=255)
                                                                                tGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCtt                                                                       >  1:515926/1‑150 (MQ=255)
                                                                                  cgcgTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTa                                                                      >  2:1419672/1‑150 (MQ=255)
                                                                                        ggTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCAGAGTATGCCTACTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTAACCAGTTGTTTGAAGGCGCGCTTAAGCgtgt                                                               >  2:1233720/1‑150 (MQ=255)
                                                                                          tCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCa                                                             <  1:1594267/150‑1 (MQ=255)
                                                                                             tttCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCt                                                          <  1:439000/150‑1 (MQ=255)
                                                                                                ccATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTg                                                                                  <  1:1263864/123‑1 (MQ=255)
                                                                                                ccATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTg                                                                                  <  1:1269793/123‑1 (MQ=255)
                                                                                                ccATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTg                                                                                  >  2:1269793/1‑123 (MQ=255)
                                                                                                ccATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAAt                                                       <  1:989126/150‑1 (MQ=255)
                                                                                                ccATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAAt                                                       <  1:1152532/150‑1 (MQ=255)
                                                                                                   ttgttgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGGCAGCTAATGcc                                                    <  2:1511953/150‑1 (MQ=255)
                                                                                                       tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGt                                                <  1:1208711/150‑1 (MQ=255)
                                                                                                       tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGt                                                >  2:126032/1‑150 (MQ=255)
                                                                                                        gACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTa                                               >  1:995439/1‑150 (MQ=255)
                                                                                                             cTGTATCACTCTACTACGGGGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAg                                          <  2:1460793/150‑1 (MQ=255)
                                                                                                              tGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAgg                                         <  2:404225/150‑1 (MQ=255)
                                                                                                                tATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCt                                       >  1:1035540/1‑150 (MQ=255)
                                                                                                                    aCTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGcc                                   <  2:515926/150‑1 (MQ=255)
                                                                                                                    aCTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGATTACCAGGcattag                                   <  1:1359203/150‑6 (MQ=255)
                                                                                                                    aCTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGGAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGcc                                   >  1:716667/1‑150 (MQ=255)
                                                                                                                     ctctACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCt                                  <  1:810848/150‑1 (MQ=255)
                                                                                                                         actacGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCaa                              >  2:1553263/1‑150 (MQ=255)
                                                                                                                             cGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                    >  2:1585267/1‑92 (MQ=255)
                                                                                                                             cGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                    <  1:1585267/92‑1 (MQ=255)
                                                                                                                               gTGAAAAAAAAGAAGTCTTAGTATGCCTTCTTTAATATGAGTAATCAGGGGTCAATTACAAATCATGAAGGAAAGTTTTATCCAGTTGTTTGAAGGAGAGCTTAAGATTGTCAGCTAATTCCTGTTAATCAGGCTTGACTTAAACGtgtg                        <  2:1599703/150‑1 (MQ=255)
                                                                                                                               gTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTg                        <  2:1202872/150‑1 (MQ=255)
                                                                                                                               gTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTg                        <  1:126032/150‑1 (MQ=255)
                                                                                                                               gTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTg                        <  2:256507/150‑1 (MQ=255)
                                                                                                                                tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGc                       <  1:629355/150‑1 (MQ=255)
                                                                                                                                tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGc                       <  2:284436/150‑1 (MQ=255)
                                                                                                                                 gAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc                      >  1:572854/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGAAGGCTGAGTATGCGTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGATAGTTTTATGCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCaa                    >  1:602088/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATc                                            <  1:553671/126‑1 (MQ=255)
                                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCaa                    >  2:20444/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCaa                    >  1:826680/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCaa                    >  1:1091565/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAac                   >  2:1108108/1‑150 (MQ=255)
                                                                                                                                      taaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAaca                  >  1:222764/2‑150 (MQ=255)
                                                                                                                                       aaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAacac                 >  2:1183411/1‑150 (MQ=255)
                                                                                                                                        aaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc                <  2:936537/150‑1 (MQ=255)
                                                                                                                                        aaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc                <  1:1191265/150‑1 (MQ=255)
                                                                                                                                        aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc                >  1:287067/1‑150 (MQ=255)
                                                                                                                                             agGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAaca                  <  2:498301/143‑1 (MQ=255)
                                                                                                                                             agGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAaca                  >  1:498301/1‑143 (MQ=255)
                                                                                                                                             agGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTcc           <  2:746182/150‑1 (MQ=255)
                                                                                                                                              ggCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa          <  1:1252583/150‑1 (MQ=255)
                                                                                                                                               gCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAg         <  1:1553263/150‑1 (MQ=255)
                                                                                                                                                  gAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACt      >  2:790071/1‑150 (MQ=255)
                                                                                                                                                   aGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTc     >  1:1438887/1‑150 (MQ=255)
                                                                                                                                                    gTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTcc    <  2:1035540/150‑1 (MQ=255)
                                                                                                                                                    gTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTcc    <  2:803966/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTG  >  NC_000913/2412495‑2412793

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