Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,6940CT45.2% 72.0 / 29.6 62F191F (TTC→TTTyafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (7/27);  new base T (28/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGA  >  NC_000913/231554‑231843
                                                                                                                                            |                                                                                                                                                        
ccTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGa                                                                                                                                               <  2:186512/151‑1 (MQ=255)
      acCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACCTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTCAATGCCTgggg                                                                                                                                              >  2:801502/1‑144 (MQ=255)
         cGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCtcggggcggag                                                                                                                                         >  1:346934/1‑139 (MQ=255)
             tGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAAAGCCTGGGGCCGGAgga                                                                                                                                      >  1:979763/1‑145 (MQ=255)
              gTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTCAATGCCCGGAGCCGGAGAAgg                                                                                                                                   >  2:935077/1‑148 (MQ=255)
              gTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAATGCCTTGGGCCGTAGaaag                                                                                                                                    >  1:801111/1‑146 (MQ=255)
                   tggtggTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAATAGCTGGGGCCGTTGAAAGagt                                                                                                                                >  1:103940/1‑144 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTTATGCCCGGAGGCGTAGaa                                                                                                                                      >  2:275985/1‑140 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgcg                                                                                                                          <  1:579158/151‑1 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGCCGATTTTCATGACCGGGGCCGCAGAAGGacggaag                                                                                                                            >  2:81324/1‑143 (MQ=255)
                       ggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGt                                                                                                                         <  2:1130137/149‑1 (MQ=255)
                          tAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTCTATCCCCGGGGCCGGAGAAGGATGaa                                                                                                                              >  1:1041537/1‑141 (MQ=255)
                           aaTATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTTATGCCTGGGGGCGCAGAAGGAAGAAAGCGtt                                                                                                                        >  2:782708/1‑146 (MQ=255)
                             tataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATTCCTTGAGGCCGAGAAGGGATGAAGGGGtttttg                                                                                                                    >  1:628706/1‑145 (MQ=255)
                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATGCCTTGGGCCGTAGAAAGATGAAAGCGTTATATCGCt                                                                                                                >  2:807060/1‑151 (MQ=255)
                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGGCGTAGGAGGAAGAAAGCGGTATATTGCt                                                                                                                >  2:732615/1‑151 (MQ=255)
                               taCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGCCGGAGAAGGATGAAAGCGGTTTATTGCtt                                                                                                               >  2:114712/1‑151 (MQ=255)
                                 cacGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTg                                                                                                             <  1:541375/151‑1 (MQ=255)
                                  acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGGGCCGGAGGAGGATGAACGCGGTATATTGCttt                                                                                                              >  2:1013079/1‑149 (MQ=255)
                                                     gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGGCGCAGAAGGGGGAACGCGGTTTATTGCTTTTCGCGGGGAAAGTTGCTg                                                                                            >  2:21996/1‑148 (MQ=255)
                                                     gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTgcgc                                                                                        <  2:500805/151‑1 (MQ=255)
                                                      ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCCGGGGCCGCAGAAGGGTGAAAACGGTTTTTCGCCTTTCGCCGGGAAAGTCg                                                                                               >  1:931690/1‑144 (MQ=255)
                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCa                                                                                      <  2:356626/151‑1 (MQ=255)
                                                         gTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg                                                                                    <  2:861435/151‑1 (MQ=255)
                                                           ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTg                                                                                  <  2:773699/151‑1 (MQ=255)
                                                            gATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAATTCCCTGAGCCCGAAAAGGAAGAACGCGGTATATAGCTTTTCGCCGGGAAAGTCGCCGCGCCCgggg                                                                                  >  1:120244/1‑148 (MQ=255)
                                                                 cTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTTAAAGCCTGGGGCCGGAGAAGGATGGACGCGGTATATCGCTTTTCGCGGGGAATGTCGCCGCGCCCGGGGCGgttt                                                                            >  1:694900/1‑151 (MQ=255)
                                                                  tataGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTTATTCCTGGAGCCGGAGAAGGATGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTTGCCGCGCCCgggggggtttt                                                                           >  1:58967/1‑142 (MQ=255)
                                                                          gCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCgatgat                                                                   <  1:873822/151‑1 (MQ=255)
                                                                           cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTCAATGCCTGGAGCCGTAGAAAGAAGAACGAGGTATATAGCTTTTCGCGGGGAAAGTCGCCGCGCCAgggggggttttcc                                                                         >  2:671828/1‑133 (MQ=255)
                                                                             gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGCCGTAGAAGGAAGAACGCGTTATATCGCCTTTCGCGGGGAATGTCGCCGCGCCAGGGGGGTTTTTCCgatga                                                                    >  2:521103/1‑147 (MQ=255)
                                                                              ttACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGGAGAAGGGTGAAAGCGTTATTTCGCTTTTCGCGGGGAATGGCGCTGCGCCCGGGGCGGTTTTCCGATGATCAg                                                                >  2:212148/1‑150 (MQ=255)
                                                                                aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAgcg                                                              <  2:891679/150‑1 (MQ=255)
                                                                                    ccTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGCCGTAGAAGGAGGAACGCGGTATATCGCTTTTCGCGGGGAAAGGCGCCGCGCCCGGGGGGGTTTTCCCATGATCAGCGcc                                                            >  1:1318/1‑148 (MQ=255)
                                                                                      tATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGATTTAATGCCCGGGGCCGTAGAAGGATGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTTGCTGCGCCAGGGGGGGTTTTCCGGTGATCAGCGCCGTCgg                                                       >  1:612458/1‑151 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAAAGCCTGGAGCCGTAGAAGGATGAAAGCGGTATATTGCTTTGCGGGGGGAAAGTTGGTGCGCCAGGGGGGGTTTTACGGTGATTAGGGCCGGCggg                                                      >  2:74346/1‑151 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCgg                                                       <  1:902182/150‑1 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCgg                                                       <  2:650974/150‑1 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCg                                                        <  2:965990/149‑1 (MQ=255)
                                                                                         tGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGc                                                     <  1:488173/150‑1 (MQ=255)
                                                                                          ggCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCg                                                    <  2:1120549/150‑1 (MQ=255)
                                                                                            cGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATACCTCGGGGCCGCGGAAGGGGAAAGCGGTTTATTCCTTTTCGCGGGGAAAGTCGCCGCGCCCGGGGCGGTTTTCCgatg                                                                     >  1:598794/1‑131 (MQ=255)
                                                                                             gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTg                                                <  1:316153/151‑1 (MQ=255)
                                                                                             gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTg                                                <  1:801474/151‑1 (MQ=255)
                                                                                                  gATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTCCATGCCTGGAGCCGTAGAAGGATGAACGCGGTATATCGCTTTGCGCGGGGAATGTCGCTGCGCCCGGGGCGTTTTTCCGCTGATCAGCGCCGTCGGGCGTTTGGTCg                                            >  2:1011915/1‑150 (MQ=255)
                                                                                                      gCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGGCGTAGAAGGATGAACGCGGTATATCGCTTTGCGCGGGGAATGTCGCTGCGCCAGGGGGGTTTTATCGATGATCAGCGCCGGCGGGCGGTTGGTCGccc                                         >  2:105852/1‑149 (MQ=255)
                                                                                                       cTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTTAAAGCCCGGGGCCGGCGAAGGGAGGAAGCGTTATATTGCTTTGCGCGGGAAAAGGCGCTGAGCCCAGGGGGGTTTTACGGAGATAAGCGCCGGCCGGGGGTTGGTCggc                                          >  2:545562/1‑145 (MQ=255)
                                                                                                         caccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAAATCCTCGAGGCCGAGAAGGATGAACGCGTTTTATCGCTTTTCGCGGGGAAAGTCGCTGCGGCAGGGGCGGTTTTCCCGTTATCAGCGGCGGCGGGCGTTTTGTCGCCCGCTc                                     >  1:24009/1‑150 (MQ=255)
                                                                                                          accacAGCGGCCCGGTCATTATGGCGGGGGGTTTTTATTCCCCGAGGCGGGGGAGGGAGAAAGCGGTATATTGCTTTGTGCGGGGAAAGTCGCTGCGCCCGGGGGGGTTTTCCGATGGTCTGCGCCGTCGGGGGGTTGGTCGccc                                         >  1:46589/1‑145 (MQ=255)
                                                                                                          accacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCTCGGGAAATGTCGCTGCACCGGGTGCGTTTTACCGATGATCAGCCCCGTCGGGCGTTTGGTCTCCCGCTCg                                    <  1:937449/150‑1 (MQ=255)
                                                                                                            cacaGCGGCCCGGTCATTATGGCGGGAGATTTTAAAGCCCCGGGGCGGAGAAGGGAGAAAGCGGTTTATCGCTTTTCGCGGGGAAAGTTGCTGCGCCCGGGGCGGTTTTCCGAAGATCAGCGCCCCCGGGGGGTTGGGCGgccccccggt                                  >  1:50928/1‑143 (MQ=255)
                                                                                                               aGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTg                              <  2:897043/151‑1 (MQ=255)
                                                                                                                  ggCCCGGTCATTATGGCGGGAGATTTTAATTCCCGGAGCCGGAGAAGGATGAACGCGTTATATTGCTTTTCGCGGGGAATGTCGCTGCGCCAGGGGCGTTTTTTCGATGATCAGCGCCGGCGGGCGGTTTGGCGCCCCCCCGAtttttttt                           >  1:959235/1‑151 (MQ=255)
                                                                                                                    cccGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGtttt                          <  1:419200/150‑1 (MQ=255)
                                                                                                                       ggTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTAc                      <  2:467552/151‑1 (MQ=255)
                                                                                                                          cATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCg                    <  1:111094/150‑1 (MQ=255)
                                                                                                                          cATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGt                   <  2:431952/151‑1 (MQ=255)
                                                                                                                           aTTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTTGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTg                  <  2:1099040/151‑1 (MQ=255)
                                                                                                                              aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTc               <  2:971176/151‑1 (MQ=255)
                                                                                                                                ggCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCt              <  1:1091390/150‑1 (MQ=255)
                                                                                                                                 gCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTg             >  2:561209/1‑150 (MQ=255)
                                                                                                                                     gagaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGt        <  1:233408/151‑1 (MQ=255)
                                                                                                                                            cAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGa  >  2:193879/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                        
CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGA  >  NC_000913/231554‑231843

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: