Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC45.1% 96.6 / 21.4 71Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (32/7);  new base C (0/32);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 9.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236
                                                                                                                                                         |                                                                                                                                                       
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTa                                                                                                                                                         <  2:412426/151‑1 (MQ=255)
       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCccgcc                                                                                                                                                 <  1:488550/151‑1 (MQ=255)
             aCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTaa                                                                                                                                            >  1:626500/1‑150 (MQ=255)
                gtggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAatc                                                                                                                                         >  1:814119/1‑150 (MQ=255)
                       aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGC‑TTCCGAAATAAAAAAGGCTACCCTCGGCT‑TGCCTTGACAAAATAGCCCTCTTCCCACGAAGACGG‑CCGCTACCCTTAGTA‑CCCCGCCGTTTAATCATCCAcc                                                                                                                                  >  2:144983/1‑150 (MQ=255)
                             ttcTCTGGGACCGGCAAACCCGCTGGAAAAATTTCTACAGGA‑TTAAAAAAAAAAAAAGGCTCCCTTCGGCT‑CCCCCTAAAAAAATCCCCCTCTTCCCCCGAAGAGGC‑CCGCAACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  1:1182900/148‑1 (MQ=255)
                             ttcTCGGGGAACGGCAAAACCGCTGAAAGAATATCTACAGGA‑TTACAAAAAAAAAAAGGCTCCCTGCGGCT‑GCCCCTGACAAAATCCCCCTCTCCCCACGAGGGGGC‑CCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAgcg                                                                                                                             <  1:216301/149‑1 (MQ=255)
                             tcctcgGGAAACGACAAACCCGCTGAAAAAATATCAAAAGTA‑TTACAAAAAAAAAACGCCTCCCTTCGCCT‑TCCCCTGACAAAATACCCCTTTCCCCAAAAAGAGGG‑CCTCTAACCTTAGCC‑CCCCGCCGTTAAATCATCCACCGa                                                                                                                                <  2:567117/144‑1 (MQ=255)
                                    ggatcaacAAACCCGTAGGAAGAATAGCAATAGAT‑TTCCAAAACAAAAAGGGCTACC‑TCGATTATTCCTTGACAAAACACCCCTTTTCCCAGGAGGAGGG‑CCGCAAACCTAAAAA‑CCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          <  1:1196772/138‑1 (MQ=255)
                                      aaCGGAAAAACCGCTGGAAGAATATCTACAGGA‑TTAAAAAAAAAAAAAGGCTCCCTTCGGCT‑TGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGG‑CCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:217791/150‑1 (MQ=255)
                                      aaCGGAAAAACCGCTGGAAAAATATCTACAGGA‑TTACAAAAAAAAAAGGGCTCCCTTCGCCT‑TCCCCGGAAAAAATGGCCCTCTTCCCACGAAGGGGG‑CCGCTACCCTAAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:96181/150‑1 (MQ=255)
                                         ggaaaacCCGCTGAAAAAATATCTAAAGAA‑TTCCAAAAAAAAAAAGGCTCCTTTCGCCT‑TCCCCAAAAAAAATCCCCCTCTCCCCACGAAGAGGC‑CCGCAACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:89762/144‑1 (MQ=255)
                                            aaacCCGCTGAAAAAATATCTACAGGA‑TTACAAAAAAAAAAAGGCTCCCTTCGGCT‑TGCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGG‑CCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            <  1:341715/151‑1 (MQ=255)
                                             aaaCCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:565597/1‑150 (MQ=255)
                                                  acTGGAAAAATACCTACAGAA‑TAACAAAAAAAAAAAGGCTACTTCCGGCT‑GCCCCTAACAAAGTAGCCCTCTCCCCACGAAGGGGG‑CCGCTCACCTTACTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                       <  2:65086/149‑1 (MQ=255)
                                                     ggAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  1:427322/1‑151 (MQ=255)
                                                     gaaaaaaTATTTAAAGGA‑TAAAAAAAAAAAAAAGGCGCCCGTCGGTT‑TCCCCTAACAAAATAGCCCTCTTCCCACAAAGGGCG‑CCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                    <  1:827066/145‑1 (MQ=255)
                                                     gaaaaaaTATCTCCAGGA‑TTACAAAAAAAAAAAGGCTCCCTTCGGCT‑TCCCCTGACAAAATACCCCTCTCCCCACGAAGAGGC‑CCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     <  1:1354050/144‑1 (MQ=255)
                                                      gaagaaTATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  2:251593/1‑151 (MQ=255)
                                                       aagaaTATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCATCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  2:518817/1‑149 (MQ=255)
                                                       aagaaTATC‑ACAAGATTTACAAAAAAAAAAGGCCTCCCTCCGCCT‑GCCCCGGAAAAAATCCCCTTCTCCCCCCG‑AGAGGGCCCGCTACCCTTTGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:425758/150‑1 (MQ=255)
                                                        agaaTATCTACAGGA‑ATACAAAAAAAAAAAGCCTACCTGCCGGT‑CGCCCTAAAAAATCCCCCCTCTCCCCAGAAGGTGGC‑CCGCTACCCAAAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:123096/149‑1 (MQ=255)
                                                        aaaaTATCAACAGGA‑TTACAAAAAAAAAAAGGCTCCCTCCTTCT‑TCCCCTGACAAAACCGCCCTCCTCCCACAAGGGGGG‑CCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  2:355305/149‑1 (MQ=255)
                                                        aaaaTATCAACAGAA‑TAACAAAAAAAAAAAGCCTCCCTCGGGCT‑TCCCCAGAAAAAATACCCCTCTCCCCACGAGGAGGG‑CCGCAACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  2:722355/149‑1 (MQ=255)
                                                          aattcatACAGGA‑TTACAAAAAAAAAAAGGCTCCCTGCGCCT‑TGCCCTAAAAAAATGCCCCTCTCCCCGCGAGGGGGC‑CGGCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              <  2:800635/145‑1 (MQ=255)
                                                          aaTATCTACAGGA‑TTACAAAAAAAAAAAGGCTCCCTTCGCCT‑GCCCCTAAAAAAATACCCCTCTCCCCACGAGGAGGC‑CCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  1:869632/150‑1 (MQ=255)
                                                            tatCTACGGTA‑TCACAAAATAAAAATGGCTCCCTGGGGCT‑CCCCTTAAAAAAACAGCCCTCTCCCCAGGAAGAGCG‑CCAATATCCGTGCAA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  1:846604/150‑1 (MQ=255)
                                                                   aGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  2:958651/1‑151 (MQ=255)
                                                                        taaaaaaaaaaaaaaGGCTCCCTTGGGCT‑GGCCCGAAAAAAATAGCCCTCTCCCCACAAAGAGGC‑CCTCAACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  <  1:697293/145‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGGCTCCCTTCGGCT‑TGCCCTAAAAAAATAGCCCTCTCCCCACGAAGGGGC‑CCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  1:880564/148‑1 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:788499/1‑151 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCACTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCAACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  2:751875/1‑151 (MQ=255)
                                                                                aaaaaagGGCTACCTTCGGCT‑TCCCCTGACAAAATACCCCTCTCCCCAGGAAGGGGG‑CCGCTACCTTGAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:587058/150‑1 (MQ=255)
                                                                                aaaaaaaGCCTCCCTGCGCCT‑TGCCCTAAAAAAACACCCCTCTCCCCACGAAGAGGG‑CCCCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:191758/150‑1 (MQ=255)
                                                                                   aaaaGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                       >  2:469323/1‑150 (MQ=255)
                                                                                    aaagcctccttgcgccT‑TCCCCTAAAAAAATACCCCTC‑TCCCCCGAAAGGGG‑CCCTCAAACTTTGCACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     <  1:528601/136‑1 (MQ=255)
                                                                                     aaGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACt                                                                    >  1:1428356/1‑151 (MQ=255)
                                                                                       gcctccttgcGGCT‑TCCCCTAAAAAAATCCCCCTCTCCCCCCGAAGGGGG‑CCCCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                                   <  2:1302921/141‑1 (MQ=255)
                                                                                         ctcccTTCGCCT‑TCCCCTGACAAAATACCCCTTTCCCCTCGAAGGGGG‑CCGTTAACTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  <  1:142595/146‑1 (MQ=255)
                                                                                            ccTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCAGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                             >  1:1239553/1‑151 (MQ=255)
                                                                                             cTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                            >  2:468353/1‑151 (MQ=255)
                                                                                              ttCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:392127/1‑151 (MQ=255)
                                                                                              ttCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:520058/1‑151 (MQ=255)
                                                                                                cGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         >  1:1108804/1‑151 (MQ=255)
                                                                                                 ggCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  2:767268/1‑151 (MQ=255)
                                                                                                 ggCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  1:457405/1‑151 (MQ=255)
                                                                                                         ccTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                                 >  1:1262088/1‑151 (MQ=255)
                                                                                                          cTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                                >  2:1264373/1‑151 (MQ=255)
                                                                                                             aaaaaaTCCCCCTCTCCCCACGGAGGGCC‑CCGTTACCCTACGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               <  2:1170230/147‑1 (MQ=255)
                                                                                                             aCAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                              >  1:963244/1‑150 (MQ=255)
                                                                                                             aCAAAATACCCCACTTCCCACGAAGAGGC‑CCGCTATCGTTAGTC‑CGCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               <  1:326846/149‑1 (MQ=255)
                                                                                                               aaaaTAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                           >  2:1358139/1‑151 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          >  2:1078837/1‑150 (MQ=255)
                                                                                                                 aaTACCCCTCTCCCCACGAAGAGGC‑CCGCTAACCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  1:288582/151‑1 (MQ=255)
                                                                                                                      cccTCTCCCCACAAAGAGGC‑CCCCTACCCTAAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:685824/148‑1 (MQ=255)
                                                                                                                       ccTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   >  1:1289806/1‑151 (MQ=255)
                                                                                                                          ctTCCCCGGAAGGGGC‑CCGTTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGa                                     <  2:1419760/146‑1 (MQ=255)
                                                                                                                             cccAGGAGGAGGC‑CCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                                <  1:639673/148‑1 (MQ=255)
                                                                                                                             cccACGAGGAGGC‑CCGCTAACCTTAGCA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                                 <  1:712398/147‑1 (MQ=255)
                                                                                                                              ccACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:1284598/1‑151 (MQ=255)
                                                                                                                                 cGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:27247/1‑151 (MQ=255)
                                                                                                                                 cGAAGAGGG‑CCGCTAACCTTAGGA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:234367/1‑151 (MQ=255)
                                                                                                                                   aaGAGGC‑CCGCTACCCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          <  1:309074/148‑1 (MQ=255)
                                                                                                                                        gg‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                  >  1:1340438/1‑151 (MQ=255)
                                                                                                                                           cccaccccttcaGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  2:582803/137‑1 (MQ=255)
                                                                                                                                           ccGCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt                 <  1:496653/150‑1 (MQ=255)
                                                                                                                                            cGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:907283/1‑150 (MQ=255)
                                                                                                                                             gCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg              <  1:1059399/151‑1 (MQ=255)
                                                                                                                                              cTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg              <  2:1325870/150‑1 (MQ=255)
                                                                                                                                              cTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             >  1:1270955/1‑151 (MQ=255)
                                                                                                                                               tAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg            >  2:1369532/1‑151 (MQ=255)
                                                                                                                                                     tAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       <  1:1252209/150‑1 (MQ=255)
                                                                                                                                                         a‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:908430/1‑151 (MQ=255)
                                                                                                                                                         a‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  2:1420799/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGA‑TTACGAAATAAAAAAGGCTACCTTCGGCT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236

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