Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,889,698 G→C 6.3% H134Q (CAC→CAG cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,6980GC 6.3% 318.2 / 9.5 96H134Q (CAC→CAGcysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (49/41);  new base C (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 = 8.70e-01

AGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGA  >  NC_000913/2889551‑2889845
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                      <  1:1188861/150‑1 (MQ=255)
             tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                      >  1:1361664/1‑150 (MQ=255)
             tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                                                                      <  2:943280/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                                                                                                                                                    >  1:748620/1‑125 (MQ=255)
                         ttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTg                                                                                                                                                    <  2:748620/125‑1 (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:1392137/150‑1 (MQ=255)
                                                     gCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGt                                                                                              >  1:1011571/1‑150 (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:156882/1‑150 (MQ=255)
                                                              cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                                                                     >  1:380926/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:133403/150‑1 (MQ=255)
                                                                                        ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCa                                                           >  1:1004464/1‑150 (MQ=255)
                                                                                          aCGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc                                                         <  2:544667/150‑1 (MQ=255)
                                                                                                 aTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAt                                                  <  1:1236560/150‑1 (MQ=255)
                                                                                                     cGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGt                                              >  2:1518853/1‑150 (MQ=255)
                                                                                                     cGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGt                                              <  2:380926/150‑1 (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                                                                          <  1:868074/102‑1 (MQ=255)
                                                                                                                         tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCa                                                                          >  2:868074/1‑102 (MQ=255)
                                                                                                                         tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTa                          >  2:361840/1‑150 (MQ=255)
                                                                                                                         tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTa                          >  1:731131/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:557615/1‑150 (MQ=255)
                                                                                                                                     gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacc              >  1:315656/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)
                                                                                                                                                 tGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAga  <  1:50082/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGA  >  NC_000913/2889551‑2889845

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: