Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,877,559 T→G 7.6% C327G (TGC→GGC)  iap → alkaline phosphatase isozyme conversion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,877,5590TG 7.6% 200.7 / 8.3 66C327G (TGC→GGC) iapalkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base G (2/3);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

AGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTT  >  NC_000913/2877411‑2877706
                                                                                                                                                    |                                                                                                                                                   
aGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTg                                                                                                                                                    >  1:468581/1‑150 (MQ=255)
      tAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGa                                                                                                                                              >  2:873907/1‑150 (MQ=255)
        aCTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACg                                                                                                                                            <  2:520269/150‑1 (MQ=255)
           ggAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttc                                                                                                                                         <  2:371309/150‑2 (MQ=255)
           ggAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttc                                                                                                                                         <  2:366024/150‑2 (MQ=255)
           ggAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttc                                                                                                                                         >  1:371309/1‑149 (MQ=255)
           ggAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttc                                                                                                                                         >  1:366024/1‑149 (MQ=255)
              aTCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttctat                                                                                                                                      <  2:547129/150‑5 (MQ=255)
                           ggATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTAcc                                                                                                                         >  1:1416018/1‑150 (MQ=255)
                               ggTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTg                                                                                                                     >  2:411865/1‑150 (MQ=255)
                                 ttATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGt                                                                                                                   <  2:565356/150‑1 (MQ=255)
                                  tATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTg                                                                                                                  <  2:71404/150‑1 (MQ=255)
                                   aTCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGa                                                                                                                 >  1:1400318/1‑150 (MQ=255)
                                      aGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAgg                                                                                                              >  1:872278/1‑150 (MQ=255)
                                        cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAg                                                                                                            <  2:1069348/150‑1 (MQ=255)
                                          aCGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGtt                                                                                                          >  1:157810/1‑150 (MQ=255)
                                           cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTg                                                                                                         >  2:322828/1‑150 (MQ=255)
                                             cgcAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAGTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGc                                                                                                       >  2:1054910/1‑150 (MQ=255)
                                              gcAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCg                                                                                                      >  2:53312/1‑150 (MQ=255)
                                                  aaaCACCTGCCTTCCCGGCGGGAAATAGTTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAgg                                                                                                  <  1:436453/150‑1 (MQ=255)
                                                  aaaCACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAgg                                                                                                  <  2:369762/150‑1 (MQ=255)
                                                  aaaCACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAgg                                                                                                  <  1:1008883/150‑1 (MQ=255)
                                                     cacCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGt                                                                                               >  2:500241/1‑150 (MQ=255)
                                                       ccTGCCTTCCCGGCGTGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCt                                                                                             <  1:723825/150‑1 (MQ=255)
                                                        cTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCtt                                                                                            <  2:348384/150‑1 (MQ=255)
                                                            cTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATg                                                                                        <  2:754019/150‑1 (MQ=255)
                                                              tCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATggg                                                                                      <  1:986488/150‑1 (MQ=255)
                                                               cccGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGt                                                                                     >  1:620136/1‑150 (MQ=255)
                                                               cccGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGt                                                                                     >  1:1388107/1‑150 (MQ=255)
                                                                      ggAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGaaaa                                                                              >  2:317526/1‑150 (MQ=255)
                                                                       gAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAAt                                                                             >  2:701390/1‑150 (MQ=255)
                                                                         aaTAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATgg                                                                           >  2:44274/1‑150 (MQ=255)
                                                                          aTAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATggg                                                                          <  2:1416018/150‑1 (MQ=255)
                                                                          aTAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATggg                                                                          <  1:701390/150‑1 (MQ=255)
                                                                            aGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAg                                                                        >  2:550047/1‑150 (MQ=255)
                                                                             gCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGc                                                                       <  2:157810/150‑1 (MQ=255)
                                                                             gCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGc                                                                       >  2:37136/1‑150 (MQ=255)
                                                                              cTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCt                                                                      >  2:470123/1‑150 (MQ=255)
                                                                              cTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCt                                                                      >  2:682721/1‑150 (MQ=255)
                                                                                 gCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCTAATGCGTCTTGATGGGTTTGAAAATGGGAGCTggg                                                                   <  2:1077937/150‑1 (MQ=255)
                                                                                 gCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTggg                                                                   >  2:827994/1‑150 (MQ=255)
                                                                                          aaGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTAcc                                                          <  1:411865/150‑1 (MQ=255)
                                                                                              cTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCag                                                      <  1:53312/150‑1 (MQ=255)
                                                                                                 gATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAgg                                                   >  1:1220770/1‑150 (MQ=255)
                                                                                                     aTCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCggg                                               >  1:1172351/1‑150 (MQ=255)
                                                                                                       cACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCggggg                                             >  2:545353/1‑150 (MQ=255)
                                                                                                         ccAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGaa                                           >  2:689866/1‑150 (MQ=255)
                                                                                                           aaCATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACt                                         <  2:1220770/150‑1 (MQ=255)
                                                                                                             cATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTcc                                       <  1:470123/150‑1 (MQ=255)
                                                                                                              atatTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCa                                      <  1:37136/150‑1 (MQ=255)
                                                                                                               tatTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCaa                                     >  2:460716/1‑150 (MQ=255)
                                                                                                               tatTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTAACGCAGAGGCGGGGGAACTCCaa                                     >  1:939250/1‑150 (MQ=255)
                                                                                                                  tGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTg                                  >  1:819416/1‑150 (MQ=255)
                                                                                                                   gATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGa                                 >  2:290052/1‑150 (MQ=255)
                                                                                                                       aGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATc                             >  2:120470/1‑150 (MQ=255)
                                                                                                                        ggCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATcc                            >  2:571009/1‑150 (MQ=255)
                                                                                                                               ccTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCg                     <  1:827994/150‑1 (MQ=255)
                                                                                                                                cTGGAAGAATAGAACGTCGCTGCCGTGACATTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGc                    >  1:715544/1‑150 (MQ=255)
                                                                                                                                      gaaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAg              <  1:120470/150‑1 (MQ=255)
                                                                                                                                        ttAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTg            <  1:1054910/149‑1 (MQ=255)
                                                                                                                                        aTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTg            <  2:1388107/150‑1 (MQ=255)
                                                                                                                                        aTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTg            >  2:886159/1‑150 (MQ=255)
                                                                                                                                         tAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTgc           >  2:1115390/1‑150 (MQ=255)
                                                                                                                                                 cGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGGCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGtt   <  2:468581/150‑1 (MQ=255)
                                                                                                                                                  gCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGttt  >  1:546860/1‑150 (MQ=255)
                                                                                                                                                  gCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGttt  <  1:463579/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTT  >  NC_000913/2877411‑2877706

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