Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,658,422 C→G 7.0% A172P (GCG→CCG)  hscA ← iron‑sulfur cluster biosynthesis chaperone HscA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,658,4220CG 7.0% 386.9 / 15.6 114A172P (GCG→CCG) hscAiron‑sulfur cluster biosynthesis chaperone HscA
Reads supporting (aligned to +/- strand):  ref base C (52/54);  new base G (3/5);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 116 total aligned reads displayed.

GCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCC  >  NC_000913/2658279‑2658571
                                                                                                                                               |                                                                                                                                                     
gcgcAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                                                                                                                                                 <  1:919195/150‑1 (MQ=255)
   cAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCttt                                                                                                                                              >  1:819049/1‑150 (MQ=255)
     gAATGGAAATATCAAACGTCCCGCCACTGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtggcg                                                                                                                                            <  1:30146/150‑5 (MQ=255)
     gAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtggcg                                                                                                                                            <  1:475926/150‑5 (MQ=255)
     gAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTgg                                                                                                                                            <  1:1351828/150‑1 (MQ=255)
      aaTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                                                                                                                                           >  2:914396/1‑150 (MQ=255)
       aTGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtggcggg                                                                                                                                          <  1:1482279/150‑7 (MQ=255)
        tGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGGCGCGTCtggcgggc                                                                                                                                         <  1:1075228/150‑8 (MQ=255)
        tGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGc                                                                                                                                         <  2:822210/150‑1 (MQ=255)
          gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       <  2:572451/150‑1 (MQ=255)
          gAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                                                                                                                                       >  1:842196/1‑150 (MQ=255)
             atatCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                                                                                                                                    <  2:1570173/150‑1 (MQ=255)
               atCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACg                                                                                                                                  >  1:137750/1‑150 (MQ=255)
               atCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACg                                                                                                                                  >  1:1023169/1‑150 (MQ=255)
                 cAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCt                                                                                                                                <  2:769550/150‑1 (MQ=255)
                  aaaCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg                                                                                                                               >  1:1080288/1‑150 (MQ=255)
                  aaaCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg                                                                                                                               <  2:1448611/150‑1 (MQ=255)
                   aaCGTCCCGCCACCGAGGTCATAAACGGCGATCGCGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTTGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg                                                                                                                              <  2:925120/150‑1 (MQ=255)
                     cGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc                                                                                                                            >  1:1497061/1‑150 (MQ=255)
                        cccGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc                                                                                                                         >  2:1302593/1‑150 (MQ=255)
                           gCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTc                                                                                                                      >  1:970408/1‑150 (MQ=255)
                            ccACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCa                                                                                                                     <  1:1584458/150‑1 (MQ=255)
                              aCCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCaaa                                                                                                                   <  1:914396/150‑1 (MQ=255)
                               ccGAGGTCATAAACGGCGATCACGCCTTCCCGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAg                                                                                                                  <  1:881788/150‑1 (MQ=255)
                                   ggTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACg                                                                                                              <  1:456465/150‑1 (MQ=255)
                                    gTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGc                                                                                                             >  2:1015801/1‑150 (MQ=255)
                                     tCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGcc                                                                                                            <  1:100418/150‑1 (MQ=255)
                                       aTAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCgg                                                                                                          <  2:1161525/150‑1 (MQ=255)
                                         aaaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGaa                                                                                                        >  1:903297/1‑150 (MQ=255)
                                         aaaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGaa                                                                                                        >  2:807878/1‑150 (MQ=255)
                                          aaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAAc                                                                                                       >  1:831956/1‑150 (MQ=255)
                                          aaCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAAc                                                                                                       <  1:477711/150‑1 (MQ=255)
                                                   cACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAAc                                                                                              >  2:250373/1‑150 (MQ=255)
                                                     cGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTa                                                                                            >  2:1340619/1‑150 (MQ=255)
                                                       ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                          >  2:1608399/1‑150 (MQ=255)
                                                       ccTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTaca                                                                                          <  2:903297/150‑1 (MQ=255)
                                                        cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTacac                                                                                         <  2:1042415/150‑1 (MQ=255)
                                                        cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTacac                                                                                         <  1:1183469/150‑1 (MQ=255)
                                                        cTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTacac                                                                                         >  1:1011041/1‑150 (MQ=255)
                                                            cTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATc                                                                                     >  1:391339/1‑150 (MQ=255)
                                                             tGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATcc                                                                                    >  1:744282/1‑150 (MQ=255)
                                                             tGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATcc                                                                                    <  1:1599440/150‑1 (MQ=255)
                                                                ccGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGc                                                                                 >  1:319108/1‑150 (MQ=255)
                                                                 cGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCt                                                                                <  2:1023169/150‑1 (MQ=255)
                                                                   gAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCg                                                                              <  2:842196/150‑1 (MQ=255)
                                                                    aaTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGc                                                                             <  1:1207232/150‑1 (MQ=255)
                                                                      tCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCt                                                                           <  2:536552/150‑1 (MQ=255)
                                                                       ccAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTg                                                                          >  2:1119131/1‑150 (MQ=255)
                                                                       ccAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTg                                                                          >  1:1311341/1‑150 (MQ=255)
                                                                       ccAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTg                                                                          >  2:454094/1‑150 (MQ=255)
                                                                        cAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGc                                                                         >  2:902375/1‑150 (MQ=255)
                                                                            ccGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAgg                                                                     <  2:693849/150‑1 (MQ=255)
                                                                             cGTAGGCGATAGCCGTAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAggg                                                                    <  1:1110219/150‑1 (MQ=255)
                                                                              gTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGc                                                                   <  2:137750/150‑1 (MQ=255)
                                                                                 ggCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTc                                                                >  1:746562/1‑150 (MQ=255)
                                                                                 ggCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTc                                                                <  2:1650100/150‑1 (MQ=255)
                                                                                  gCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCa                                                               <  1:1244365/150‑1 (MQ=255)
                                                                                   cGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAg                                                              <  1:1006113/150‑1 (MQ=255)
                                                                                        gCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGCCGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGcc                                                         >  1:911986/1‑150 (MQ=255)
                                                                                         ccGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGccc                                                        >  1:1257903/1‑150 (MQ=255)
                                                                                          cGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCcg                                                       >  1:35283/1‑150 (MQ=255)
                                                                                           gCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCcgc                                                      >  1:1453152/1‑150 (MQ=255)
                                                                                             aGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCgcgc                                                    <  1:1016577/150‑1 (MQ=255)
                                                                                              gCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGcgcc                                                   >  1:33322/1‑150 (MQ=255)
                                                                                               cggtcggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGcgccg                                                  >  2:1644853/1‑150 (MQ=255)
                                                                                                 gtcggtTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCgccgcc                                                <  2:1438897/150‑1 (MQ=255)
                                                                                                        cGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                         >  2:164142/1‑150 (MQ=255)
                                                                                                        cGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCtt                                         >  2:18982/1‑150 (MQ=255)
                                                                                                         gTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCttt                                        >  1:425817/1‑150 (MQ=255)
                                                                                                          ttAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTg                                       <  1:1015801/150‑1 (MQ=255)
                                                                                                           taagtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGa                                      <  1:1196919/150‑1 (MQ=255)
                                                                                                             agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAgg                                    <  2:695769/150‑1 (MQ=255)
                                                                                                             agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAgg                                    >  1:1160381/1‑150 (MQ=255)
                                                                                                             agtaagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAgg                                    <  1:1183368/150‑1 (MQ=255)
                                                                                                                aagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATg                                 >  1:388537/1‑150 (MQ=255)
                                                                                                                aagCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCATCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATg                                 <  1:273242/150‑1 (MQ=255)
                                                                                                                 agCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGt                                >  2:1128116/1‑150 (MQ=255)
                                                                                                                 agCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGt                                >  1:108105/1‑150 (MQ=255)
                                                                                                                   cgcAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTcc                              <  2:819049/150‑1 (MQ=255)
                                                                                                                    gcAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCg                             <  2:1168602/150‑1 (MQ=255)
                                                                                                                     cAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGc                            <  1:902375/150‑1 (MQ=255)
                                                                                                                       ggACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAg                          >  2:1518036/1‑150 (MQ=255)
                                                                                                                          cGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGaaa                       >  1:772286/1‑150 (MQ=255)
                                                                                                                            tGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAAcg                     <  1:1128116/150‑1 (MQ=255)
                                                                                                                                 gCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAc                >  2:943621/1‑150 (MQ=255)
                                                                                                                                  cccGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc               >  1:1070412/1‑150 (MQ=255)
                                                                                                                                  caccaaaGACGCGGTGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAAATACACCATCCAGCTCGCCTGCCAGGGCTGCAGTTGCCCGCGCCGCCAGTGCGTTGAGGATGTGCGCAGAAACGCGCAcc               >  2:1109851/7‑150 (MQ=255)
                                                                                                                                  caccaaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCAcc               >  1:335275/7‑150 (MQ=255)
                                                                                                                                    ccaaaGACGCGGCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCgg             >  1:1091282/5‑150 (MQ=255)
                                                                                                                                     gacAGACGCTCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCggg            >  2:445557/3‑150 (MQ=255)
                                                                                                                                      ccAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGt           <  2:388537/150‑1 (MQ=255)
                                                                                                                                        aGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTc         <  1:1361354/150‑1 (MQ=255)
                                                                                                                                         gACGCGCCTCGTCTTTGGTGCCCTGACGATGGGAATCGTCAAAGTACGCCTGAACGGTGATAACTACACCATCCAGCTCGACTGACAGGGATTCAGTTTCCCGCGCCTCCAGTTCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTca        <  2:1497061/150‑1 (MQ=255)
                                                                                                                                         gACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTca        <  2:108105/150‑1 (MQ=255)
                                                                                                                                           cgcgCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGTTAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTcagc      <  2:319108/150‑1 (MQ=255)
                                                                                                                                           cgcgCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACgcgc                  <  1:1135338/138‑1 (MQ=255)
                                                                                                                                           cgcgCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACgcgc                  >  2:1135338/1‑138 (MQ=255)
                                                                                                                                           cgcgCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTcagc      >  2:425892/1‑150 (MQ=255)
                                                                                                                                              gccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCagcagc   <  1:1119131/150‑1 (MQ=255)
                                                                                                                                               ccgcGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGcc  >  1:1467948/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTCAAAGTACGCCGGAACGGTGATAACTACACCATCCAGCTCGCCTGCCAGGGCTTCAGTTGCCCGCGCCGCCAGTGCTTTGAGGATGTCCGCAGAAACGCGCACCGGGTTCAGCAGCC  >  NC_000913/2658279‑2658571

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: