Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC45.5% 159.7 / 13.8 88Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (43/5);  new base C (0/40);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.74e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236
                                                                                                                                                            |                                                                                                                                                      
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                                                        <  2:790362/151‑1 (MQ=255)
    ttGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                                                    >  2:129011/1‑151 (MQ=255)
       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 >  2:1213420/1‑151 (MQ=255)
        tAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgccg                                                                                                                                                >  1:519577/1‑151 (MQ=255)
         aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgccg                                                                                                                                                <  2:735580/150‑1 (MQ=255)
         aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGt                                                                                                                                               >  1:807908/1‑151 (MQ=255)
             aCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTaa                                                                                                                                            >  2:1349571/1‑150 (MQ=255)
             aCAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTaa                                                                                                                                            >  1:329301/1‑150 (MQ=255)
               aGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatc                                                                                                                                         >  1:778871/1‑151 (MQ=255)
                 tggtgAAATTCATTCTCTGGGAACGGCAAAGCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatcat                                                                                                                                       >  2:83279/1‑151 (MQ=255)
                      aaaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCAcc                                                                                                                                  >  1:768077/1‑151 (MQ=255)
                          tctttcccggggAACGCAAAAGCCGGGGGAAAAATATCTAAAGAATAACAAAAAAAAAAAGATTCCCTTGGCTTTCCCTAAACAAAAGCGCCCCCTACCCACAAAGAGGCCGCCTAACCTAAGTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                               <  2:1456684/141‑1 (MQ=255)
                            attcTCTGGGAACGGAAAAACCGCTGGAAAAATATCTACAGAATAACAAAATAAAAAAGGCTCCCTTGGGCTTGCCCTGAAAAAATGCCCCTCTCCCCACGAAGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                            <  2:575592/151‑1 (MQ=255)
                              ttctcgggAAAGGAAAAACCCGCGGGAAAAATATCTACAGGATCACAAAAAAAAAAAGCCTCCCTTCGCCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAGGAGGCCCGCTAACTTAAGTCCCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                            <  1:1206054/145‑1 (MQ=255)
                               ctctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                         >  1:861048/1‑151 (MQ=255)
                                 cggggAACGAAAAACCCGCGGAAAAAATATCAACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATCCCCCTCCCCCCAAGAGGAGGGCCCATACCTTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                       <  1:370935/149‑1 (MQ=255)
                                  tGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      >  2:445404/1‑151 (MQ=255)
                                    ggtagcgaaaaccccGCTGAAAAAATGTCTAAAATATTAAAAAAAAAAAACGCCTCCCTGCGTCTACCCCTAAAAAATTACCCTTCCCCCCGGGAGGGGGCCCGTACACCTTAGCACGCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     <  1:1382619/138‑1 (MQ=255)
                                       aCGGAAAACCCGCTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  2:770797/151‑1 (MQ=255)
                                        cGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                                >  2:756420/1‑151 (MQ=255)
                                            aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  1:457632/1‑151 (MQ=255)
                                              aacccGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATCCCCCTCTCCCCACAAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTa                                                                                                           <  2:494738/148‑1 (MQ=255)
                                                 ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         >  1:22545/1‑150 (MQ=255)
                                                 ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         >  1:1450100/1‑150 (MQ=255)
                                                     gggAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCAGGAGGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                    <  1:319967/150‑1 (MQ=255)
                                                           aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              >  1:153322/1‑151 (MQ=255)
                                                            atatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             >  1:168892/1‑151 (MQ=255)
                                                            atatCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCAGACAAAATACCCCTCTCCCCAGGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  2:1253225/151‑1 (MQ=255)
                                                            atatCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  2:1326528/151‑1 (MQ=255)
                                                             tatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                            >  1:512800/1‑151 (MQ=255)
                                                              atCTACAGGTTTAAAAAAAAAAAAAGGCTCCTTTCGGCTTCCCCTGAAAAAATCGCCCTCTCCCCAAGAAGGGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                            <  1:1563489/150‑1 (MQ=255)
                                                              atCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                           >  1:1408764/1‑151 (MQ=255)
                                                              atCTACAGGATTAAAAAAAAAAAAAGCCTCCCTTCGCCTCCCCCTGAAAAAATCCCCCTCTTCCCACGAAGAGGCCCCTAACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                            <  2:1404270/150‑1 (MQ=255)
                                                                  aCAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        >  2:1616360/1‑150 (MQ=255)
                                                                   cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:1457589/1‑151 (MQ=255)
                                                                   cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:972339/1‑151 (MQ=255)
                                                                    aGGATAACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  2:1244166/151‑1 (MQ=255)
                                                                     gaaTTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAAGCCCCCTCTCCCCACAAAGGGGCCCGCTACCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                                    <  1:1333443/149‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGGCTACCTGCGGCTTGCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCAAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:715187/148‑1 (MQ=255)
                                                                          acaaaaaaaaaaaGCCTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                <  2:624845/147‑1 (MQ=255)
                                                                           aaaaaTAAAAAAGGCTCCCTTCGGTTTGCCCTAAAAAAATACCCCTCTCCCCAAGAGGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              <  1:1379809/149‑1 (MQ=255)
                                                                            gAAATAAAAAAGGCTACCTTGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:21647/1‑151 (MQ=255)
                                                                             aaaaaaaaaaGGCCTCCCTTCGGCTGGCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCCCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  1:477172/151‑1 (MQ=255)
                                                                             aaaaaaaaaaGGCTCCCTTCGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            <  1:1338677/151‑1 (MQ=255)
                                                                             aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  2:91200/1‑151 (MQ=255)
                                                                              aaaaaaaaaGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           <  2:1536075/148‑1 (MQ=255)
                                                                              aaTAAAAACGGCTCCCTTCGCCCCCCCCTGAAAAAACGCCCCTCTCCCCCAGAAGGCGCCCACTATCTTTCGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           <  1:1032059/151‑1 (MQ=255)
                                                                                 aaaaaaGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCAACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                        <  2:1477862/151‑1 (MQ=255)
                                                                                    aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     >  1:105509/1‑151 (MQ=255)
                                                                                       ggCTCCCTGCGCCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAGGGGGGCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  <  1:474524/151‑1 (MQ=255)
                                                                                       ggCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  >  1:724695/1‑151 (MQ=255)
                                                                                        gCTCCCTTCGCCTTCCCCTAAAAAAAAACCCCTCTCCCCACGAAGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                 <  2:970071/151‑1 (MQ=255)
                                                                                         gcccccTGCGCCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTgtc                                                                 <  2:1119831/147‑3 (MQ=255)
                                                                                          cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:740363/1‑151 (MQ=255)
                                                                                          cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:1469486/1‑151 (MQ=255)
                                                                                           tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:760869/1‑150 (MQ=255)
                                                                                           tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:1366503/1‑150 (MQ=255)
                                                                                            accctgggCTTTCCCCTAAAAAAAAGCCCCTCTCCCCACAAGGGGGCCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCa                                                               <  1:1287519/145‑1 (MQ=255)
                                                                                            aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                              >  1:124932/1‑151 (MQ=255)
                                                                                            aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                              >  2:368084/1‑151 (MQ=255)
                                                                                             ccTGCGGCTTCCCCTAAAAAAATCGCCCTCTCCCAACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                             <  1:1473455/151‑1 (MQ=255)
                                                                                                 tcggcctccccGAAAAAATCACCCTCCTCCCCACGAAGGGGGCCGTAACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                            <  1:46445/142‑1 (MQ=255)
                                                                                                  ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  1:346856/1‑151 (MQ=255)
                                                                                                   gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                       >  2:488937/1‑151 (MQ=255)
                                                                                                     ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     >  2:1268899/1‑151 (MQ=255)
                                                                                                      cccccTGAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    <  1:474615/149‑1 (MQ=255)
                                                                                                       gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  1:202032/1‑150 (MQ=255)
                                                                                                                aaaaTCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                           <  2:1333998/151‑1 (MQ=255)
                                                                                                                  aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         >  1:181996/1‑151 (MQ=255)
                                                                                                                  aaTACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  2:1252829/151‑1 (MQ=255)
                                                                                                                     ccccccTCTCCCCGCGGAGGGGCCCCCTACCTTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:118452/148‑1 (MQ=255)
                                                                                                                        cccTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAc                                    <  2:542446/151‑1 (MQ=255)
                                                                                                                               ccccACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                                <  2:1057560/148‑1 (MQ=255)
                                                                                                                                cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACAATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                             >  1:1245351/1‑151 (MQ=255)
                                                                                                                                  caggaggagGGCCGCAACCCTAAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           <  1:1531391/145‑1 (MQ=255)
                                                                                                                                   aCGAAGAGGCCCCCCAACCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          <  1:1234369/151‑1 (MQ=255)
                                                                                                                                     gaggggggCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                        <  2:1453963/146‑1 (MQ=255)
                                                                                                                                     gAAGAGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                        <  1:1610697/151‑1 (MQ=255)
                                                                                                                                         gggcccGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                    <  2:806033/147‑1 (MQ=255)
                                                                                                                                         agGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                      >  2:573722/1‑149 (MQ=255)
                                                                                                                                          gggcccccaacccTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                    <  2:1099619/140‑1 (MQ=255)
                                                                                                                                          ggcccGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  2:1050729/148‑1 (MQ=255)
                                                                                                                                              cGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:122645/1‑150 (MQ=255)
                                                                                                                                                 tAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg            >  1:91383/1‑151 (MQ=255)
                                                                                                                                                  aaCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             <  2:1119393/149‑1 (MQ=255)
                                                                                                                                                   acccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc             <  2:679773/147‑1 (MQ=255)
                                                                                                                                                      cTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc        >  2:543771/1‑151 (MQ=255)
                                                                                                                                                       ttAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       >  1:82320/1‑151 (MQ=255)
                                                                                                                                                         aGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCgg     <  2:510083/151‑1 (MQ=255)
                                                                                                                                                            aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:1317794/1‑151 (MQ=255)
                                                                                                                                                            aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:347092/1‑151 (MQ=255)
                                                                                                                                                            aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:1034537/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 < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: