Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,492,8340TG48.4% 74.4 / 17.4 64G89G (GGT→GGG) tsgAputative transporter
Reads supporting (aligned to +/- strand):  ref base T (9/23);  new base G (30/0);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 5.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAC  >  NC_000913/3492689‑3492977
                                                                                                                                                 |                                                                                                                                                
cTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGAt                                                                                                                                             >  1:562789/1‑151 (MQ=255)
  gTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTggt                                                                                                                                             >  2:990494/1‑147 (MQ=255)
      ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGtt                                                                                                                                          >  2:170330/1‑148 (MQ=255)
        aGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGGTTTGTGTTCAGcc                                                                                                                                     >  2:545860/1‑151 (MQ=255)
           aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGGTGTTTAGc                                                                                                                                      >  2:441402/1‑147 (MQ=255)
              aGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAAcc                                                                                                                               >  1:386954/1‑151 (MQ=255)
                  aTACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTgg                                                                                                                            >  1:276978/1‑150 (MQ=255)
                     ccttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGgttgttgtttcgccccaaccctg                                                                                                                            >  2:968982/1‑133 (MQ=255)
                     ccttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGgcgc                                                                                                                         <  2:1039932/150‑1 (MQ=255)
                          accttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCACAGCCTGGCGCTGTtct                                                                                                                   >  1:473339/1‑151 (MQ=255)
                           ccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTtctc                                                                                                                  <  1:746776/151‑1 (MQ=255)
                               ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTATGTTTAGCCCCAGCCTGGGGCTGTTTTcgggg                                                                                                              >  1:503917/1‑149 (MQ=255)
                                  cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTGCGGTTTGATGTTTAGCCACAGCCTGGCGGGGTTTTcggcgg                                                                                                             >  2:104900/1‑149 (MQ=255)
                                  cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcggcgg                                                                                                             <  2:426885/149‑1 (MQ=255)
                                  cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGgcggcg                                                                                                           >  1:61513/1‑151 (MQ=255)
                                    aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGAt                                                                                                         <  1:824101/151‑1 (MQ=255)
                                     cGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGTTTTGTGTTTAGCCACAACCTGGCGCTGTTTTccgcgg                                                                                                             >  2:1020163/1‑146 (MQ=255)
                                       ccGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTTGGGCTGTTTTCGGCGGCGGTGtt                                                                                                      >  2:107906/1‑151 (MQ=255)
                                        cGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACCACCCGGCGCTGTTCTCGgcggcggggttt                                                                                                     >  2:748954/1‑145 (MQ=255)
                                          gCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑t                                                                                                  <  1:487924/151‑1 (MQ=255)
                                           cATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGGGGGCGGTGTTCA‑tt                                                                                                 >  1:430136/1‑151 (MQ=255)
                                               ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑Ttctc                                                                                              <  1:872160/150‑1 (MQ=255)
                                                 aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCgg                                                                                            <  1:270519/150‑1 (MQ=255)
                                                 aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTATGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGGGGGGTTTA‑TTTTCggg                                                                                           >  1:595973/1‑151 (MQ=255)
                                                 aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGACGGTTTGGTGTTTCGCCACAGCCTGGCGCTGTTCTCGGCGGGGATGTTCC‑TTCTCggg                                                                                           >  1:950269/1‑151 (MQ=255)
                                                  aTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTCAGGCACAACCTGGCGCTGGTCTCGGCGGGGGGGTTTT‑TTTTCgggg                                                                                          >  2:524411/1‑151 (MQ=255)
                                                    ctctATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGGCGGGTTGTTGGTTTGCCCCAGCCTGGCGCTGTTCTCCGCGGCGGTGTTTC‑TTTTCgggggg                                                                                        >  2:451285/1‑149 (MQ=255)
                                                      ctATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGggt                                                                                         <  1:92901/148‑1 (MQ=255)
                                                        aTCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTTATGTTCAGCCACAACCCGGCGGTGGTTTCTGCGGCCGTGGTTC‑TTCTCGgggggg                                                                                       >  2:708397/1‑148 (MQ=255)
                                                                cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATc                                                                            <  1:416783/151‑1 (MQ=255)
                                                                  aCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCAc                                                                          <  1:1040279/151‑1 (MQ=255)
                                                                  aCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTCTGCCACAGCCTGGCGCTGTTCTCTGGGGGGATGGTTC‑TTTTCGGGGGGGGCAGCGGccacac                                                                          >  2:351148/1‑146 (MQ=255)
                                                                           tGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGGTGCCGGCTTGGTGTTAAGCCCCCGCCCGGCGCTCTTTTCGGTGGCGCTGGTTC‑TTCTCGGGGGGGGCCGCGGGATCGCCATggctg                                                                  >  1:365124/1‑146 (MQ=255)
                                                                            gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGAtt                                                                <  2:130547/151‑1 (MQ=255)
                                                                            gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGGTGCCGGGTTGGTTTTTTGCCACAACCAGGCGCTGTTTTCGGGGGCGAGGTTCACTTTTCGGGG‑GGGCCGCGGCATCACCATg                                                                      >  2:743782/1‑145 (MQ=255)
                                                                             aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTTATGTTTAGCCACAGCCTGGCGCTGTTCTCCGCGGCGGTGTTTA‑TTTTCGGGGGGGGCAGCGGCACCACCATTTTGAtt                                                                >  1:236276/1‑150 (MQ=255)
                                                                             aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTg                                                               <  2:468053/151‑1 (MQ=255)
                                                                               ggAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTgg                                                              <  2:944741/150‑1 (MQ=255)
                                                                                 aaaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAc                                                           <  2:950775/151‑1 (MQ=255)
                                                                                 aaaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGGGATGTTTA‑TTCTCGGGGGGGGCCGCGGCATCACCATGTCGGTTgg                                                              >  2:156474/1‑148 (MQ=255)
                                                                                  aaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTCCGGTTTGTTGTTTAGCCCCAACCTGGCGCTGTTCTTGGCGGGGGTGTTCA‑TTCTCTGGGGGGGCCGCGGCATCACCACGTCggt                                                                 >  1:464064/1‑142 (MQ=255)
                                                                                      gTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGGTGTTCAGCCACAGCCTGGCGCCGTTCTCGGCGGCGGTGTTCC‑TTTTCGGGGGGGGCAGCGGCATCACCCTGTCGATTGGGACATTcc                                                      >  1:497179/1‑151 (MQ=255)
                                                                                       tCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCCGGCGCTGTTCTCGGCGGCGGTGTTCA‑TTCTCGGGGGGGTCAGCGGCATCACCCTGTCGATTTGGACAtt                                                        >  1:49304/1‑148 (MQ=255)
                                                                                        cccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGTCGGGTTGGTGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCT‑TTTTCGGGGGGGGCAGCGGGCTCACCCCGTTGGTTGGTACATTcc                                                      >  1:61587/1‑149 (MQ=255)
                                                                                               aaaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTTTTCTCGGCGGCGGTGTTTA‑TTCTCGGGGGGGGCAGCGGCCTCCCCATGTCGATTTGTACAtt                                                        >  1:96350/1‑140 (MQ=255)
                                                                                                  aCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACaaa                                          <  1:33147/151‑1 (MQ=255)
                                                                                                    gCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGATGTTTAGCCACCGCCCGGCGCTGTTCTCGGGGGGGGTGTTTT‑TTTTCGGGGGGGGCAGCGGCATCACCATGTCGGTTGGGTCATTCCTGGTAACACaaa                                          >  2:716288/1‑149 (MQ=255)
                                                                                                         tACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATg                                    <  1:189500/150‑1 (MQ=255)
                                                                                                         tACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATg                                    <  2:457842/150‑1 (MQ=255)
                                                                                                         tACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGGTTGTTGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCA‑TTTTCGGGGGGGGCAGCGGCATCACCCTGTCGGTTGGTACATTTCTTGTAACACAAAAGTAt                                     >  1:464079/1‑149 (MQ=255)
                                                                                                          aCGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTTCGGGTTGTTGTTTTGCCACAACCTGGGGCTGCTCTCTGCGGGGAGGTTGA‑TTTTCGGGGGGGGGCACGCGGTCACCCTGACGCTTGGTTATATCCTGGTTACACAAAAGTATg                                    >  2:588275/1‑149 (MQ=255)
                                                                                                          aCGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGaa                                  <  1:465725/151‑1 (MQ=255)
                                                                                                                  gCTTTCTCCTGATGGTGCTGGCGGTTGTCGGGTTTATTTTCAGCCACAGCCTGGCGCTGTTCTCGGGGGGGGTGTTTA‑TTCTCGGGGGGGGCAGGGGCCTCACCATGTCGATTGGGATATTCCCTGTAACACAAATGTATGaa                                  >  2:139789/1‑143 (MQ=255)
                                                                                                                   cTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTTATGTTTAGCCACAGCCTGGCGCTGTTTTCGGCGGCGGTGGTTA‑TTCTCGGGGGGGGCAGCGGCCTCACCCTGTCGGTTGGTACATTTCTGGTAACACAAAAGTATGAAggggg                             >  1:496281/1‑145 (MQ=255)
                                                                                                                   cTTTCTCCTGATGGTGCTGGCGGGTGCCCGATTGGGGTTCAGCCACACCCTGGCGGTGTTCTCCGCGGCGGCGTTTA‑TTATCCGGGGGGTAAACGGGATCACCCAGTCGCTTGGTTTATTCCTTGTAACAACAAAGTATGCAgg                                >  1:346416/1‑144 (MQ=255)
                                                                                                                     ttCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAg                         <  2:532537/149‑1 (MQ=255)
                                                                                                                     ttCTCCTGATGGTGCTGGCGGTTGCCGGGTTTGTGTTTAGCCACAGCCTTGCGCTGTTCTCGGCGGCGGTGTTCA‑TTCTCGGGGGGGTCAGCGGCATCACCATGTCGATTGGTACATTTCTGGTAAAACAAAAGTATGAAGGGCGTCa                          >  1:396594/1‑148 (MQ=255)
                                                                                                                      tctcCTGATGGTGCTGGCGGTTGCCGGTTTGATGTGCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTCGTCAGCGACATCGCCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAgcg                       <  2:657623/150‑1 (MQ=255)
                                                                                                                       ctcCTGATGGTGCTGGCGGTTTTCGGGTTTTTGTTTAGCCACAGCCTGGCGCCGTTCTCGGGGGGGATGTTCA‑TTTTCTGGGGGGGCAGCGGCCTTACCCTGTCGGTTGGGATATTTCTGGTAAAACAAATGTATGAAGGGGGGCAgcg                       >  2:105543/1‑149 (MQ=255)
                                                                                                                           tGATGGTGCTGGCGGTTGCCGGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTTA‑TTTTCGGGGGGGGCAGCGGGATCACCATGTCGGTTGGGACATTTCTGGGAAAACAAATGTATGAAGGGGGGCAGCGCgggt                  >  1:243760/1‑148 (MQ=255)
                                                                                                                            gATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTc                 <  1:555253/150‑1 (MQ=255)
                                                                                                                                     ggCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGCGTTGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGtt        >  2:571193/1‑150 (MQ=255)
                                                                                                                                      gCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAt      >  2:893141/1‑151 (MQ=255)
                                                                                                                                         gTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTa   <  2:276963/151‑1 (MQ=255)
                                                                                                                                           tGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc  <  1:98246/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCA‑TTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAC  >  NC_000913/3492689‑3492977

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