Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,889,695 C→A 6.2% Q135H (CAG→CAT cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,6950CA 6.2% 323.2 / 7.7 96Q135H (CAG→CATcysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (49/41);  new base A (3/3);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGG  >  NC_000913/2889546‑2889841
                                                                                                                                                     |                                                                                                                                                   
gTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTAGAGGTGCAGAGTACGTTACGGTTCATGTCGTTACCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATc                                                                                                                                                     >  2:1232605/1‑150 (MQ=255)
 tACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCTCATCAAGACCGACCGAGTGCAGCATCt                                                                                                                                                    <  1:781383/150‑1 (MQ=255)
     aGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTg                                                                                                                                                >  1:544667/1‑150 (MQ=255)
     aGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTg                                                                                                                                                <  1:309578/150‑1 (MQ=255)
       gCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGAAGCTGCGACTAGTAAGGGTTCGAGGTGCAGAGTACGTTACGTTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCa                                                                                                                                              >  2:742302/1‑150 (MQ=255)
         aaCAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                                                                                                           <  2:17016/150‑1 (MQ=255)
          aCAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCg                                                                                                                                          <  1:586560/150‑1 (MQ=255)
          aCAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCACCAGATGCTGCActc                                                                                                                                           <  2:797114/150‑3 (MQ=255)
           cAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCACCAGATGCTGCActcg                                                                                                                                          <  2:1232534/150‑4 (MQ=255)
             gATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt                                                                                                                                       <  1:264472/150‑1 (MQ=255)
             gATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCACCAGATGCTGCActcggt                                                                                                                                        <  1:1123856/150‑6 (MQ=255)
             gATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCACCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt                                                                                                                                       >  1:1155865/1‑150 (MQ=255)
               tGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTg                                                                                                                                     <  1:740215/150‑1 (MQ=255)
                  tCAGAAATCTTCTTCGCCCACTCGTACGCTTCGGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTGGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                  <  1:603551/150‑1 (MQ=255)
                  tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                  <  2:943280/150‑1 (MQ=255)
                  tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                  >  1:1361664/1‑150 (MQ=255)
                  tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                  <  1:1188861/150‑1 (MQ=255)
                      aaaTCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCt                                                                                                                              >  2:83422/1‑150 (MQ=255)
                       aaTCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                                                                                                                             >  1:1311641/1‑150 (MQ=255)
                          cttcttCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCttttt                                                                                                                          >  2:1468118/1‑150 (MQ=255)
                            tcttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCa                                                                                                                        <  1:1534111/150‑1 (MQ=255)
                              ttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTg                                                                                                                                                <  2:748620/125‑1 (MQ=255)
                              ttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTg                                                                                                                                                >  1:748620/1‑125 (MQ=255)
                                cGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAAt                                                                                                                    >  2:1604410/1‑150 (MQ=255)
                                   ccACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGcc                                                                                                                 <  1:1138706/150‑1 (MQ=255)
                                       tCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg                                                                                                             <  1:631845/150‑1 (MQ=255)
                                        cGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGCGTGCAGCATCTGGTGCCCCGGGTTGACGTTCTTTTTCAGAATGCCGTGGa                                                                                                            >  1:311746/1‑150 (MQ=255)
                                         gTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa                                                                                                           <  2:132834/150‑1 (MQ=255)
                                                tCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAc                                                                                                    >  2:408762/1‑150 (MQ=255)
                                                 ccGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACg                                                                                                   >  2:628487/1‑150 (MQ=255)
                                                      tGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTgg                                                                                              >  2:1145860/1‑150 (MQ=255)
                                                       gCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGTGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGc                                                                                             >  1:453160/1‑150 (MQ=255)
                                                         aGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCgg                                                                                           >  1:886851/1‑150 (MQ=255)
                                                          tcttctaCTCGTAAGGGTTCGATGTGCTTAGTCCGTTACTGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGt                                                                                          <  1:698060/144‑1 (MQ=255)
                                                          gCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGt                                                                                          >  1:1011571/1‑150 (MQ=255)
                                                          gCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGt                                                                                          <  1:1392137/150‑1 (MQ=255)
                                                             gCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTgg                                                                                       >  2:1236560/1‑150 (MQ=255)
                                                               gACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTc                                                                                     <  2:212283/150‑1 (MQ=255)
                                                               gACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTc                                                                                     <  1:1432308/150‑1 (MQ=255)
                                                                  tCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAgg                                                                                  <  1:52341/150‑1 (MQ=255)
                                                                   cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                                                                 >  1:380926/1‑150 (MQ=255)
                                                                   cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                                                                 >  1:156882/1‑150 (MQ=255)
                                                                    gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCg                                                                                >  2:493978/1‑150 (MQ=255)
                                                                     tAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGa                                                                               >  1:1322033/1‑150 (MQ=255)
                                                                       aGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAt                                                                             <  2:1311641/150‑1 (MQ=255)
                                                                             cGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGcc                                                                       >  1:1181846/1‑150 (MQ=255)
                                                                               aGGTTCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATAAAGACCGTCCGAGTGCAGCATCTGGTGCACCGTTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACTTCTGGCGGTTTGTCAGGCTAATGCTGCCAt                                                                     <  2:606064/150‑1 (MQ=255)
                                                                                  tGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGa                                                                  >  1:323618/1‑150 (MQ=255)
                                                                                         taggTTTCGGTTCATGTCGTTATCTGTCGCCAGCGCATGAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTGAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGtt                                                           <  2:1102268/147‑1 (MQ=255)
                                                                                          aCGTTACGGTTCATGTCGTTAGCTGTCTCAAGCGCATCAAGTCCGACCGAGTGCAGCATCTTGTTAACCGGTTTGAAGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGGGGTTGGTCAGGCGAATGCTGCCATATATTGTGttt                                                          <  2:1322033/150‑1 (MQ=255)
                                                                                          aCGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                                                          >  1:1583305/1‑150 (MQ=255)
                                                                                           cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCGATAGATGGTGtttt                                                         >  1:608375/1‑150 (MQ=255)
                                                                                           cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGtttt                                                         <  2:156882/150‑1 (MQ=255)
                                                                                             ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCa                                                       >  1:1004464/1‑150 (MQ=255)
                                                                                             ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCa                                                       <  1:133403/150‑1 (MQ=255)
                                                                                               aCGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc                                                     <  2:544667/150‑1 (MQ=255)
                                                                                                      aTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAt                                              <  1:1236560/150‑1 (MQ=255)
                                                                                                          cGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGt                                          <  2:380926/150‑1 (MQ=255)
                                                                                                          cGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGt                                          >  2:1518853/1‑150 (MQ=255)
                                                                                                           gTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTc                                         <  1:612028/150‑1 (MQ=255)
                                                                                                              aGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGAt                                      >  1:1546358/1‑150 (MQ=255)
                                                                                                                ttGTCGTCAGCGCATCAAGACCGACCGAGTGGAGCATCTGGTGCAGCGGTTTGACGTTCTTTTTCAGAATGGCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGGCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCg                                    <  2:323618/149‑1 (MQ=255)
                                                                                                                 tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCATAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGc                                   >  2:400691/1‑150 (MQ=255)
                                                                                                                 tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGc                                   >  1:1442839/1‑150 (MQ=255)
                                                                                                                 tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGc                                   >  2:1367698/1‑150 (MQ=255)
                                                                                                                  ggcGCCAGCGCATCAAGCCCGACCGAGGGCAGCATCTGGTGCACCGGTTTGACGTGCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGcc                                  <  1:984399/148‑1 (MQ=255)
                                                                                                                   tCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCt                                 <  2:93900/150‑1 (MQ=255)
                                                                                                                    cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTg                                >  1:1103765/1‑150 (MQ=255)
                                                                                                                    cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTg                                >  2:352929/1‑150 (MQ=255)
                                                                                                                      ccAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGcc                              >  1:1428163/1‑150 (MQ=255)
                                                                                                                        aGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGACAAATTTGTCGATCGCCTGCCAc                            >  2:1134354/1‑150 (MQ=255)
                                                                                                                         gcgcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACt                           <  1:573915/150‑1 (MQ=255)
                                                                                                                          cgcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTg                          >  1:431015/1‑150 (MQ=255)
                                                                                                                           gcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGt                         <  2:442006/150‑1 (MQ=255)
                                                                                                                           gcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGt                         >  2:802744/1‑150 (MQ=255)
                                                                                                                              tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCa                                                                      >  2:868074/1‑102 (MQ=255)
                                                                                                                              tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCa                                                                      <  1:868074/102‑1 (MQ=255)
                                                                                                                              tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTa                      >  1:731131/1‑150 (MQ=255)
                                                                                                                              tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTa                      >  2:361840/1‑150 (MQ=255)
                                                                                                                                aaGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGt                    >  2:1505113/1‑150 (MQ=255)
                                                                                                                                aaGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGt                    <  1:1145860/150‑1 (MQ=255)
                                                                                                                                  gaccgaccGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTgg                  <  2:1181846/150‑1 (MQ=255)
                                                                                                                                    ccgaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTc                                         <  1:834220/125‑1 (MQ=255)
                                                                                                                                    ccgaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTc                                         >  2:834220/1‑125 (MQ=255)
                                                                                                                                     cgaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTaa               >  2:1446668/1‑150 (MQ=255)
                                                                                                                                       accgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTAAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATc             >  2:170659/1‑150 (MQ=255)
                                                                                                                                       aaccggTGCACCAGATGCTGCACTCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAAt              >  1:1196005/6‑150 (MQ=255)
                                                                                                                                       aaccggTGCACCAGATGCTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATc             >  2:1223728/6‑150 (MQ=255)
                                                                                                                                         ccggTGCACCAGATGCTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcac           >  1:691052/4‑150 (MQ=255)
                                                                                                                                          gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCTAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacc          <  1:83422/150‑1 (MQ=255)
                                                                                                                                          gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacc          >  1:315656/1‑150 (MQ=255)
                                                                                                                                          gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacc          >  1:557615/1‑150 (MQ=255)
                                                                                                                                             tGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacccca       >  2:108772/1‑150 (MQ=255)
                                                                                                                                               cagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcaccccacc     <  2:315656/150‑1 (MQ=255)
                                                                                                                                                agcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcaccccaccc    >  1:444475/1‑150 (MQ=255)
                                                                                                                                                  caTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCgg  >  1:1162311/1‑150 (MQ=255)
                                                                                                                                                  caTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCgg  <  1:408762/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGG  >  NC_000913/2889546‑2889841

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: