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

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584938‑3585236
                                                                                                                                                         |                                                                                                                                                       
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                                                       <  1:899213/150‑1 (MQ=255)
  ttGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                                                    >  2:838665/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 >  1:438153/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 >  1:758750/1‑151 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 >  2:726494/1‑151 (MQ=255)
        aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                              >  2:383720/1‑151 (MQ=255)
          gACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTa                                                                                                                                             >  2:497175/1‑150 (MQ=255)
           aCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTaa                                                                                                                                            >  1:230519/1‑150 (MQ=255)
            cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAat                                                                                                                                          >  1:168904/1‑151 (MQ=255)
            cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAat                                                                                                                                          >  1:597213/1‑151 (MQ=255)
             aGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatc                                                                                                                                         >  2:346404/1‑151 (MQ=255)
                    aaaTTTATTCCTCGGGGAACGCAAAACCCGTGGAAAAAATATCTACAGGATTACAAAAAAAAAAAGCCTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACAAAGAGGGCCGCAACCTTTAGTCCCCCGCCGTTAAATCATCCAc                                                                                                                                   <  1:223934/151‑1 (MQ=255)
                    aaaTTAATTCTCGGGGAACGGAAAACCCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCTTGAAAAAATAGCCCTCTTCCCACGAAGGGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCAc                                                                                                                                   <  1:814069/150‑1 (MQ=255)
                           ttcTCGGGAAACGGAAAACCCGTGGGAAAATTATCAACAGGATTACAAAAAAAAAAGGGCTACCTGCGGCTTCCCATGAAAAAATAGCCCTCTCCCCACAAAGAGGCCCGCTAACCTTATTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  2:570615/148‑1 (MQ=255)
                               ctcgggAAAGGCAAAAACCGCTGAAAAAATATATACAGAATGAGAAAAAAAAAAAGGCTCCCTCGGCCTTCCCCGGAAAAAACAGCCCTCTCCCCACAAAGAGGGCCGTTACTCTTAGACCCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                            <  1:756847/145‑1 (MQ=255)
                                  gggAACGAAAAACCCGCTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      <  2:971847/150‑1 (MQ=255)
                                   ggAAAGGCAAAAACACTGGAAAAATAAATAAAAGATTAAGAACTAAAACAAGCTACCTTAGGATTGCACAAACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     >  2:812170/1‑150 (MQ=255)
                                     aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                  >  2:751600/1‑151 (MQ=255)
                                            aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  1:801782/1‑151 (MQ=255)
                                             aacccGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTa                                                                                                           <  2:538200/148‑1 (MQ=255)
                                             aaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:968869/1‑150 (MQ=255)
                                               accgggggaaagaaTCTCAACAGAATCAAAAAAAAAAAAGGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTTTCCCCACGAAGGGGCCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                          <  2:66366/142‑1 (MQ=255)
                                                  cgtggaaaaaaTATTTAAAGGATTACAAAAAAAAAAAGCCTCCCTCCGGCTTCCCCGAAAAAAATCCCCCTCCCCCCACGAAGGGGGCGGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                          <  1:1209433/139‑1 (MQ=255)
                                                  cgcggaaagaaTATCTGGGGGATAAAAAAAAAAAAAAGCCTCCCTTCGGCTTGCCCAAAAAAAATAGCCCTCTCCCCGCGAAGAGGCCCGCTAACCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                       <  2:768457/145‑1 (MQ=255)
                                                        aaaaaTATCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTTCCCCTAACAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:1054674/147‑1 (MQ=255)
                                                         agaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  2:1089340/1‑149 (MQ=255)
                                                         agaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                >  2:56525/1‑151 (MQ=255)
                                                           aaTATCTACAGTATTACAAAAAAAAAAAGGTTACTTCCGGCTACCCCAGACAAAAGCCCCCCCCCCCCACGAGGGGGCCCCCAAACTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              <  1:487740/151‑1 (MQ=255)
                                                           aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              >  1:739708/1‑151 (MQ=255)
                                                           aaTATCTAAAGGATAACAAAAAAAAAAAGGCTCCCTCCGCCTTCCCCTGAAAAAATAGCCCTCTCCCCACAAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  2:595794/150‑1 (MQ=255)
                                                                  aCAGGATAACAAAAAAAAAAGGCCTCCCTTGGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        <  1:742444/150‑1 (MQ=255)
                                                                  aCAGAATAACAAAAAAAAAAAGGCTCCCTTCGCCTTCCCCTAACAAAATCCCCCCCTCCCCAGGAGGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         <  1:491140/149‑1 (MQ=255)
                                                                     ggATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  2:816612/1‑150 (MQ=255)
                                                                        ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  >  2:162853/1‑150 (MQ=255)
                                                                         tACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                >  2:858050/1‑151 (MQ=255)
                                                                         tACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                >  2:507161/1‑151 (MQ=255)
                                                                          acaaaaaaaaaaaGCCTCCCTTCGGCTCCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGCCCGCAAACTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:788223/148‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGCCTCCCTGCGTCTGCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCCCAAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:971941/148‑1 (MQ=255)
                                                                           cGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              >  1:1201697/1‑151 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:218551/1‑151 (MQ=255)
                                                                              aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  1:938256/1‑151 (MQ=255)
                                                                              aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  2:58382/1‑151 (MQ=255)
                                                                              aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  1:1140274/1‑151 (MQ=255)
                                                                                aaaaaaaGGCCTCCCTCCGCCTTGCCCTAAAAAAATACCCCTCTCCCCACGAGGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                          <  1:3069/150‑1 (MQ=255)
                                                                                aaaaaaaGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGGCCGCTACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:1050364/150‑1 (MQ=255)
                                                                                aaaaaaaGCCTCCCTCCGGCTTCCCCAAAAAAAATCCCCCTCTCCCCAGAAAGGGGCCCGCAACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  1:15687/150‑1 (MQ=255)
                                                                                  aaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                       >  2:1177016/1‑151 (MQ=255)
                                                                                   aaaaGGCTACCTCCGGCTTGCCTTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCTCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTa                                                                      <  2:28593/151‑1 (MQ=255)
                                                                                    aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     >  1:548722/1‑151 (MQ=255)
                                                                                          cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTGCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:673866/1‑151 (MQ=255)
                                                                                               ttgggtttccccgaaCAAAATGCCCCTCTTCCCACAAGGAGGCCCGCTACCTTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                <  2:1160255/132‑1 (MQ=255)
                                                                                                  ggCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                            <  2:438154/147‑1 (MQ=255)
                                                                                                     ttccccTAAAAAAATCGCCCTCTCCCCCCGAAGAGGGCCGCTCCCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                       <  1:1140110/146‑1 (MQ=255)
                                                                                                      tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  1:21232/1‑151 (MQ=255)
                                                                                                       ccctaaaaaaaaaaCCCCTCTCCCCAGGAAGAGGGCCACCACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    <  2:162581/142‑1 (MQ=255)
                                                                                                        cccTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  2:184220/1‑149 (MQ=255)
                                                                                                        cccTAAAAAAATGCCCCTCTCCCCACGAAGGGGCCCTCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:33341/144‑1 (MQ=255)
                                                                                                        cccTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  2:886973/148‑1 (MQ=255)
                                                                                                           tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               >  2:861328/1‑151 (MQ=255)
                                                                                                           tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               >  1:1138933/1‑151 (MQ=255)
                                                                                                              cAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                            >  2:652497/1‑151 (MQ=255)
                                                                                                                aaaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          >  2:1075316/1‑151 (MQ=255)
                                                                                                                 aaatgCCCCCCTCCCCCGGAGGGGGCCGGCTAACCTATTACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          <  2:582844/147‑1 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         >  2:145968/1‑151 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         >  1:444854/1‑151 (MQ=255)
                                                                                                                    aGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                                      >  1:261907/1‑151 (MQ=255)
                                                                                                                      cccTCTCCCCACGAAGGGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGa                                     <  2:953842/150‑1 (MQ=255)
                                                                                                                               cACGAAGGGGCCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  2:351169/150‑1 (MQ=255)
                                                                                                                                   aaGAGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                       <  2:168904/151‑1 (MQ=255)
                                                                                                                                     gggggCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           <  1:838662/143‑1 (MQ=255)
                                                                                                                                     gagGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca                     <  1:1172730/151‑1 (MQ=255)
                                                                                                                                           cGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa               <  1:381868/151‑1 (MQ=255)
                                                                                                                                           cGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa               <  2:490436/151‑1 (MQ=255)
                                                                                                                                                aaCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt           >  1:1031451/1‑151 (MQ=255)
                                                                                                                                                  ccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             <  1:1257747/147‑1 (MQ=255)
                                                                                                                                                     tAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       <  2:729634/150‑1 (MQ=255)
                                                                                                                                                      aGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCgg     >  1:727038/1‑151 (MQ=255)
                                                                                                                                                        tACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCggg    >  1:556037/1‑150 (MQ=255)
                                                                                                                                                         aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  2:52095/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  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: