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

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236
                                                                                                                                                         |                                                                                                                                                       
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                                                         <  2:412426/151‑1 (MQ=255)
       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 <  1:488550/151‑1 (MQ=255)
             aCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTaa                                                                                                                                            >  1:626500/1‑150 (MQ=255)
                gtggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatc                                                                                                                                         >  1:814119/1‑150 (MQ=255)
                       aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGCTTCCGAAATAAAAAAGGCTACCCTCGGCTTGCCTTGACAAAATAGCCCTCTTCCCACGAAGACGGCCGCTACCCTTAGTACCCCGCCGTTTAATCATCCAcc                                                                                                                                  >  2:144983/1‑150 (MQ=255)
                             ttcTCTGGGACCGGCAAACCCGCTGGAAAAATTTCTACAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTCCCCCTAAAAAAATCCCCCTCTTCCCCCGAAGAGGCCCGCAACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  1:1182900/148‑1 (MQ=255)
                             ttcTCGGGGAACGGCAAAACCGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTGCGGCTGCCCCTGACAAAATCCCCCTCTCCCCACGAGGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcg                                                                                                                             <  1:216301/149‑1 (MQ=255)
                             tcctcgGGAAACGACAAACCCGCTGAAAAAATATCAAAAGTATTACAAAAAAAAAACGCCTCCCTTCGCCTTCCCCTGACAAAATACCCCTTTCCCCAAAAAGAGGGCCTCTAACCTTAGCCCCCCGCCGTTAAATCATCCACCGa                                                                                                                                <  2:567117/144‑1 (MQ=255)
                                    ggatcaacAAACCCGTAGGAAGAATAGCAATAGATTTCCAAAACAAAAAGGGCTACCTCGATTATTCCTTGACAAAACACCCCTTTTCCCAGGAGGAGGGCCGCAAACCTAAAAACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          <  1:1196772/138‑1 (MQ=255)
                                      aaCGGAAAAACCGCTGGAAGAATATCTACAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:217791/150‑1 (MQ=255)
                                      aaCGGAAAAACCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAGGGCTCCCTTCGCCTTCCCCGGAAAAAATGGCCCTCTTCCCACGAAGGGGGCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:96181/150‑1 (MQ=255)
                                         ggaaaacCCGCTGAAAAAATATCTAAAGAATTCCAAAAAAAAAAAGGCTCCTTTCGCCTTCCCCAAAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCAACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                   <  1:89762/144‑1 (MQ=255)
                                            aaacCCGCTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            <  1:341715/151‑1 (MQ=255)
                                             aaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:565597/1‑150 (MQ=255)
                                                  acTGGAAAAATACCTACAGAATAACAAAAAAAAAAAGGCTACTTCCGGCTGCCCCTAACAAAGTAGCCCTCTCCCCACGAAGGGGGCCGCTCACCTTACTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                       <  2:65086/149‑1 (MQ=255)
                                                     ggAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  1:427322/1‑151 (MQ=255)
                                                     gaaaaaaTATTTAAAGGATAAAAAAAAAAAAAAGGCGCCCGTCGGTTTCCCCTAACAAAATAGCCCTCTTCCCACAAAGGGCGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                    <  1:827066/145‑1 (MQ=255)
                                                     gaaaaaaTATCTCCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     <  1:1354050/144‑1 (MQ=255)
                                                      gaagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  2:251593/1‑151 (MQ=255)
                                                       aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCATCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  2:518817/1‑149 (MQ=255)
                                                       aagaaTATCACAAGATTTACAAAAAAAAAAGGCCTCCCTCCGCCTGCCCCGGAAAAAATCCCCTTCTCCCCCCGAGAGGGCCCGCTACCCTTTGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:425758/150‑1 (MQ=255)
                                                        agaaTATCTACAGGAATACAAAAAAAAAAAGCCTACCTGCCGGTCGCCCTAAAAAATCCCCCCTCTCCCCAGAAGGTGGCCCGCTACCCAAAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:123096/149‑1 (MQ=255)
                                                        aaaaTATCAACAGGATTACAAAAAAAAAAAGGCTCCCTCCTTCTTCCCCTGACAAAACCGCCCTCCTCCCACAAGGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  2:355305/149‑1 (MQ=255)
                                                        aaaaTATCAACAGAATAACAAAAAAAAAAAGCCTCCCTCGGGCTTCCCCAGAAAAAATACCCCTCTCCCCACGAGGAGGGCCGCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  2:722355/149‑1 (MQ=255)
                                                          aattcatACAGGATTACAAAAAAAAAAAGGCTCCCTGCGCCTTGCCCTAAAAAAATGCCCCTCTCCCCGCGAGGGGGCCGGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              <  2:800635/145‑1 (MQ=255)
                                                          aaTATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGCCTGCCCCTAAAAAAATACCCCTCTCCCCACGAGGAGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  1:869632/150‑1 (MQ=255)
                                                            tatCTACGGTATCACAAAATAAAAATGGCTCCCTGGGGCTCCCCTTAAAAAAACAGCCCTCTCCCCAGGAAGAGCGCCAATATCCGTGCAACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  1:846604/150‑1 (MQ=255)
                                                                   aGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  2:958651/1‑151 (MQ=255)
                                                                        taaaaaaaaaaaaaaGGCTCCCTTGGGCTGGCCCGAAAAAAATAGCCCTCTCCCCACAAAGAGGCCCTCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  <  1:697293/145‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGGCTCCCTTCGGCTTGCCCTAAAAAAATAGCCCTCTCCCCACGAAGGGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  1:880564/148‑1 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:788499/1‑151 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCACTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCAACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  2:751875/1‑151 (MQ=255)
                                                                                aaaaaagGGCTACCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCAGGAAGGGGGCCGCTACCTTGAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:587058/150‑1 (MQ=255)
                                                                                aaaaaaaGCCTCCCTGCGCCTTGCCCTAAAAAAACACCCCTCTCCCCACGAAGAGGGCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:191758/150‑1 (MQ=255)
                                                                                   aaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                       >  2:469323/1‑150 (MQ=255)
                                                                                    aaagcctccttgcgccTTCCCCTAAAAAAATACCCCTCTCCCCCGAAAGGGGCCCTCAAACTTTGCACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     <  1:528601/136‑1 (MQ=255)
                                                                                     aaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACt                                                                    >  1:1428356/1‑151 (MQ=255)
                                                                                       gcctccttgcGGCTTCCCCTAAAAAAATCCCCCTCTCCCCCCGAAGGGGGCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                                   <  2:1302921/141‑1 (MQ=255)
                                                                                         ctcccTTCGCCTTCCCCTGACAAAATACCCCTTTCCCCTCGAAGGGGGCCGTTAACTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  <  1:142595/146‑1 (MQ=255)
                                                                                            ccTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCAGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                             >  1:1239553/1‑151 (MQ=255)
                                                                                             cTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                            >  2:468353/1‑151 (MQ=255)
                                                                                              ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:392127/1‑151 (MQ=255)
                                                                                              ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:520058/1‑151 (MQ=255)
                                                                                                cGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         >  1:1108804/1‑151 (MQ=255)
                                                                                                 ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  2:767268/1‑151 (MQ=255)
                                                                                                 ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  1:457405/1‑151 (MQ=255)
                                                                                                         ccTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                                 >  1:1262088/1‑151 (MQ=255)
                                                                                                          cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                                >  2:1264373/1‑151 (MQ=255)
                                                                                                             aaaaaaTCCCCCTCTCCCCACGGAGGGCCCCGTTACCCTACGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               <  2:1170230/147‑1 (MQ=255)
                                                                                                             aCAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                              >  1:963244/1‑150 (MQ=255)
                                                                                                             aCAAAATACCCCACTTCCCACGAAGAGGCCCGCTATCGTTAGTCCGCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                               <  1:326846/149‑1 (MQ=255)
                                                                                                               aaaaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                           >  2:1358139/1‑151 (MQ=255)
                                                                                                                 aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                          >  2:1078837/1‑150 (MQ=255)
                                                                                                                 aaTACCCCTCTCCCCACGAAGAGGCCCGCTAACCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  1:288582/151‑1 (MQ=255)
                                                                                                                      cccTCTCCCCACAAAGAGGCCCCCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:685824/148‑1 (MQ=255)
                                                                                                                       ccTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   >  1:1289806/1‑151 (MQ=255)
                                                                                                                          ctTCCCCGGAAGGGGCCCGTTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGa                                     <  2:1419760/146‑1 (MQ=255)
                                                                                                                             cccAGGAGGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                                <  1:639673/148‑1 (MQ=255)
                                                                                                                             cccACGAGGAGGCCCGCTAACCTTAGCACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                                 <  1:712398/147‑1 (MQ=255)
                                                                                                                              ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:1284598/1‑151 (MQ=255)
                                                                                                                                 cGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:27247/1‑151 (MQ=255)
                                                                                                                                 cGAAGAGGGCCGCTAACCTTAGGACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:234367/1‑151 (MQ=255)
                                                                                                                                   aaGAGGCCCGCTACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          <  1:309074/148‑1 (MQ=255)
                                                                                                                                        ggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                  >  1:1340438/1‑151 (MQ=255)
                                                                                                                                           cccaccccttcaGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  2:582803/137‑1 (MQ=255)
                                                                                                                                           ccGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt                 <  1:496653/150‑1 (MQ=255)
                                                                                                                                            cGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:907283/1‑150 (MQ=255)
                                                                                                                                             gCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg              <  1:1059399/151‑1 (MQ=255)
                                                                                                                                              cTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg              <  2:1325870/150‑1 (MQ=255)
                                                                                                                                              cTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             >  1:1270955/1‑151 (MQ=255)
                                                                                                                                               tAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg            >  2:1369532/1‑151 (MQ=255)
                                                                                                                                                     tAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       <  1:1252209/150‑1 (MQ=255)
                                                                                                                                                         aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:908430/1‑151 (MQ=255)
                                                                                                                                                         aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  2:1420799/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  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: