Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT45.9% 78.3 / 33.4 61N156K (AAC→AAAyiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (24/9);  new base T (0/28);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 7.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCC  >  NC_000913/4111597‑4111892
                                                                                                                                                        |                                                                                                                                                  
caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                                    >  2:432625/1‑151 (MQ=255)
 aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                   <  2:508271/151‑1 (MQ=255)
        ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTc                                                                                                                                             <  2:491265/149‑1 (MQ=255)
         tgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcccc                                                                                                                                          >  1:735993/1‑151 (MQ=255)
          gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccg                                                                                                                                         <  1:543447/151‑1 (MQ=255)
             ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                       <  1:768565/150‑1 (MQ=255)
              tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                                                                                     >  1:569027/1‑151 (MQ=255)
                tGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGtt                                                                                                                                    >  1:496674/1‑150 (MQ=255)
                   cATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                                                                                                 >  1:425982/1‑150 (MQ=255)
                    aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                                                                                                 >  2:580297/1‑149 (MQ=255)
                     tttATTCAACCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGTTTCAACAGCGAAGGTCACCCCCTTTTCGCCTTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTGGGTGTAATAACTTTTGTTTTCCCCCGCCGTtagtag                                                                                                                               <  1:120946/150‑1 (MQ=255)
                         ttCAACCCTTCAAGCGATTCAGATAGCCCCAGCTAAATCGGTTCAACAGCGAGGGTCCCCCCTTTTCCGCCGTTGTCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCg                                                                                                                             <  1:769248/148‑1 (MQ=255)
                            aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                                                                                        >  1:728622/1‑150 (MQ=255)
                            aaaCCTACAAGCGATTCAGAGAGCGCCAGTTTAATCGGTCCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGAAAACAAAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                                                                                       <  2:230553/151‑1 (MQ=255)
                             aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                                                                                        >  1:334410/1‑149 (MQ=255)
                                ctccAAGCGATCCAGATAGCGCCAGCTTAATCGGTTCAACGGCGAAGGTCACCCCTTTTTCGCCGTTGCCCGCGACAACATAACGCGGTGCCCTTTCTGTCTCGGTTTAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                    <  2:82210/147‑1 (MQ=255)
                                     aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                                                                                              >  2:251724/1‑151 (MQ=255)
                                        atttCAGATGGCACCAGTTAAATCTCTGCAACAGCGAAGGCCGCCCCCTTTCCGCTGTTGTCCGGGAAAATATAACGCAGTGCGCTTTATGTATCGGTGTAAAATCTTTTGTTTTACCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                   <  1:14856/141‑1 (MQ=255)
                                          ttCAGATAGCGCCAGTTAAAGGGGTTAAACAGGGAGGGTAACCCCCTTTTCCCCTTTGTCCGCAAAAAAATACCGCTTTCCCCCTTGTGTCCCGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                           <  2:255716/149‑1 (MQ=255)
                                             aGATAGCGCCATTTTAATCGTTAAAACAGCGAAGGCCCCCCCTTTTGCCCTTTTCTCCGCCAAAACATACCGCAGTCCTCCTTCTGTCGGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                      <  1:458492/151‑1 (MQ=255)
                                             aGATAGCGCCAGCTTAATCGTTCCAACAGCGAAGGTCACCCCCTTTCCGCCGTTGTCCGCGACAAAATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                                        <  1:263980/149‑1 (MQ=255)
                                               atggcgcCAGTTAAACCGGTTAAACAGCGAAGGCCACCCCTTTTTCGCCGTTTTCCGCAACAAAAAACCGCAGTCCCCCTTCTGTCTCGGTTTAATACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                    <  2:702709/148‑1 (MQ=255)
                                               aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                    >  1:680234/1‑151 (MQ=255)
                                                 aGCACCAGCTTAACCGGTACAACAGGGAAGCCCGCCCCCTTTTCGCGGTTGCCCGCAACAAAATAGCGCATTGCGCCTTCTGTCGGGGTGAAATCTCTTTTTTTTTCCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                      <  2:595572/147‑1 (MQ=255)
                                                   cgcCAGTTTAATCGGTTAAACAGCGAAGGCCACCCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCGGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                     <  2:590512/146‑1 (MQ=255)
                                                   cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                                >  2:322116/1‑151 (MQ=255)
                                                               cgttaaaaCAGCGAAGGACACCCCCTTTTCGCCGTGGCCCGCAACAACATAACGCAGTGCCCCTTTTGTCTGGGTTTAATAATGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                                                                                      <  1:342508/143‑1 (MQ=255)
                                                               cGGTACAACAGCGAGGCCCCCCCCTTTTGCGCCTTTGTCCGCGAAAACATAACGCAGTCCCCCTTCTGTCGCGGTTTAATAGTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                     <  1:732325/150‑1 (MQ=255)
                                                                ggTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                                   >  1:595059/1‑151 (MQ=255)
                                                                ggTCCAACAGGGAAGGCCACCCCCTTTCCGCCGTTGCCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                     <  2:343349/149‑1 (MQ=255)
                                                                     aaaaGCGAAGGTCACCCCCTTTTCTCTTTTTTCCGCGACAACAAAACGCAGTGCCCTTTCTTTCTGGGTTTATTACTGTTTTTTTTTCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                              <  1:319425/148‑1 (MQ=255)
                                                                     aaCAGCGAAGGTCCCCCCCTTTCCGCCTTTGCCCGCAACAACATACCGCAGTCCCCCTTCTGTCTCGGTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgg                                                                              <  2:153013/151‑2 (MQ=255)
                                                                      aCAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                             >  1:114976/1‑151 (MQ=255)
                                                                             aggccccccccTTTTCCGCGTTTGCCGGCGAAAACAAAGCGCGTTGCCCCTTTTGTCTGGTTAAAAAACGTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg                                                                         <  2:547065/142‑1 (MQ=255)
                                                                                 cacccccTTTCCTCCGTTGTCCCGAACAACAAACCGCGGTGCCCTTTTTGTTTCGTTTTAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCa                                                                    <  1:198232/146‑1 (MQ=255)
                                                                                    ccctTTTGCGCGTTTGCCCGCGAAAACATACCGCGGTGCCCCTTTTGTCTCGTTGTAATAGTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               <  1:113902/151‑1 (MQ=255)
                                                                                    ccccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               >  2:313794/1‑151 (MQ=255)
                                                                                    ccccTTTCCGCCGTTGTCCGCAACAACATAACGCGTTGCCCCTTCTTTCGGGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               <  2:736067/151‑1 (MQ=255)
                                                                                             cgccttggcccGCAAAAACATAACGCGGTCCCCTTTCTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTa                                                       <  1:314493/142‑1 (MQ=255)
                                                                                                  ttGCCCGCAAAAACAAACCGCATTCCACTTTTTGCCTGTGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                        <  2:194392/145‑1 (MQ=255)
                                                                                                   tGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 >  1:766593/1‑151 (MQ=255)
                                                                                                    gTCCGCAAAAACATAACGCGGTGCACCTTCTGTCTCGGTTTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 <  1:712936/150‑1 (MQ=255)
                                                                                                     tCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 >  1:762023/1‑149 (MQ=255)
                                                                                                      ccGCGACAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                                <  2:342440/149‑1 (MQ=255)
                                                                                                           acaacaTACCGCTGTCCACCTTCTGTCTCGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                          <  2:139038/150‑1 (MQ=255)
                                                                                                           acaacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                         >  1:455023/1‑151 (MQ=255)
                                                                                                               caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                                     >  1:761433/1‑151 (MQ=255)
                                                                                                                 tAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                                   >  2:348269/1‑151 (MQ=255)
                                                                                                                    cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCgg                                >  1:511605/1‑151 (MQ=255)
                                                                                                                     gagttgccctttCTGTCTGGTTGTATTTCGGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                                 <  1:498200/139‑1 (MQ=255)
                                                                                                                     gCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                               >  2:169552/1‑151 (MQ=255)
                                                                                                                         tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCtttt                           >  2:250125/1‑151 (MQ=255)
                                                                                                                           cACCTTCTGTCTCGTTGTAATACCTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                         <  1:63139/151‑1 (MQ=255)
                                                                                                                            cctttCTGTCTGTGTTTAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                          <  1:395699/146‑1 (MQ=255)
                                                                                                                              cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc                      >  2:479532/1‑151 (MQ=255)
                                                                                                                                tctgtctgGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                     <  2:558546/150‑1 (MQ=255)
                                                                                                                                     ctcGGTATAAAATTGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                 <  2:27968/149‑1 (MQ=255)
                                                                                                                                       cGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg             >  1:156766/1‑151 (MQ=255)
                                                                                                                                            taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg        >  1:122987/1‑151 (MQ=255)
                                                                                                                                             aatatCGGTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTGGCTCTGGCCCCTGGAt          <  1:489718/148‑1 (MQ=255)
                                                                                                                                             aataaCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt          <  2:728696/148‑1 (MQ=255)
                                                                                                                                              atatCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt          <  2:687906/147‑1 (MQ=255)
                                                                                                                                                aaCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa           <  2:115008/144‑1 (MQ=255)
                                                                                                                                                     ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGcc  <  2:458864/149‑1 (MQ=255)
                                                                                                                                                     ttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTa     <  1:545774/146‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCC  >  NC_000913/4111597‑4111892

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: