Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4680AG46.9% 88.7 / 19.2 64D55G (GAT→GGT) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (3/31);  new base G (30/0);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 3.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079320‑1079616
                                                                                                                                                    |                                                                                                                                                              
ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                                                                                              <  1:587302/151‑1 (MQ=255)
ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                                                                                              <  2:385681/151‑1 (MQ=255)
     gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAg                                                                                                                                                          <  2:753488/150‑1 (MQ=255)
     gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAg                                                                                                                                                          <  1:1199074/150‑1 (MQ=255)
     gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                                         <  1:356501/151‑1 (MQ=255)
     gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTTTGAgg                                                                                                                                                         >  1:365617/1‑150 (MQ=255)
           gACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggctgg                                                                                                                                                  <  1:935107/151‑1 (MQ=255)
              aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCgggcgg                                                                                                                                                   >  2:543378/1‑144 (MQ=255)
                 ttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGtt                                                                                                                                          <  1:56509/150‑1 (MQ=255)
                   gtgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAAt                                                                                                                                      <  2:333555/151‑1 (MQ=255)
                     gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATg                                                                                                                                     <  2:888190/150‑1 (MQ=255)
                         atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATggg                                                                                                                                   <  1:959960/148‑1 (MQ=255)
                                tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTTGGCCATTCGTGACGGCATTATCGGCGGGGGGGGGGGGGATGAGGGGCTGGGTGGTTATAGGGGTGGCgg                                                                                                                          >  2:1110925/1‑149 (MQ=255)
                                tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCggg                                                                                                                         <  2:524424/150‑1 (MQ=255)
                                 ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGgcg                                                                                                                      <  1:11246/151‑1 (MQ=255)
                                   cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGGGGGGCCGGGTGTTTATtgggggggcggggg                                                                                                                        >  1:889881/1‑139 (MQ=255)
                                   cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGTGTAGGCGGGCGGGTGTTTATGGGGGTGCCggggg                                                                                                                       >  2:912874/1‑147 (MQ=255)
                                    aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGTCGGGTATGAGCGGCTGGCTGGTAATTGGGGTGGCCgg                                                                                                                          >  1:312389/1‑146 (MQ=255)
                                      gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAAt                                                                                                                                      <  1:1196340/132‑1 (MQ=255)
                                       atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGGATTGGGGGGTGGGTGTTTATTGGGGTGGCGGGGGCTg                                                                                                                  >  1:886953/1‑148 (MQ=255)
                                         attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttttt                                                                                                            <  2:20378/151‑1 (MQ=255)
                                             attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttcttt                                                                                                        <  2:444686/151‑1 (MQ=255)
                                             attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTTTGAGGGGGCGGGTGTTTATTGGGGTGCCGGGGGCCGGtttttttt                                                                                                        >  1:849605/1‑151 (MQ=255)
                                              ttgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTAGAGGGGGGGGGGTTATTGGGGTGTCGGGGGGTGTtttttttt                                                                                                        >  2:965360/1‑150 (MQ=255)
                                                gggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTcc                                                                                                      <  1:280033/150‑1 (MQ=255)
                                                        cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGGGGGTGGGTGTTTATGGGGGTGGCGGGGGCGGTTTTTTTTTTCGGGGttt                                                                                               >  1:894602/1‑149 (MQ=255)
                                                         gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCg                                                                                            <  2:418518/151‑1 (MQ=255)
                                                         gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGGGGGGGGGGGGGTTTATGGGGGTGGCGGGGGGGGTTTTTTTTTTCgggggttt                                                                                              >  2:389858/1‑144 (MQ=255)
                                                          ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGCGGGGCCTTATCGGGTGGGGGTTTTTTGGGTTGCCGGGGCGCGGTTTTTTTTCCgggggttt                                                                                              >  1:810985/1‑143 (MQ=255)
                                                               cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGCCGGGTATTAGGGGCTGGGTGTGTATTGGGGGGGCGGGGGGCGTTTTTTTTTTCGGGATTTCCGaaa                                                                                         >  2:381187/1‑148 (MQ=255)
                                                                gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTg                                                                                     <  2:789946/151‑1 (MQ=255)
                                                                     aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGTGGGGTATTAGGGGCCGGCTGTTTATTGGGTTGCTGGGGGCTGTTTTTTTTTTCGGGGtttttcgaaaacgg                                                                                    >  1:827753/1‑136 (MQ=255)
                                                                         aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGGGGGGCTGGGGGGTTATGGGGGTGGCCGGCGGTGTTTTTTTTTTCCGGGTTTTCCGAAACTGGATCGccc                                                                             >  2:449421/1‑149 (MQ=255)
                                                                          aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGCTGGCTGTTTATTGGGGTGGCGGGGGCTGGTTTTTTTTTCCGGGTTTTCCAAAGCTGGGTCGCCCtt                                                                           >  2:235962/1‑151 (MQ=255)
                                                                           aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                         <  2:608173/151‑1 (MQ=255)
                                                                           aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGCGGGTGGGCGGTTATGGGGGTGGCGGGGGCCGTTTTTTTTTTTCGGGTTTTCGGAAGCTGGGTCGCCCtt                                                                           >  2:608300/1‑150 (MQ=255)
                                                                           aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGGGCGGGTGGGGGGTTATTGGGGTGTGGGGGGGCGGTTTTTTTTTCTGGGGTTTCGGGAACTGGGTCGcccct                                                                           >  1:381954/1‑147 (MQ=255)
                                                                            aCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGGGGGTGTTAATTGGGGTTGCGGGGGCCGGTTTTTTTTTCCGGGTTTCCCGAAACGGGATCggcctt                                                                           >  1:255792/1‑144 (MQ=255)
                                                                                tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGTGGGGGTTTATGGGGGTGGCGGGGGGGGTTTTTTTTTTCGGGGTTTTCGGAAGCTGGATCGCCCTTGGcc                                                                      >  2:805449/1‑149 (MQ=255)
                                                                                     acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACAtt                                                               <  2:240566/151‑1 (MQ=255)
                                                                                     acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGGGGCTGGGTGTTTATGGGGGTGGCGGGGGCCGTTTTTTTTTTCGGGGTTTCCGCAAGCTGGATCGCCCTTGGCCTGACAtt                                                               >  2:444700/1‑151 (MQ=255)
                                                                                          atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTAGGGGGTGGCTGTTAATTGGGGTGGCGGGGGGTGGTTTTTTTTCCGGGGTTTTCGACAGGTGGGTGGCCCTTGGGCTGAAATTAgggg                                                          >  2:443143/1‑149 (MQ=255)
                                                                                          atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGCTGGGTATGGGCGGCTGGGTGGTTAAAGGGGTGGCGGGCGGTGTTTTTTTTTTTGGGGTTTTCCGAAGGTCGCGTGTCATTAGCCTGACATTac                                                             >  2:844025/1‑147 (MQ=255)
                                                                                            tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                         <  1:1204740/150‑1 (MQ=255)
                                                                                                 gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGGGGGTGTTTATGGGGGTGGCGGGGGGTGTTTTTTTTTCCGGGGGTTTCGAAAGCTGGGTCGCCCTTGGCCTGACATTAGGCGCGTggg                                                    >  2:106355/1‑149 (MQ=255)
                                                                                                   ggTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGGGGGTTGAGCGGGTGGCGGGTATTGGGGGGGGCGGGCGGCGGTTTTTTTTTCGGGGGTTTCCGAAACCGGGTCGCCCTTGGGCCTGCATTTGgcgggggg                                                     >  1:659791/1‑142 (MQ=255)
                                                                                                      cGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGTGGGTATGAGCGGGTGGGTGTTTAAGGGGGTTGCGGGCGCTGGTTTTTTTTTCGGGGGTTTCGGAAGACGGATCGCCCTTTGCCCGTCCTTAGGCGCGGGGGTTaa                                                >  2:931223/1‑149 (MQ=255)
                                                                                                       gTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTg                                             <  2:1111406/151‑1 (MQ=255)
                                                                                                       gTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATTAGCGGCTGGGTGGTTATTGGGGTTCCGGGGGCTGGTTTTTTTTTTCGGGTTTTCGAAAGCTGGAGCGCCCTTGGCCTGACATTAGGGGCGTGGGTTAACt                                              >  2:226135/1‑150 (MQ=255)
                                                                                                         aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGa                                           <  1:1002426/151‑1 (MQ=255)
                                                                                                            cTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGCGGCTGTTTATGGGGGTGCCGGGCGCCGGTTTTTTTTTCGGGGTTTTCGAAAAGTGGATCGCCATTGGCCCGACATTAAGCGCGGGGGTTAACTgg                                            >  2:373730/1‑147 (MQ=255)
                                                                                                                 gCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtg                                   <  1:1152149/151‑1 (MQ=255)
                                                                                                                    aTTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGTGGGTGTTTATGGGGGTGGCGGGGGCTGGTTTTTTTTTCGGGGGTTTCGAAAACTGGGACGCCCCTTGGCTGACATTAGGCGCGGGGGTTaaattggaaccggggg                                  >  2:1094783/1‑135 (MQ=255)
                                                                                                                      tCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCg                               <  2:213104/150‑1 (MQ=255)
                                                                                                                        gTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGc                            <  2:304573/151‑1 (MQ=255)
                                                                                                                           aCGGCATTATCGGCGGGTGCGGCGGGTATGAGGGGGTGGCTGGTTAATGGGTTTGCGGGGGCCGTTTTTTTTTTCCGGGGTTTTCAAAGCTGGATCGCCCTTTGGCCGAACTTTTGCGCGGGGGTTTAATTGAAGCGGGGGGgcgggcgg                          >  1:1157224/1‑150 (MQ=255)
                                                                                                                            cGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGtt                        <  2:628229/151‑1 (MQ=255)
                                                                                                                                  tATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtg                   <  1:885529/150‑1 (MQ=255)
                                                                                                                                   aTCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtgt                  <  1:93584/150‑1 (MQ=255)
                                                                                                                                     cggcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGc                <  2:499951/150‑1 (MQ=255)
                                                                                                                                        cggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAc            >  2:13957/1‑151 (MQ=255)
                                                                                                                                                 cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaa   <  1:959442/151‑1 (MQ=255)
                                                                                                                                                 cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATaca    <  2:774209/150‑1 (MQ=255)
                                                                                                                                                 cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTAGTGGCCGGGCGGTTGCGTGTGCATACCGAATaca    >  1:1127261/1‑150 (MQ=255)
                                                                                                                                                  ggATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  <  2:144961/151‑1 (MQ=255)
                                                                                                                                                   gATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGCATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  >  2:288375/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                              
CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079320‑1079616

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