Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,877,557 G→A 6.2% R326H (CGC→CAC)  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,5570GA 6.2% 353.3 / 6.3 111R326H (CGC→CAC) iapalkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base G (53/51);  new base A (5/2);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGC  >  NC_000913/2877409‑2877699
                                                                                                                                                    |                                                                                                                                              
gAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGc                                                                                                                                               <  1:1843333/150‑1 (MQ=255)
gAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGc                                                                                                                                               <  2:1547776/150‑1 (MQ=255)
gAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGc                                                                                                                                               <  1:1027776/150‑1 (MQ=255)
 aaGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCt                                                                                                                                              >  2:286955/1‑150 (MQ=255)
  aGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTg                                                                                                                                             >  1:1291963/1‑150 (MQ=255)
       cTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                                                                                                                                        >  2:1463478/1‑150 (MQ=255)
       cTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                                                                                                                                        >  1:1126731/1‑150 (MQ=255)
          aCTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACg                                                                                                                                     >  1:1293762/1‑150 (MQ=255)
          aCTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACg                                                                                                                                     >  1:1735458/1‑150 (MQ=255)
           cTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                                                                                                                                    >  1:1834155/1‑150 (MQ=255)
            tGGAATCTTTGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTTCCGTGACGtt                                                                                                                                   >  2:543530/1‑150 (MQ=255)
               aaTCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttcta                                                                                                                                <  2:794310/150‑4 (MQ=255)
                aTCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGc                                                                                                                               <  1:299034/150‑1 (MQ=255)
                  cTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCACGGCAGCGACGttctattc                                                                                                                             <  1:767166/150‑7 (MQ=255)
                   ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGa                                                                                                                            >  1:1008999/1‑150 (MQ=255)
                   ttGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGa                                                                                                                            >  1:1334950/1‑150 (MQ=255)
                      ggAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATaa                                                                                                                         >  1:220893/1‑150 (MQ=255)
                        aaTAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGTGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATg                                                                                                                       <  2:1069119/150‑1 (MQ=255)
                        aaTAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGACAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATg                                                                                                                       <  2:35751/150‑1 (MQ=255)
                            aGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTAc                                                                                                                   >  2:1456710/1‑150 (MQ=255)
                             ggATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTAcc                                                                                                                  <  2:1293637/150‑1 (MQ=255)
                             ggATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTAcc                                                                                                                  <  2:1334950/150‑1 (MQ=255)
                             ggATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTAcc                                                                                                                  <  2:220893/150‑1 (MQ=255)
                                 ggTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTg                                                                                                              >  1:546690/1‑150 (MQ=255)
                                   ttATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGt                                                                                                            <  2:1225998/150‑1 (MQ=255)
                                     aTCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGa                                                                                                          <  2:1119881/150‑1 (MQ=255)
                                          cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAg                                                                                                     >  1:1849383/1‑150 (MQ=255)
                                             cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGc                                                                                                                      >  2:137918/1‑130 (MQ=255)
                                             cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGc                                                                                                                      <  1:137918/130‑1 (MQ=255)
                                             cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTg                                                                                                  >  2:749810/1‑150 (MQ=255)
                                             cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTg                                                                                                  >  2:1659241/1‑150 (MQ=255)
                                              gcgcAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGTATAATCACCAACATATTGATAAGGCTCTTCCTGCAATAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTgg                                                                                                 <  2:66430/150‑1 (MQ=255)
                                               cgcAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGc                                                                                                >  1:1040015/1‑150 (MQ=255)
                                                 cAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTg                                                                                                           >  1:1866121/1‑137 (MQ=255)
                                                 cAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTg                                                                                                           <  2:1866121/137‑1 (MQ=255)
                                                 cAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGa                                                                                              >  2:170669/1‑150 (MQ=255)
                                                    aaaCACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAgg                                                                                           <  1:1659241/150‑1 (MQ=255)
                                                      acacCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCg                                                                                         <  2:1522628/150‑1 (MQ=255)
                                                        acCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTc                                                                                       <  2:311673/150‑1 (MQ=255)
                                                          cTGCCTTCCCGGCGGGAAATAGCTGTCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCtt                                                                                     <  2:1691990/150‑1 (MQ=255)
                                                          cTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCtt                                                                                     <  1:119887/150‑1 (MQ=255)
                                                           tGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTg                                                                                    >  2:84829/1‑150 (MQ=255)
                                                            gCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGa                                                                                   <  2:1735458/150‑1 (MQ=255)
                                                            gCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGa                                                                                   >  1:1822394/1‑150 (MQ=255)
                                                            accTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCg                                                                                                                              <  2:942400/106‑1 (MQ=255)
                                                            accTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCg                                                                                                                              >  1:942400/2‑107 (MQ=255)
                                                             ccTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGAt                                                                                  >  2:4628/1‑150 (MQ=255)
                                                              cTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATg                                                                                 >  1:1751506/1‑150 (MQ=255)
                                                              cTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATg                                                                                 <  1:1501361/150‑1 (MQ=255)
                                                                        ggAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGaaaa                                                                       >  2:1813576/1‑150 (MQ=255)
                                                                              aGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAg                                                                 <  2:1823251/150‑1 (MQ=255)
                                                                               gCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGc                                                                >  2:1355377/1‑150 (MQ=255)
                                                                                cTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCt                                                               >  2:1677498/1‑150 (MQ=255)
                                                                                   gCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTggg                                                            <  2:546690/150‑1 (MQ=255)
                                                                                    cATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGa                                                           >  1:1387350/1‑150 (MQ=255)
                                                                                      tGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGt                                                         <  1:1184947/150‑1 (MQ=255)
                                                                                       gACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGtt                                                        >  2:1058663/1‑150 (MQ=255)
                                                                                        aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTc                                                       >  1:1310278/1‑150 (MQ=255)
                                                                                        aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTc                                                       <  2:1040015/150‑1 (MQ=255)
                                                                                         cGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCt                                                      >  2:653269/1‑150 (MQ=255)
                                                                                          gTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTa                                                     <  1:36833/150‑1 (MQ=255)
                                                                                          gTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTa                                                     >  2:1001593/1‑150 (MQ=255)
                                                                                              gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGc                                                 <  2:367068/150‑1 (MQ=255)
                                                                                                  ggATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCagag                                             >  2:559356/1‑150 (MQ=255)
                                                                                                  ggATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCagag                                             >  1:1024080/1‑150 (MQ=255)
                                                                                                        tCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCgggg                                       <  1:1001593/150‑1 (MQ=255)
                                                                                                        tCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCgggg                                       <  2:1588841/150‑1 (MQ=255)
                                                                                                         cACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCggggg                                      <  1:286955/150‑1 (MQ=255)
                                                                                                         cACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCggggg                                      >  1:921750/1‑150 (MQ=255)
                                                                                                          aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGa                                     >  1:340072/1‑150 (MQ=255)
                                                                                                          aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGa                                     >  2:1837294/1‑150 (MQ=255)
                                                                                                          aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGa                                     <  2:1953721/150‑1 (MQ=255)
                                                                                                             aaCATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACt                                  <  2:1849383/150‑1 (MQ=255)
                                                                                                              aCATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTTCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTc                                 <  1:543530/150‑1 (MQ=255)
                                                                                                                  atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGTATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGATCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAg                             <  2:206402/150‑1 (MQ=255)
                                                                                                                  atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAg                             <  2:1310278/150‑1 (MQ=255)
                                                                                                                  atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAg                             <  1:170669/150‑1 (MQ=255)
                                                                                                                  atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAg                             >  1:866829/1‑150 (MQ=255)
                                                                                                                   ttGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGt                            >  1:1704668/1‑150 (MQ=255)
                                                                                                                   ttGATAAGGATCTTACTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGt                            >  2:602756/1‑150 (MQ=255)
                                                                                                                    tGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTg                           <  1:624322/150‑1 (MQ=255)
                                                                                                                    tGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTg                           <  2:340072/150‑1 (MQ=255)
                                                                                                                     gATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGa                          >  1:716660/1‑150 (MQ=255)
                                                                                                                       tAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGata                        <  2:866829/150‑1 (MQ=255)
                                                                                                                         aGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATc                      >  2:1168053/1‑150 (MQ=255)
                                                                                                                         aGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATc                      >  1:342296/1‑150 (MQ=255)
                                                                                                                          ggCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATcc                     <  2:1217547/150‑1 (MQ=255)
                                                                                                                            ctctTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCcat                   <  1:749810/150‑1 (MQ=255)
                                                                                                                            ctctTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCcat                   >  1:648386/1‑150 (MQ=255)
                                                                                                                                  cTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAGGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGc             >  2:504814/1‑150 (MQ=255)
                                                                                                                                   tGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCa            >  2:1694779/1‑150 (MQ=255)
                                                                                                                                      aataaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTTCCTCTGGTGAAGGAGTTGGCGTAGGCGTCTTGATGGGTTTGAAAATGGGATCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATcc         <  2:1694531/147‑1 (MQ=255)
                                                                                                                                        gaaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAg       <  2:730922/150‑1 (MQ=255)
                                                                                                                                         cgcataaCGTCACGGCAGCGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGt      >  1:1089662/6‑150 (MQ=255)
                                                                                                                                         cgcataaCGTCACGGCAGCGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGt      >  1:1554115/6‑150 (MQ=255)
                                                                                                                                         cgcataaCGTCACCGCAGCGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGt      >  2:921224/6‑150 (MQ=255)
                                                                                                                                          gcataaCGTCACGGCAGCGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTg     >  1:845176/5‑150 (MQ=255)
                                                                                                                                          gcataaCGTCACGGCAGCGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTg     >  2:1030323/5‑150 (MQ=255)
                                                                                                                                             gAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTgcgc  >  1:1429633/1‑150 (MQ=255)
                                                                                                                                             gAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTgcgc  <  2:1704668/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GAAGCGACTAACTGGAATCTTGGGAATAAGGATGGTTATCAGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGC  >  NC_000913/2877409‑2877699

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