Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC41.6% 132.2 / 29.1 78Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (40/5);  new base C (0/32);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 9.45e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGATAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584938‑3585236
                                                                                                                                                         |                                                                                                                                                       
tattGATAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑c                                                                                                                                                       <  1:899213/150‑1 (MQ=255)
  ttGATAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑cccc                                                                                                                                                    >  2:838665/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCccgcc                                                                                                                                                 >  1:438153/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCccgcc                                                                                                                                                 >  1:758750/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCccgcc                                                                                                                                                 >  2:726494/1‑151 (MQ=255)
        aaGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGtt                                                                                                                                              >  2:383720/1‑151 (MQ=255)
          gACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTa                                                                                                                                             >  2:497175/1‑150 (MQ=255)
           aCAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTaa                                                                                                                                            >  1:230519/1‑150 (MQ=255)
            cAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAat                                                                                                                                          >  1:168904/1‑151 (MQ=255)
            cAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAat                                                                                                                                          >  1:597213/1‑151 (MQ=255)
             aGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAatc                                                                                                                                         >  2:346404/1‑151 (MQ=255)
                    aaaTTTATTCCTCGGGGAAC‑GC‑AAAACCCGTGGAAAAAATATCTACAGGATTACAAAAAAAAAAAGC‑CTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACAAAGAGGGCC‑GCAACCTTTAGTC‑CCCCGCCGTTAAATCATCCAc                                                                                                                                   <  1:223934/151‑1 (MQ=255)
                    aaaTTAATT‑CTCGGGGAACGGA‑AAACCCGCTGG‑AAAAATATCTACAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTCCCTTGAAAAAATAGCCCTCTTCCCACGAAGGGGGCC‑GCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCAc                                                                                                                                   <  1:814069/150‑1 (MQ=255)
                           tt‑cTCGGGAAACGGA‑AAACCCGTGGG‑AAAATTATCAACAGGATTACAAAAAAAAAAGGG‑CTACCTGCGGCTTCCCATGAAAAAATAGCCCTCTCCCCACAAAGAGGCCC‑GCTAACCTTATTA‑CCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  2:570615/148‑1 (MQ=255)
                               ctcgggAAAGGCAAAAACCGCTGA‑AAAAATATATACAGAATGAGAAAAAAAAAAAGG‑CTCCCTCGGCCTTCCCCGGAAAAAACAGCCCTCTCCCCACAAAGAGGGCC‑GTTACTCTTAGAC‑CCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                            <  1:756847/145‑1 (MQ=255)
                                  gggAACGAA‑AAACCCGCTGA‑AAAAATATCTACAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTCCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGGCC‑GCTACCCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      <  2:971847/150‑1 (MQ=255)
                                   ggAAAGGC‑AAAAACACTGG‑AAAAATAAATAAAAGATTAAGAACTAAAACAAG‑CTACCTTAGGATTGCACAAACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     >  2:812170/1‑150 (MQ=255)
                                     aaCGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                  >  2:751600/1‑151 (MQ=255)
                                            aaaaCCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  1:801782/1‑151 (MQ=255)
                                             aacccGCTGG‑AAAAATATCTACAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCC‑GCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTa                                                                                                           <  2:538200/148‑1 (MQ=255)
                                             aaaCCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:968869/1‑150 (MQ=255)
                                               accgggggaaagaaTCTCAACAGAATCAAAAAAAAAAAAGGG‑CTCCCTTCGGCTTCCCCTGAAAAAATACCCCTTTCCCCACGAAGGGGCCC‑GCTACCCTAAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                          <  2:66366/142‑1 (MQ=255)
                                                  cgtggaaaaaaTATTTAAAGGATTACAAAAAAAAAAAGC‑CTCCCTCCGGCTTCCCCGAAAAAAATCCCCCTCCCCCCACGAAGGGGGCG‑GCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                          <  1:1209433/139‑1 (MQ=255)
                                                  cgcggaaagaaTATCTGGGGGATAAAAAAAAAAAAAAGC‑CTCCCTTCGGCTTGCCCAAAAAAAATAGCCCTCTCCCCGCGAAGAGGCCC‑GCTAACCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                       <  2:768457/145‑1 (MQ=255)
                                                        aaaaaTATCTACAGGATTACAAAAAAAAAAAGG‑CTCCCTCCGGCTTCCCCTAACAAAATACCCCTCTCCCCACGAAGGGGGCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:1054674/147‑1 (MQ=255)
                                                         agaaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  2:1089340/1‑149 (MQ=255)
                                                         agaaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                >  2:56525/1‑151 (MQ=255)
                                                           aaTATCTACAGTATTACAAAAAAAAAAAGG‑TTACTTCCGGCTACCCCAGACAAAAGCCCCCCCCCCCCACGAGGGGGCCC‑CCAAACTTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              <  1:487740/151‑1 (MQ=255)
                                                           aaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              >  1:739708/1‑151 (MQ=255)
                                                           aaTATCTAAAGGATAACAAAAAAAAAAAGG‑CTCCCTCCGCCTTCCCCTGAAAAAATAGCCCTCTCCCCACAAAGAGGCCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  2:595794/150‑1 (MQ=255)
                                                                  aCAGGATAACAAAAAAAAAAGGC‑CTCCCTTGGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        <  1:742444/150‑1 (MQ=255)
                                                                  aCAGAATAACAAAAAAAAAAAGG‑CTCCCTTCGCCTTCCCCTAACAAAATCCCCCCCTCCCCAGGAGGGGGCCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         <  1:491140/149‑1 (MQ=255)
                                                                     ggATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  2:816612/1‑150 (MQ=255)
                                                                        ttACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  >  2:162853/1‑150 (MQ=255)
                                                                         tACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                >  2:858050/1‑151 (MQ=255)
                                                                         tACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                >  2:507161/1‑151 (MQ=255)
                                                                          acaaaaaaaaaaaGC‑CTCCCTTCGGCTCCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGCCC‑GCAAACTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:788223/148‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGC‑CTCCCTGCGTCTGCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCC‑CCAAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:971941/148‑1 (MQ=255)
                                                                           cGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              >  1:1201697/1‑151 (MQ=255)
                                                                             aaaTAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:218551/1‑151 (MQ=255)
                                                                              aaTAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  1:938256/1‑151 (MQ=255)
                                                                              aaTAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  2:58382/1‑151 (MQ=255)
                                                                              aaTAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  1:1140274/1‑151 (MQ=255)
                                                                                aaaaaaaGGCCTCCCTCCGCCTTGCCCTAAAAAAATACCCCTCTCCCCACGAGGGGGCCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                          <  1:3069/150‑1 (MQ=255)
                                                                                aaaaaaaGG‑CTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGGCC‑GCTACCCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:1050364/150‑1 (MQ=255)
                                                                                aaaaaaaGC‑CTCCCTCCGGCTTCCCCAAAAAAAATCCCCCTCTCCCCAGAAAGGGGCCC‑GCAACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  1:15687/150‑1 (MQ=255)
                                                                                  aaaaaGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                       >  2:1177016/1‑151 (MQ=255)
                                                                                   aaaaGG‑CTACCTCCGGCTTGCCTTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCC‑TCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTa                                                                      <  2:28593/151‑1 (MQ=255)
                                                                                    aaaGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     >  1:548722/1‑151 (MQ=255)
                                                                                          cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTG‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:673866/1‑151 (MQ=255)
                                                                                               ttgggtttccccgaaCAAAATGCCCCTCTTCCCACAAGGAGGCCC‑GCTACCTTTATTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                <  2:1160255/132‑1 (MQ=255)
                                                                                                  ggCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCC‑GCTACCCTTATTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                            <  2:438154/147‑1 (MQ=255)
                                                                                                     ttccccTAAAAAAATCGCCCTCTCCCCCCGAAGAGGGCC‑GCTCCCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                       <  1:1140110/146‑1 (MQ=255)
                                                                                                      tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  1:21232/1‑151 (MQ=255)
                                                                                                       ccctaaaaaaaaaaCCCCTCTCCCCAGGAAGAGGGCC‑ACCACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    <  2:162581/142‑1 (MQ=255)
                                                                                                        cccTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  2:184220/1‑149 (MQ=255)
                                                                                                        cccTAAAAAAATGCCCCTCTCCCCACGAAGGGGCCC‑TCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:33341/144‑1 (MQ=255)
                                                                                                        cccTAAAAAAATACCCCTCTCCCCACGAAGGGGGCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  2:886973/148‑1 (MQ=255)
                                                                                                           tGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               >  2:861328/1‑151 (MQ=255)
                                                                                                           tGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               >  1:1138933/1‑151 (MQ=255)
                                                                                                              cAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                            >  2:652497/1‑151 (MQ=255)
                                                                                                                aaaTAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          >  2:1075316/1‑151 (MQ=255)
                                                                                                                 aaatgCCCCCCTCCCCCGGAGGGGGCCGGCTAACC‑TATTACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          <  2:582844/147‑1 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         >  2:145968/1‑151 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         >  1:444854/1‑151 (MQ=255)
                                                                                                                    aGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                                      >  1:261907/1‑151 (MQ=255)
                                                                                                                      cccTCTCCCCACGAAGGGGCCC‑GCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGa                                     <  2:953842/150‑1 (MQ=255)
                                                                                                                               cACGAAGGGGCCC‑GCTACCCTTATTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  2:351169/150‑1 (MQ=255)
                                                                                                                                   aaGAGGCCC‑GCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                       <  2:168904/151‑1 (MQ=255)
                                                                                                                                     gggggCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           <  1:838662/143‑1 (MQ=255)
                                                                                                                                     gagGGCC‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca                     <  1:1172730/151‑1 (MQ=255)
                                                                                                                                           c‑GCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa               <  1:381868/151‑1 (MQ=255)
                                                                                                                                           c‑GCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa               <  2:490436/151‑1 (MQ=255)
                                                                                                                                                aaCCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt           >  1:1031451/1‑151 (MQ=255)
                                                                                                                                                  ccTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             <  1:1257747/147‑1 (MQ=255)
                                                                                                                                                     tAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       <  2:729634/150‑1 (MQ=255)
                                                                                                                                                      aGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCgg     >  1:727038/1‑151 (MQ=255)
                                                                                                                                                        tA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCggg    >  1:556037/1‑150 (MQ=255)
                                                                                                                                                         a‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  2:52095/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
TATTGATAAAGACAGTGGTGAAATTCATT‑CTCTGGGAACGGC‑AAAACCGCTGG‑AAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCC‑GCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584938‑3585236

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