Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,889,691 G→C 6.4% L137V (CTG→GTG)  cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,6910GC 6.4% 314.1 / 9.6 93L137V (CTG→GTG) cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base C (3/3);  total (51/42)
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.36e-01

GCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCC  >  NC_000913/2889542‑2889839
                                                                                                                                                     |                                                                                                                                                     
gcgGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcag                                                                                                                                                       <  2:820738/150‑1 (MQ=255)
    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                                                                                                                                              >  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: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: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: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                                 >  2:1367698/1‑150 (MQ=255)
                                                                                                                     tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGc                                 >  1:1442839/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                    >  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)
                                                                                                                                                     |                                                                                                                                                     
GCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCC  >  NC_000913/2889542‑2889839

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: