Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,397,8300TC42.6% 80.1 / 20.2 61I265V (ATT→GTT) uspEstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (29/6);  new base C (0/26);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACG  >  NC_000913/1397680‑1397978
                                                                                                                                                      |                                                                                                                                                    
cttcttCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAAt                                                                                                                                                      >  1:522392/1‑151 (MQ=255)
    ttcGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATcccc                                                                                                                                                  >  1:397612/1‑151 (MQ=255)
              aGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCATTCACCTGTTCCGCCGTGTTGCCGTGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCTCCGCCTGTAAGt                                                                                                                                        >  1:298376/1‑151 (MQ=255)
              aGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGt                                                                                                                                        >  2:1206341/1‑151 (MQ=255)
                ttcagggggTGTTTGATACTGGTCAGTTTTAATGACCAGCAGGTCGCAGCGAAGATGACCAATCCCCTGTCCCGCCGTTTTGCCGAGAAATGCTGCTGAAATCCCGGTGCGTCCTACCGGGCCCAGACCCAAAACCCCCGCCTGTAAGt                                                                                                                                        <  1:374549/143‑1 (MQ=255)
                ttCACCGGGTGTCTGATTGTGGTCAGGTTAAATAACCAGCGGGTCGCAGCGAAGATGATCAACCCCCTTTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGTGCGCCCTACCGTCCCAAGACCCACAACCCCCGCCTGTAAGTg                                                                                                                                       <  1:1172139/150‑1 (MQ=255)
                ttCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTg                                                                                                                                       >  2:1190274/1‑150 (MQ=255)
                   aaCGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTcc                                                                                                                                   >  1:968288/1‑151 (MQ=255)
                     cGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGc                                                                                                                                 >  2:1170497/1‑151 (MQ=255)
                       ggTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAAATGATCAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATCCCGGTGCCCCCTACCGTGCCCAAAACCACACTCCCCGCCTGTAAGTGCTCCGCCa                                                                                                                               <  1:1207815/151‑1 (MQ=255)
                           tCTGATACTGGTCAGTTTAAATGCCCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCTGCTGAAATCCCGGGGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCaa                                                                                                                              <  2:1170592/148‑1 (MQ=255)
                               aTACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAgg                                                                                                                        <  2:983949/150‑1 (MQ=255)
                                 aCTGGTCAGGTTTGATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAAt                                                                                                                     >  1:254012/1‑151 (MQ=255)
                                           tttATTGACAAGCGGGTCGCGGGAAAGAGAACAAATCCCTTTTCCCGCCGGGTTGCGAAGGATTGCTGCTAAAATCCCGGCGGGCCCGCCGGCCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtcttc                                                                                                            <  1:570713/150‑1 (MQ=255)
                                               aTGACCAGCGGGTCGAACCAAAGATGATCAACCCCCTGTCCCGCCGTTTGGCCGTGAAATGCCGCTAAAATACCGGTGCCCCCTCCCGGGCCCAAACCCAAAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCa                                                                                                       <  1:201836/151‑1 (MQ=255)
                                                 gACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                                                                        >  1:207169/1‑148 (MQ=255)
                                                   ccAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATCCCGGTGCGTCCTACCGTGCCCAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGAcc                                                                                                   <  2:334017/151‑1 (MQ=255)
                                                    cagcagGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCt                                                                                                  >  2:387286/1‑151 (MQ=255)
                                                      gcagGTCGCAGCGAAGATGATCAACCCCCTTTTCCGCCGTGTTGCCGAGGAATGCTGCTAAAATCCCGGGGCGCCCTACCGTCCCCAGAACCACATCCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCtt                                                                                                 <  1:707855/150‑1 (MQ=255)
                                                        agGTCGCAGGGAAGAAGTTCAATCCCCTTTTCCGCCGTGTTGCCGAGGATTGCTTCTAAAATCCCGTTGGTTCCTGCCGTGCCCAGACCCCCAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGa                                                                                                     <  1:980918/144‑1 (MQ=255)
                                                         ggTCGCAGCGAAGATGATCATTCCCCTGTCCCGCCGTGTTGCGGAGGAATGCTGCTGAAATCCCGGTGCGTCCTCCCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTc                                                                                             <  2:188550/151‑1 (MQ=255)
                                                         ggTCGCAGCGAAGAGGATCAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCtt                                                                                                 <  1:325897/147‑1 (MQ=255)
                                                          gccGCAGCGAAGATGACCAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATCCGGGTGCGTCCTACCGTCCCAAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCt                                                                                            <  1:88233/149‑1 (MQ=255)
                                                             gCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACg                                                                                         >  2:1136293/1‑151 (MQ=255)
                                                              cAGCGAAGATGATCAACCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATCCCGGTGCGTCCTCCCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACg                                                                                         <  1:260631/150‑1 (MQ=255)
                                                                  gAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCg                                                                                     >  2:688048/1‑150 (MQ=255)
                                                                   aaGATGATCAATCCCCTTTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATCCCGGTGCGCCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTg                                                                                       <  1:360975/147‑1 (MQ=255)
                                                                   aaGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGAAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAAAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGt                                                                                    <  2:269448/150‑1 (MQ=255)
                                                                        gatCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGGGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGt                                                                               <  2:1123691/150‑1 (MQ=255)
                                                                            aatcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTc                                                                           >  1:1017364/1‑150 (MQ=255)
                                                                                 ccTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAAt                                                                     >  2:924906/1‑151 (MQ=255)
                                                                                 ccTGTCCCGCCGTGTTGCCGAGAAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTa                                                                       <  1:417066/149‑1 (MQ=255)
                                                                                       ccgccgTTTTGCCGAGAAATGCTGCTAAAATACCGGTGCGCCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGa                                                                <  1:335173/150‑1 (MQ=255)
                                                                                       ccgccgTGTTGCCGAGGATTGCTGCTGAAATCCCGGTGCGCCCTCCCGTCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGaa                                                               <  2:216095/151‑1 (MQ=255)
                                                                                       ccgccgTGTTGCCGAGGAATGCTGCTAAAATCCGGGTGGGTCCTCCCGTCCCCAAACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGaa                                                               <  2:969929/151‑1 (MQ=255)
                                                                                       ccgccgTGTTGCCGAGAAATGCTGCTAAAATCCCGGTCCGTCCTACCGTCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCg                                                                 <  1:977029/149‑1 (MQ=255)
                                                                                        cgccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAt                                                              >  2:615211/1‑151 (MQ=255)
                                                                                         gccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAtt                                                             >  2:730362/1‑151 (MQ=255)
                                                                                         gccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAtt                                                             >  2:432116/1‑151 (MQ=255)
                                                                                          ccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAttt                                                            >  1:216738/1‑151 (MQ=255)
                                                                                                  ccGAGGAATGCTGCTGAAATCCCGGGGCGCCCTCCCGTCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGcc                                                       <  1:976172/148‑1 (MQ=255)
                                                                                                  ccGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGc                                                     >  1:1074428/1‑150 (MQ=255)
                                                                                                    gAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAgg                                                  >  1:1196459/1‑151 (MQ=255)
                                                                                                      ggAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGc                                                >  1:825408/1‑151 (MQ=255)
                                                                                                          tgctgcTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTc                                            >  1:794867/1‑151 (MQ=255)
                                                                                                            ctgctgAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCAt                                          >  2:767520/1‑151 (MQ=255)
                                                                                                             tgctgAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCAtt                                         >  2:1052161/1‑151 (MQ=255)
                                                                                                              gctgAAATCCCGGTGCGCCCTCCCGTCCCCAAACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTc                                            <  1:1210218/147‑1 (MQ=255)
                                                                                                                tgAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGcc                                      >  2:48827/1‑151 (MQ=255)
                                                                                                                  aaaTACCGGTGCGTCCTACCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTg                                        <  2:272871/147‑1 (MQ=255)
                                                                                                                     tACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCa                                  >  1:1044083/1‑150 (MQ=255)
                                                                                                                            gCGTCCTCCCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTg                           <  2:433627/150‑1 (MQ=255)
                                                                                                                                ccTACCGTGCCCAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCAc                      <  1:541597/151‑1 (MQ=255)
                                                                                                                                ccTACCGTCCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGccc                        <  1:745060/149‑1 (MQ=255)
                                                                                                                                         cccAGACCCACACTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCAt             <  2:205976/151‑1 (MQ=255)
                                                                                                                                         cccAGAACCCCAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATg                 <  1:1146385/147‑1 (MQ=255)
                                                                                                                                               aCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTtat       >  2:177990/1‑151 (MQ=255)
                                                                                                                                                ccACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTtata      <  2:23095/151‑1 (MQ=255)
                                                                                                                                                  acaATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATaa     >  2:406327/1‑150 (MQ=255)
                                                                                                                                                  acaACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATaaa    <  2:716983/151‑1 (MQ=255)
                                                                                                                                                    aaTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACg  >  2:1165379/1‑151 (MQ=255)
                                                                                                                                                    aaTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACg  >  1:721135/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACG  >  NC_000913/1397680‑1397978

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