Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,582,141 T→G N93K (AAT→AAG yhhZ → putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,582,1410TG98.1% 422.0 / ‑1.7 107N93K (AAT→AAGyhhZputative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (49/56);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 110 total aligned reads displayed.

TAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCG  >  NC_000913/3581996‑3582286
                                                                                                                                                 |                                                                                                                                                 
tAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAgaga                                                                                                                                                <  2:1274334/149‑1 (MQ=255)
tAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGAtt                                                                                                                                              >  1:732345/1‑151 (MQ=255)
tAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGAtt                                                                                                                                              >  1:1455/1‑151 (MQ=255)
tAAATAATATTAATAATACGGGGCTTGGTGCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGAt                                                                                                                                               <  2:737483/150‑1 (MQ=255)
 aaaTAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGAttt                                                                                                                                             <  1:881856/151‑1 (MQ=255)
 aaaTAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGAtt                                                                                                                                              <  1:1196645/150‑1 (MQ=255)
     aataTTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGt                                                                                                                                          <  2:867885/150‑1 (MQ=255)
      atatTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAg                                                                                                                                        <  2:507171/151‑1 (MQ=255)
        atTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGAt                                                                                                                                      <  2:895393/151‑1 (MQ=255)
        atTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGAt                                                                                                                                      <  2:78839/151‑1 (MQ=255)
          taataaTACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTg                                                                                                                                    >  1:276775/1‑151 (MQ=255)
           aataatACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTgg                                                                                                                                   <  1:95582/151‑1 (MQ=255)
           aataatACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTg                                                                                                                                    <  1:300464/150‑1 (MQ=255)
            ataataCGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGa                                                                                                                                  <  2:1026265/151‑1 (MQ=255)
             taataCGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGa                                                                                                                                  >  2:199901/1‑150 (MQ=255)
              aataCGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGaaa                                                                                                                                <  1:512240/151‑1 (MQ=255)
              aataCGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGaa                                                                                                                                 >  2:868480/1‑150 (MQ=255)
               ataCGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGaaaa                                                                                                                               <  2:759603/151‑1 (MQ=255)
                  cGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGta                                                                                                                            <  1:824801/151‑1 (MQ=255)
                  cGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGt                                                                                                                             >  1:587761/1‑150 (MQ=255)
                   ggggCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCGTTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtat                                                                                                                           >  2:1194541/1‑151 (MQ=255)
                    gggCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtatt                                                                                                                          <  2:872783/151‑1 (MQ=255)
                    gggCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtatt                                                                                                                          >  1:1015419/1‑151 (MQ=255)
                     ggCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtatta                                                                                                                         <  1:197120/151‑1 (MQ=255)
                      gCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtattat                                                                                                                        <  2:1018148/151‑1 (MQ=255)
                      gCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtattat                                                                                                                        <  2:861689/151‑1 (MQ=255)
                        ttGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtattatt                                                                                                                       >  2:976945/1‑150 (MQ=255)
                         tGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGtattatt                                                                                                                       <  1:3713/149‑1 (MQ=255)
                          ggTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATtata                                                                                                                    <  2:520206/151‑1 (MQ=255)
                           gTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATtatat                                                                                                                   >  1:1088920/1‑151 (MQ=255)
                           gTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATtatat                                                                                                                   <  2:847630/151‑1 (MQ=255)
                             tCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATtatata                                                                                                                  <  2:262873/150‑1 (MQ=255)
                              cacaGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATAc                                                                                                                 >  2:481289/1‑150 (MQ=255)
                              cacaGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACa                                                                                                                <  1:410511/151‑1 (MQ=255)
                               acaGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACaa                                                                                                               <  2:214447/151‑1 (MQ=255)
                                 aGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACt                                                                                                             <  2:386097/151‑1 (MQ=255)
                                 aGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACt                                                                                                             >  2:240615/1‑151 (MQ=255)
                                  gTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTa                                                                                                            <  2:44164/151‑1 (MQ=255)
                                  gTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTa                                                                                                            <  2:248788/151‑1 (MQ=255)
                                      cATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAagag                                                                                                        <  1:922405/151‑1 (MQ=255)
                                      cATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAagag                                                                                                        <  2:278279/151‑1 (MQ=255)
                                       aTGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAgg                                                                                                       <  2:719703/151‑1 (MQ=255)
                                           tataACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGAGCt                                                                                                   >  1:255172/1‑151 (MQ=255)
                                            ataACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCtt                                                                                                  >  1:316833/1‑151 (MQ=255)
                                            ataACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCtt                                                                                                  >  2:300912/1‑151 (MQ=255)
                                              aaCATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCtttt                                                                                                <  1:1260111/151‑1 (MQ=255)
                                               aCATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCttttt                                                                                               <  1:1021548/151‑1 (MQ=255)
                                                  ttttGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTAt                                                                                            >  2:334997/1‑151 (MQ=255)
                                                    ttGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCg                                                                                          >  2:1150118/1‑151 (MQ=255)
                                                      gTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGc                                                                                        >  1:646314/1‑151 (MQ=255)
                                                            taattGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtc                                                                                   >  2:436666/1‑150 (MQ=255)
                                                              attGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtcatc                                                                                >  2:1060434/1‑151 (MQ=255)
                                                                 gATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCAc                                                                              <  2:67235/150‑1 (MQ=255)
                                                                      aaGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCAcc                                                                             >  1:836614/1‑146 (MQ=255)
                                                                         cACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAtt                                                                     >  1:1170944/1‑151 (MQ=255)
                                                                         cACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAtt                                                                     >  2:1192823/1‑151 (MQ=255)
                                                                         cACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCAtt                                                                     <  2:88786/151‑1 (MQ=255)
                                                                                ttatttatTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAAcc                                                              >  1:1135360/1‑151 (MQ=255)
                                                                                     tattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACt                                                          <  2:1062516/150‑1 (MQ=255)
                                                                                      attaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCATCTgg                                                        <  1:440849/151‑1 (MQ=255)
                                                                                      attaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTgg                                                        >  1:120685/1‑151 (MQ=255)
                                                                                      attaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTgg                                                        >  1:643217/1‑151 (MQ=255)
                                                                                         aaTTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGAt                                                      <  1:206345/150‑1 (MQ=255)
                                                                                         aaTTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGAt                                                      <  2:1455/150‑1 (MQ=255)
                                                                                             ccATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAg                                                  <  1:937495/150‑1 (MQ=255)
                                                                                             ccATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAg                                                  <  1:1126310/150‑1 (MQ=255)
                                                                                              cATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGaa                                                <  2:377137/151‑1 (MQ=255)
                                                                                              cATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGa                                                 >  2:136771/1‑150 (MQ=255)
                                                                                                 tAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAAc                                              >  1:331275/1‑150 (MQ=255)
                                                                                                 tAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAAc                                              >  1:1292784/1‑150 (MQ=255)
                                                                                                  aaCAATAATGAACAATTATTTATGGGATTTTACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACaa                                            <  2:725391/151‑1 (MQ=255)
                                                                                                   aCAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAAt                                           <  2:1169882/151‑1 (MQ=255)
                                                                                                       taatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACt                                       >  1:772683/1‑151 (MQ=255)
                                                                                                        aatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACt                                       <  1:161156/150‑1 (MQ=255)
                                                                                                        aatGAACAATTAGTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATTCAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATTCAGAAACAATAACTa                                      <  1:1204327/151‑1 (MQ=255)
                                                                                                          tGAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAtt                                    >  1:221246/1‑151 (MQ=255)
                                                                                                           gAACAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTa                                   <  1:594825/151‑1 (MQ=255)
                                                                                                            aaCAATTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAACCCAACTGGATACAGAAACAATAACTAtt                                    >  2:725534/1‑149 (MQ=255)
                                                                                                               aaTTATTTATGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTa                               <  2:720663/151‑1 (MQ=255)
                                                                                                                  tatttatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATg                             >  1:74199/1‑150 (MQ=255)
                                                                                                                   atttatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGa                            >  1:787891/1‑150 (MQ=255)
                                                                                                                     ttatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAAtt                         >  2:782712/1‑151 (MQ=255)
                                                                                                                     ttatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAAt                          <  2:1297465/150‑1 (MQ=255)
                                                                                                                     ttatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAAt                          <  2:1279151/150‑1 (MQ=255)
                                                                                                                      tatGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAAttt                        >  1:501694/1‑151 (MQ=255)
                                                                                                                        tGGGATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTAt                      >  2:904679/1‑151 (MQ=255)
                                                                                                                         gggATTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTAt                      <  1:27604/150‑1 (MQ=255)
                                                                                                                         gggATTTGACTTCTATCGAATAAAGAGATTGGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATc                     <  2:734349/151‑1 (MQ=255)
                                                                                                                          ggATTTGACTTCTATCGAATAAAGAGATTTGGGAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATcc                    <  1:282246/151‑1 (MQ=255)
                                                                                                                            aTTTGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCctc                  <  2:203067/151‑1 (MQ=255)
                                                                                                                              ttGACTTCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTg                >  1:68071/1‑151 (MQ=255)
                                                                                                                                    tCTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAAc           >  1:1012997/1‑150 (MQ=255)
                                                                                                                                     cTATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACAt         >  2:106210/1‑151 (MQ=255)
                                                                                                                                      tATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATc        >  2:506604/1‑151 (MQ=255)
                                                                                                                                      tATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATc        >  1:159165/1‑151 (MQ=255)
                                                                                                                                      tATCGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATc        >  1:676281/1‑151 (MQ=255)
                                                                                                                                         cGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCtt      <  2:855481/150‑1 (MQ=255)
                                                                                                                                         cGAATAAAGAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCtt      <  1:220756/150‑1 (MQ=255)
                                                                                                                                           aataaagaGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATc   >  1:26986/1‑151 (MQ=255)
                                                                                                                                            ataaagaGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATAGAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCg  <  2:176514/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAATTAATTGATAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCG  >  NC_000913/3581996‑3582286

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