Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3490CA41.5% 120.7 / 36.7 65V69F (GTT→TTT) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (6/32);  new base A (27/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 9.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAT  >  NC_000913/266208‑266494
                                                                                                                                                |                                                                                                                                                   
tCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAgg                                                                                                                                             <  2:1314079/150‑1 (MQ=255)
      cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑Ac                                                                                                                                      <  1:799178/150‑1 (MQ=255)
         ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGa                                                                                                                                  <  1:566832/151‑1 (MQ=255)
         ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGGAAAAACAGCAAGGGGGG‑Gcc                                                                                                                                     >  2:821058/1‑148 (MQ=255)
          gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGAt                                                                                                                                 <  2:747824/151‑1 (MQ=255)
           tATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGGAAAAACCGCAAggggggggcgg                                                                                                                                    >  2:305706/1‑143 (MQ=255)
            aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATa                                                                                                                                <  1:1312052/150‑1 (MQ=255)
            aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑AGGAAAAACCGCAAAGGGGG‑GCCCGaga                                                                                                                                >  1:561292/1‑148 (MQ=255)
             ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑AGGGAAAAACCTCAAAGGGGGG‑CCCgggga                                                                                                                                >  1:378043/1‑147 (MQ=255)
               ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑TGTAAAAACTCCAAAGGGGG‑GGCGGAGaa                                                                                                                               >  1:1297816/1‑148 (MQ=255)
                cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑AGTAAAAACTGCAAGGGGGG‑ACCGGATa                                                                                                                                >  1:763100/1‑146 (MQ=255)
                     tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCaa                                                                                                                       <  2:76439/150‑1 (MQ=255)
                      gATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAAcc                                                                                                                     <  1:453990/151‑1 (MQ=255)
                           tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCtgtg                                                                                                                 <  1:1169239/150‑1 (MQ=255)
                            aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAAACTGCAAGGGGGG‑GCCGGATAAGAAGGCAACCTGTGAc                                                                                                               >  2:220240/1‑151 (MQ=255)
                             tCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAAACAGCCAAGGGGG‑GCCGGGTAAAAAGTCAAACCGTGGcc                                                                                                              >  1:1057395/1‑151 (MQ=255)
                              cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGGAAAAACCGCAAAGGGGG‑GCCGGGGAAGTAGGCAACCTGTGAccc                                                                                                             >  2:817747/1‑150 (MQ=255)
                              cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGAAAAAACCACAAGGGGGG‑GCCGAGTAAGTGGGGAAACTGTGagc                                                                                                              >  2:480287/1‑148 (MQ=255)
                                   aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCAttc                                                                                                          <  2:1224372/149‑1 (MQ=255)
                                              aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAACCGCAAGGGGGG‑GCCGGATAAGTAGTCAACCTGTGACCATTTTTC‑GCTGAcatc                                                                                             >  2:596051/1‑150 (MQ=255)
                                              aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAAACTGCAAGGGGGG‑GCCGGGGAAGAAGTGAAACTGTGAACAttcttttctggca                                                                                                >  2:45942/1‑140 (MQ=255)
                                                tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAA‑CGGAAAACTTCCTCGGGGGG‑GCCGGGAAAG‑AGTCAAACCCTGTGCACACTTC‑TCTGAcatcctt                                                                                          >  1:1072962/1‑148 (MQ=255)
                                                 catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTTAAAACTGCAAAGGGGG‑GCCGGGGAAATAGTCAAACCGTGACCaattttttct                                                                                                    >  1:267918/1‑132 (MQ=255)
                                                  atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CAAAAAAAACTCAAGGGGGG‑GCCGGGTAAATAGTCAAACTGTGACCCTTCTTC‑GCTGAcataaaccc                                                                                        >  2:743448/1‑145 (MQ=255)
                                                       tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCAGCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCtggt                                                                                     <  1:637968/149‑1 (MQ=255)
                                                               cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACcac                                                                           <  2:71097/151‑1 (MQ=255)
                                                                 gTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGGAAAAACCGCAAGGGGGG‑GCCGGGGAAAAAGTCAACCCGTGACCCTTCTTC‑TCTGACATCATCCTGGGGGGTGACcac                                                                           >  2:931008/1‑149 (MQ=255)
                                                                   tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CCGAAAAACCGCCAGGGGGG‑GCCGGGTAAAAAGGTAAACCGTGACCATTCTTT‑CCCGACATCATCCTTGGGGCTGGcccc                                                                           >  2:204996/1‑145 (MQ=255)
                                                                     tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAAATGCAAGGGGGG‑ACCGGATAAATAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGGGGGTGAcc                                                                             >  2:410753/1‑143 (MQ=21)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGT‑A‑C‑TAAAACCCGCAAGGGGGGAACCGGA‑AAAAAATCAAACTGTTAACCTTTTTT‑TCTTTCATCATCCtgctggtctgccccaaccctt                                                                    >  2:74457/1‑136 (MQ=255)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGTCTGACCACATCACt                                                                     <  1:1289039/150‑1 (MQ=255)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACt                                                                     <  2:258286/150‑1 (MQ=255)
                                                                       tCCTACATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTg                                                                    <  1:76264/150‑1 (MQ=255)
                                                                        ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                  <  2:1258681/151‑1 (MQ=255)
                                                                           gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCAGTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGaatc                                                                <  1:177423/150‑3 (MQ=255)
                                                                            cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                              <  1:495806/151‑1 (MQ=255)
                                                                             atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑AGGAAAAACTGCAAAGAGGG‑GCCGGATAAATAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGGGGGTGACCCCATCCCTCAAAcac                                                              >  1:1201521/1‑148 (MQ=18)
                                                                              tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAACTGCAAGGGGGG‑GCCGGGAAAGTAGTCAACCTGTGGCCCTTCTTT‑GCCTACCTCATTCTTGTGGGTGACCACATCCCTGaaaaaggg                                                            >  2:131218/1‑146 (MQ=255)
                                                                                 ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgg                                                         <  2:1096822/151‑1 (MQ=255)
                                                                                 ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAAACCGCAAAGGGGG‑GGCGGGGAAAAAGTCAACCCGTGGACATTCTTT‑GCTGACATCATCCCTGGGGGTGACCAAAACACTGAAAAGCGgggg                                                         >  2:1297641/1‑151 (MQ=255)
                                                                                   aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑AGTAAAAACAGCAAAGGGGG‑AACGGGTAAGAAAGCAACCTTTGAACATTTTTCTGCTGCAATCAT‑CTTGTGGCTGACCACCTTACTAAAAAACGgggggg                                                       >  1:898842/1‑151 (MQ=255)
                                                                                     gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAACCCCACGGGGGG‑GGCGGATAAGAAGTCAACCCGTGAACCTTCTTC‑GCTGTCATCATCCTTGTGGCTGGCCACACCACTGAAAAGCGgggggggg                                                     >  1:971876/1‑149 (MQ=255)
                                                                                     gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑C‑GGTAAAAACTGAAAGTAGGG‑ACCCGGAAAGCGATAAACCCCTCACCATTTTTC‑TCTTCGACAATCCCGGGGGCTGGCCCCATCCCTGAAAACCGGTggggg                                                      >  2:1110177/1‑149 (MQ=255)
                                                                                      aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                    <  2:125851/151‑1 (MQ=255)
                                                                                       ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAg                                                   <  2:852972/151‑1 (MQ=255)
                                                                                       ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAACTGCAAGGGGGG‑ACCGGAAAAAAAGTCAAACTGTGACCATTCTTC‑CCTGACAACATCCCGGTGGCTGACCACATCACTGAAACGCGgggggggg                                                     >  2:907149/1‑147 (MQ=255)
                                                                                          agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTc                                                <  2:139417/151‑1 (MQ=255)
                                                                                            tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGGAAAAACTGCAAGGGGGG‑GCCGGGGAAGTAGGCAAACCGTGACCCTTCTTC‑GCCGACATCATCCCGGGGGGTGGGCCCATCACTGAAACGCGggggggggggggcc                                               >  1:1121858/1‑142 (MQ=255)
                                                                                              gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAAACTGCAAAGGGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGGGGCTGACCACATCACTTAAAAGCGGGGGGGTGGGGGCCtt                                             >  2:514947/1‑150 (MQ=18)
                                                                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                           <  2:635894/150‑1 (MQ=255)
                                                                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                           <  1:200497/150‑1 (MQ=255)
                                                                                                     ccATCCTCGCAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAACCTGCAAGGGGGG‑GCCGGGTAAGAAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGCGGGGGGAGGGCCTTTAAttt                                       >  2:26040/1‑149 (MQ=18)
                                                                                                       aTCCTCGAAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTg                                   <  2:1087161/151‑1 (MQ=255)
                                                                                                          cTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTCAAACCTGCAAGGTGGG‑GCCGGATAAGTAGTCAACCCGTGACCATTCTTT‑GCTGACATCATCCTGGTGGGTGACCACAACACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACt                                    >  2:504208/1‑147 (MQ=255)
                                                                                                              cAGCCATCGGGCGACTGGC‑GGGG‑A‑CGTAAAACCCGCAAGGGGGG‑AGCGGATAAGAAGTCAACCCGTGACCCTTCTTC‑GCTTACCTCCTCCCGGTGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGGCTTTAATTTACTGATgggc                             >  1:849572/1‑148 (MQ=255)
                                                                                                                gCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                          <  1:1197946/151‑1 (MQ=255)
                                                                                                                      gggCGACTGGC‑GGGT‑A‑CGTAAAAACCGCAAAGGGGG‑GCCGGGTAAAGAGGCAAACCGTGGCCATTCTTT‑GCTGGCATCATCCCGGGGGGTGGCCCCAACACTTAAAAGGGGGGGGGGGGGGGGCTTTAATTTTCTGATgggccttagg                       >  1:1265569/1‑140 (MQ=255)
                                                                                                                       ggCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                   <  1:408560/151‑1 (MQ=255)
                                                                                                                           aCTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt               <  2:1130843/151‑1 (MQ=255)
                                                                                                                             tGGC‑GGGT‑A‑CGTAAAAACCGCAAGGGGGG‑ACCCGGTAAGTAGGTAACCTGCGGCCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGGGGTGGGGTGTGGTCCTTAAATTTACTGATCGCCTTCCGGtcatc                  >  1:797539/1‑146 (MQ=18)
                                                                                                                              ggC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt            <  2:428753/151‑1 (MQ=255)
                                                                                                                               gC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg           <  1:784059/151‑1 (MQ=255)
                                                                                                                               gC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg           <  1:1254423/151‑1 (MQ=255)
                                                                                                                                   ggT‑A‑CGTAAAACCTGCAAGGTGTG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa        <  1:707042/151‑1 (MQ=255)
                                                                                                                                            aaaaCCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGAGCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt  <  1:1063727/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAT  >  NC_000913/266208‑266494

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