Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,271,8270AT47.3% 63.2 / 61.7 74V8E (GTA→GAA) chaAcalcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (2/37);  new base T (35/0);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 8.49e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACACTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCATTGG  >  NC_000913/1271677‑1271977
                                                                                                                                                       |                                                                                                                                                        
acaacaCTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTa                                                                                                                                                          <  2:1131623/151‑1 (MQ=255)
    cacTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTcccc                                                                                                                                                      >  2:604953/1‑149 (MQ=255)
      cTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGcc                                                                                                                                                     >  2:587697/1‑150 (MQ=255)
        aaaGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCctc                                                                                                                                                   >  2:313837/1‑150 (MQ=255)
          aGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTg                                                                                                                                                <  2:901789/151‑1 (MQ=255)
          aGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCtt                                                                                                                                                 >  1:638028/1‑150 (MQ=255)
          aGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTATCGCctcct                                                                                                                                                 >  2:959273/1‑148 (MQ=255)
               cTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCa                                                                                                                                            <  2:901527/150‑1 (MQ=255)
                tACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCGCCTCTTTAGCAtt                                                                                                                                          >  2:380180/1‑151 (MQ=255)
                    tAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCTCTTGAGCAttt                                                                                                                                         >  1:9116/1‑148 (MQ=255)
                     aaaaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACa                                                                                                                                     <  2:682108/151‑1 (MQ=255)
                     aaaaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCTCCTTGGCAttt                                                                                                                                         >  1:451632/1‑147 (MQ=255)
                      aaaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCCCCTCtttcggattt                                                                                                                                         >  2:37405/1‑138 (MQ=255)
                         tACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAGCATTTTTCAtt                                                                                                                                   >  2:225465/1‑149 (MQ=255)
                          aCCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTAt                                                                                                                                 <  1:650446/150‑1 (MQ=255)
                           ccAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGCTTTTTCCGCCCCTTGTGCCtttttcttttt                                                                                                                                 >  1:640881/1‑141 (MQ=255)
                             aaTAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCGCCTCTTTAGCATTTGAAATTATTGtt                                                                                                                             >  2:378477/1‑151 (MQ=255)
                              aTAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTa                                                                                                                            <  2:1090553/151‑1 (MQ=255)
                                aaGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGAAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTAt                                                                                                                           <  1:543890/150‑1 (MQ=255)
                                   cgcAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCCCCCCTTTTGCATTTTATAtttttggttttcc                                                                                                                        >  1:872275/1‑139 (MQ=255)
                                     cAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCtt                                                                                                                     <  1:396504/150‑1 (MQ=255)
                                     cAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCTCCCCTTTTGCATTTTTCATTTTTGTTTTTCCttt                                                                                                                    >  1:694643/1‑151 (MQ=255)
                                     cAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTTGCATTTTTCATTATTGtttttccctt                                                                                                                     >  2:1143948/1‑143 (MQ=255)
                                     cAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCTCTTGAGCATTTGAAATTATGGTTTTCCCctt                                                                                                                    >  2:787690/1‑151 (MQ=255)
                                      aagaagATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCttt                                                                                                                    <  2:221930/150‑1 (MQ=255)
                                      aagaagATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTg                                                                                                                   <  2:1081819/151‑1 (MQ=255)
                                         aagaTTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAg                                                                                                                <  2:403967/151‑1 (MQ=255)
                                         aagaTTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCTCTTTTGCATTTTACATTTTTGTTATTCCctt                                                                                                                    >  1:317943/1‑147 (MQ=255)
                                             ttGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAGAAATTTACATTTTTGtttttccctt                                                                                                                     >  2:566371/1‑135 (MQ=255)
                                              tGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGa                                                                                                            <  2:1003642/150‑1 (MQ=255)
                                                aTGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTAGGATTTTACATTTTTGTTTTTCCCTTTCAGGTGaa                                                                                                           >  2:527604/1‑149 (MQ=255)
                                                  ggCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAAttt                                                                                                        <  1:1089067/150‑1 (MQ=255)
                                                   gCAATGACCACTGGTGGTGTCTGGCTGCTTCCCCCCAGGAACAGCACTACCGGCGCAAGAACCGGGAAAATAAGCGAAGTGTCATTGTGGCGGGTTTTTAGCTCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAAttt                                                                                                        <  1:997366/149‑1 (MQ=255)
                                                          ccACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAAt                                                                                              <  1:390709/151‑1 (MQ=255)
                                                            aCTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGt                                                                                            <  1:1062370/151‑1 (MQ=255)
                                                            aCTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTTGCATTTTACATTATTGTTTTTCCTTTTCAGATGAATTTATTTaaaa                                                                                               >  2:1109418/1‑148 (MQ=255)
                                                             cTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCTCTTTAGCATTTTTTATTATTGTTTTTCCTTTGCAGATGAATTTATTGaaaa                                                                                               >  2:102078/1‑147 (MQ=255)
                                                              tGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCCTTAGCAATTTTCATTTTTGTTTTTCCTTTTCCGATGGAtttttttaaaa                                                                                                >  2:242660/1‑138 (MQ=255)
                                                               ggTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTTGCATTTGTCATTTTGGTTTTTCCTTTTCAGATGAAtttttttaaaa                                                                                                >  1:589456/1‑137 (MQ=255)
                                                                gTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTaaaa                                                                                        <  2:1066579/151‑1 (MQ=255)
                                                                 tAGTGTCTGGCTGCTTCCCCACAGGAACGGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTaaaa                                                                                        <  1:806497/150‑1 (MQ=255)
                                                                    tgtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGCGTTTTTTCGCCCCTTTTGGACTTTTCATTTTTGTTTTCCCTTTTCAGATGAATTTTTTGAAAATATaaaaaa                                                                                      >  1:894510/1‑149 (MQ=255)
                                                                         ggCTGCTTCCCCACAGGAACAGCGCTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAAtaggta                                                                                <  1:223973/150‑1 (MQ=255)
                                                                         ggCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAAtaggta                                                                                <  2:461831/150‑1 (MQ=255)
                                                                         ggCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAAtaggta                                                                                <  1:587692/150‑1 (MQ=255)
                                                                             gctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCCTTGGCATTTTTCATTTTTGTTTTTCCTTTTCAGATGAAATTTTTTAAAAATTAAAAAAAAGGGAGGaa                                                                            >  1:475684/1‑150 (MQ=255)
                                                                               ttCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGCTTTTTTCCCCCCTTTTGACTTTTACATTTTTGTTTTCTCCTTTTAAGGGAATTTTTCCGAAAAGTAAAAAAAAAGTAGGaaaa                                                                          >  1:345527/1‑150 (MQ=255)
                                                                                tCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGaaaaa                                                                         <  2:643358/150‑1 (MQ=255)
                                                                                tCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGaaaa                                                                          <  1:70910/149‑1 (MQ=255)
                                                                                    cacaGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTAGGATTTTTCATTATTGTTATTCCCTTTCAAGTGAAATTTTTGAAAAAGTAAAAAAAAAGGAGGaaaaaaaa                                                                      >  2:679193/1‑149 (MQ=255)
                                                                                      caGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCTCCTTAGCATTTTACATTTTGGGTTTTCCTTTTCCGATGAATTTTTTGTAAAATTAAAAAAAAAGGAAGaaaaaaaa                                                                       >  2:548596/1‑145 (MQ=255)
                                                                                           aCAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAgtg                                                             <  1:180087/151‑1 (MQ=255)
                                                                                           aCAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCTTTTGCATTTGACATTTTTGTTTTTCCTTTTCAGATGAATTTTTTGAAAATGTAAAAAAAAGGTAGGGAAAAAAAAAGaaa                                                                >  1:760055/1‑148 (MQ=255)
                                                                                            cAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAgtgt                                                            <  1:673035/151‑1 (MQ=255)
                                                                                                  tACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTACGCCCCTTGAGGAATTTTAATTTTGGTTTTTCCTTTGCCGATGGATTTTTCGAAAAAGTAAAAAAAAAGGAGGaaaaaaaaaaaaaa                                                                >  1:399792/1‑141 (MQ=255)
                                                                                                      aGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCg                                                  <  2:1054369/151‑1 (MQ=255)
                                                                                                      aGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAGGATTTTACATTTTTGTTTTTCCTTTTCTGATGAATTTTTTGAAAATGTAAAAAAAAAGGAGGAAAAAAAAAAAAAAAAGTCTgggt                                                      >  1:31248/1‑145 (MQ=255)
                                                                                                             gAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAGCATTTTACATTTTTGTTTTTCCCTTTCAGGTGAATTTTTTGAAAAAGTAAAAAAAAAGTAGGGAAAAAAAAAAAAAAAGGCTGGGTATTGGGAACa                                            >  2:411431/1‑150 (MQ=255)
                                                                                                                ccGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACAtttt                                        <  1:37405/151‑1 (MQ=255)
                                                                                                                 cGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCGCCCCTTTAGCATTTGACATTATTGTTATTCCTTTGCAGATGAATTTTTTTAAAAAGTAAAAAAAAAGTAGGAAAAAAAAAAAAAAGAGTCTGGATATTGGTAAAAattt                                        >  1:22625/1‑150 (MQ=255)
                                                                                                                 cGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCTTTAGCATTTTACATTATGGTTATCCCTTTTCAGATGAATTTATTGAAAATGTAAAAAATAGGGAGGAAAAAAAAAAAAAAGAGGCTGGGTATCGGGAAAAtttt                                        >  2:440847/1‑150 (MQ=255)
                                                                                                                  gggAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCGCCCCCTGTGCATTTGTAATTTTTGTTATTCCTTTTCCGGTGAATTTTTTGAAAATGTAAAAAATAAGTAGGGAAAATAAAAAAAAGTGGTTGGGTATTGGTAACAtttt                                        >  2:256205/1‑149 (MQ=255)
                                                                                                                   ggAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTACCCCCCTTGAGCATTTGACATTATGGGTTTCCCTTTTCAGATGAATTTATTGAAAATATAAAAAAAAGGTAGGAAAAATAAAAAAAAAGGTCTGGATATCGGTAAAAattt                                        >  2:157291/1‑148 (MQ=255)
                                                                                                                   ggAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACg                                     <  1:378474/151‑1 (MQ=255)
                                                                                                                     aaaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGaa                                   <  1:527599/151‑1 (MQ=255)
                                                                                                                        aTAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCTCTTTAGGATTTTACATTTTGGTTTTTCCTTTTCAGATGAATTTTTTGAAAAATTAAAAAAAAAGGAGGGAAAAAAAAAAAAAAAGTCTGGGTATCGGTAAAAATTTATGAAatt                                >  1:518777/1‑151 (MQ=255)
                                                                                                                          aaGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAAttttt                              <  1:1174453/151‑1 (MQ=255)
                                                                                                                             cGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTAc                            <  2:622473/150‑1 (MQ=255)
                                                                                                                                     ccTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCt                      <  1:1100841/148‑1 (MQ=255)
                                                                                                                                     ccTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTc                   <  2:1128517/151‑1 (MQ=255)
                                                                                                                                        tgtgGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTggag                  <  2:839074/149‑4 (MQ=255)
                                                                                                                                         gtgGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGAt                <  1:1001637/150‑1 (MQ=255)
                                                                                                                                         gtgGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTCTACCCCGCTGTCGAt                <  1:242658/150‑1 (MQ=255)
                                                                                                                                                   ttttACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTc       <  2:605252/150‑1 (MQ=255)
                                                                                                                                                    tttACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCa      <  1:987055/150‑1 (MQ=255)
                                                                                                                                                       aCCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCATTgg  <  1:853646/151‑1 (MQ=255)
                                                                                                                                                       aCCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCATTg   <  1:242679/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                        
ACAACACTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAACATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCATTGG  >  NC_000913/1271677‑1271977

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