Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,877,562 C→A 7.7% R328S (CGT→AGT) ‡ 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,5620CA 7.7% 214.0 / 6.3 65R328S (CGT→AGT) ‡iapalkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base A (2/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTT  >  NC_000913/2877417‑2877706
                                                                                                                                                 |                                                                                                                                                
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)
                                                                                                                                                 |                                                                                                                                                
TAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTT  >  NC_000913/2877417‑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: ‑