Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3580TG44.6% 83.0 / 23.2 67T66P (ACC→CCC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (4/32);  new base G (29/0);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.83e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506
                                                                                                                                                   |                                                                                                                                                      
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                                                       <  1:200497/150‑1 (MQ=255)
                                                                                          aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                       <  2:635894/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:1254423/151‑1 (MQ=255)
                                                                                                                         gC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                       <  1:784059/151‑1 (MQ=255)
                                                                                                                             ggT‑A‑CGTAAAACCTGCAAGGTGTG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                    <  1:707042/151‑1 (MQ=255)
                                                                                                                                      aaaaCCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGAGCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt              <  1:1063727/150‑1 (MQ=255)
                                                                                                                                             gCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg        >  2:1216246/1‑149 (MQ=255)
                                                                                                                                                aGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaaa   <  1:907992/151‑1 (MQ=255)
                                                                                                                                                aGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa    >  1:1332301/1‑150 (MQ=255)
                                                                                                                                                  gTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:560895/1‑150 (MQ=255)
                                                                                                                                                   tGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:254591/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGT‑A‑CGTAAAACCTGCAAGGTGGG‑ACCGGATAAGTAGTCAACCTGTGACCATTCTTC‑GCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506

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: ‑