Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,257 G→C 6.6% L640L (CTG→CTC fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2570GC 6.6% 313.6 / 8.2 102L640L (CTG→CTCfhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (44/51);  new base C (5/2);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 104 total aligned reads displayed.

ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGG  >  NC_000913/2856109‑2856405
                                                                                                                                                    |                                                                                                                                                    
aCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTgg                                                                                                                                                     <  2:772356/150‑1 (MQ=255)
   tGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAgg                                                                                                                                                  <  1:1486762/150‑1 (MQ=255)
   tGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAgg                                                                                                                                                  <  1:811137/150‑1 (MQ=255)
      aGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCTCCTGCCGCAACGGTTGTCGCCCTGGAGGGCg                                                                                                                                               >  2:58618/1‑150 (MQ=255)
       gTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGa                                                                                                                                              >  1:434910/1‑150 (MQ=255)
         ggCCGGGTAACGTACGCGAACTGGAAAACGGCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAg                                                                                                                                            <  1:109252/150‑1 (MQ=255)
             gggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGa                                                                                                                                        <  2:1654443/150‑1 (MQ=255)
              ggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                       >  2:353114/1‑150 (MQ=255)
              ggTAACGTACGAGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                       <  2:1385541/150‑1 (MQ=255)
                  aCGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATc                                                                                                                                   <  2:905135/150‑1 (MQ=255)
                  aCGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATc                                                                                                                                   >  1:715341/1‑150 (MQ=255)
                     tACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGt                                                                                                                                >  1:977803/1‑150 (MQ=255)
                         cgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAt                                                                                                                            >  2:1667585/1‑150 (MQ=255)
                         cgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAt                                                                                                                            >  1:1425387/1‑150 (MQ=255)
                               ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                      >  1:1216440/1‑150 (MQ=255)
                               ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                      >  2:697210/1‑150 (MQ=255)
                               ggAAAACGTCATAGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGAAGAATATCAGTTGATTGTgcg                                                                                                                      <  2:165025/150‑1 (MQ=255)
                                gAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcgc                                                                                                                     >  2:48903/1‑150 (MQ=255)
                                 aaaaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCACAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGcgcg                                                                                                                    <  2:338636/150‑1 (MQ=255)
                                   aaCGTCATTGAGCGCGCGGTATTGCTACCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGCACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTg                                                                                                                  <  2:1060517/150‑1 (MQ=255)
                                   aaCGTCATTGAGAGCGCGGTTTTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCAGCCTGCCGCAACGGTTGTCGACCTGGAGGGCGAAGATGATTATCAGTTGCTTGTGCGCGTg                                                                                                                  >  2:1628710/1‑150 (MQ=255)
                                     cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                                <  1:852183/150‑1 (MQ=255)
                                     cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                                <  2:434910/150‑1 (MQ=255)
                                        cATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                             <  1:1151109/150‑1 (MQ=255)
                                         aTTGAGCGCGCGGTATTGCTAACCCGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                            <  2:52939/150‑1 (MQ=255)
                                         aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                            >  1:1396081/1‑150 (MQ=255)
                                          ttGAGCGCGCGGTATTGCTAACTCGCGGTATCGTGCTGCCGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAAAGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAATGTGCGCATGCTgaaag                                                                                                           <  2:144135/150‑1 (MQ=255)
                                          ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaag                                                                                                           <  2:715341/150‑1 (MQ=255)
                                           tGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaaga                                                                                                          <  1:48903/150‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATTAATATCAGTAGATTGTGCGCGTGCTGAAAGAAAc                                                                                                       <  1:284092/150‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAc                                                                                                       <  2:1483275/150‑1 (MQ=255)
                                                gcgcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCa                                                                                                     >  2:1204041/1‑150 (MQ=255)
                                                    ggTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                                 <  1:626148/150‑1 (MQ=255)
                                                     gTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGc                                                                                                >  1:1601654/1‑150 (MQ=255)
                                                       aTTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGt                                                                                              >  2:773821/1‑150 (MQ=255)
                                                       aTTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGt                                                                                              <  1:609829/150‑1 (MQ=255)
                                                       aTTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGt                                                                                              >  1:544941/1‑150 (MQ=255)
                                                          gCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGt                                                                                           <  2:373716/150‑1 (MQ=255)
                                                          gCTAACACGCGGTAACGTGATGCAGCTGTCATAGCCAGATATACTTATACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGt                                                                                           >  1:1288811/1‑150 (MQ=255)
                                                                 cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACCCCGCCTGCCGCAACGGTTGTGGCCCAGGAGGCCGCAGATGAAATTCCGTAGATTGTGCGCTTGCTGAAAGAACCCACCGGGGTGGTTGCCggg                                                                                    >  2:289337/1‑150 (MQ=255)
                                                                  gcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCTGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGc                                                                                   <  1:734338/150‑1 (MQ=255)
                                                                  gcgGTAACGTGCTGCAGCTGGCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGc                                                                                   <  2:250720/150‑1 (MQ=255)
                                                                     gTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTa                                                                                >  1:849869/1‑150 (MQ=255)
                                                                       aaCGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTaaa                                                                              >  2:897103/1‑150 (MQ=255)
                                                                       aaCGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTaaa                                                                              <  2:652704/150‑1 (MQ=255)
                                                                        aCGGGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAg                                                                             <  2:1214852/150‑1 (MQ=255)
                                                                         cGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACTCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                            >  2:7200/1‑150 (MQ=255)
                                                                         cGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                            >  1:1267194/1‑150 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaag                                                                                                           <  2:385502/118‑1 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaag                                                                                                           >  1:385502/1‑118 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaag                                                                                                           <  2:168280/118‑1 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgc                                                                           <  1:851366/150‑1 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGATTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaag                                                                                                           >  1:168280/1‑118 (MQ=255)
                                                                            gctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgcgc                                                                         <  1:353114/150‑1 (MQ=255)
                                                                              tgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTg                                                                       >  1:175858/1‑150 (MQ=255)
                                                                                cAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCgcggcg                                                                     <  2:678323/150‑1 (MQ=255)
                                                                                     gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACg                                                                >  1:1652170/1‑150 (MQ=255)
                                                                                     gTCATTGCCAGATATTGTATTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTCGTTGCCGGGCCTAAAGGCGCTGCGCAACg                                                                >  2:420113/1‑150 (MQ=255)
                                                                                          tGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                           <  1:1223293/150‑1 (MQ=255)
                                                                                           gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                          <  1:773821/150‑1 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                         <  1:920040/150‑1 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGGCTgggg                                                         <  1:682561/150‑1 (MQ=255)
                                                                                             cAGATATTGTTTTCCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGTCGTGCTTGCCGGGCCTAAAGGCGATGCGCAACGTCTGGGGc                                                        >  2:1055653/1‑150 (MQ=255)
                                                                                               gATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTg                                                      <  2:1216440/150‑1 (MQ=255)
                                                                                                 tatTGTTTTCCCGGAACCTGAAACACCGCCTGCCGCAACGGCTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGAGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAACGCGCTGCGCAACGTCTGGGGCTGaa                                                    <  1:1145366/150‑1 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCa                                               <  1:1509796/150‑1 (MQ=255)
                                                                                                        ttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACg                                             <  1:1679046/150‑1 (MQ=255)
                                                                                                          aCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTATGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                           <  2:1543318/150‑1 (MQ=255)
                                                                                                          aCCGGAACCTGAAACGCCGCCTGCCGCAACGGTGGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTATGCGCGGGCTGACAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                           <  2:1550045/150‑1 (MQ=255)
                                                                                                           ccGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAcc                                          >  1:595633/1‑150 (MQ=255)
                                                                                                             ggAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGCCCct                                        >  2:1076518/1‑150 (MQ=255)
                                                                                                              gAACGTGACACGCCCCCTGCCGCAACGGTTGTCGCCCTGGAGGGCTATGATGAATATCAGTTGATTGGGCGCGTGCTGAATGAAACCAACGGCGTGGTGGCCGGGGCTAATGGCGCGGCGCAACGTCTGGGGCTGAAACGCACGACCctg                                       <  1:1420065/150‑1 (MQ=255)
                                                                                                              gAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctg                                       <  1:417293/150‑1 (MQ=255)
                                                                                                              gAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGCGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctg                                       >  2:1692758/1‑150 (MQ=255)
                                                                                                                 ccTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                    >  1:1151775/1‑150 (MQ=255)
                                                                                                                  cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                   >  2:25090/1‑150 (MQ=255)
                                                                                                                  cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                   <  1:1262345/150‑1 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                  >  1:601606/1‑150 (MQ=255)
                                                                                                                      aaCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACg                               >  1:890649/1‑150 (MQ=255)
                                                                                                                         gccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGCAACGCACGACCCTGCTGTCACGGAt                            >  1:1550/1‑150 (MQ=255)
                                                                                                                           ggggtGCCGCCACGGGTGTCGCCCTCGAGGGCGACGATGTCTCGCAGTGGTGTGGCCGCGTGCTGCCAGAAACCTACGGCGGGGTTGCCGGGCCTATAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGa                          <  1:289337/146‑1 (MQ=255)
                                                                                                                           cgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGa                          <  2:856381/150‑1 (MQ=255)
                                                                                                                              cTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGc                       <  1:58618/150‑1 (MQ=255)
                                                                                                                               tGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCg                      <  2:1267194/150‑1 (MQ=255)
                                                                                                                                  cGCACCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCt                   <  2:950734/150‑1 (MQ=255)
                                                                                                                                  cGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCt                   >  1:1207710/1‑150 (MQ=255)
                                                                                                                                   gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTg                  >  2:456537/1‑150 (MQ=255)
                                                                                                                                       cGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACACACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa              <  2:1601654/150‑1 (MQ=255)
                                                                                                                                        catcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGACTGGGAAt             >  2:220348/6‑150 (MQ=255)
                                                                                                                                         gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAGCGGCGTGGTTGCCGGGCCTAAAGGCCCTGCGCAACGTCTGGGGCTGAAACGCACGACCATGCTGTCAAGTATGAAGCGGCTGGGAAtt            >  1:544280/1‑150 (MQ=255)
                                                                                                                                         gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTATAGGCGCTGCGCACCGTCTGGGGCTGAAACGCACGACCCTGCTGTCAAGGATGAAGCGGCTGGGAAtt            >  1:866297/1‑150 (MQ=255)
                                                                                                                                         atcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGCCCCTGCTGTCACGGATGCAGCGGGTGGGAAtt            >  1:131676/5‑150 (MQ=255)
                                                                                                                                         atcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt            >  2:46599/5‑150 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACGCTGCTGTCACGGATGAAGCGGCTGGGAATTGa          >  1:1421541/3‑150 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa          >  2:476871/3‑150 (MQ=255)
                                                                                                                                             tCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTTTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATa        <  2:1425387/150‑1 (MQ=255)
                                                                                                                                             tCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGACTGGGAATTGATa        <  1:220348/150‑1 (MQ=255)
                                                                                                                                                cccTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAAt     >  1:1444299/1‑150 (MQ=255)
                                                                                                                                                 ccTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGTATGAAGCGGCTGGGAATTGATAAATc    <  1:260167/150‑1 (MQ=255)
                                                                                                                                                   tggtggGCTAATATGAATTTCTGTTGATTGTGAGCGTGATGAAAGAAAACAACGGCGTGGTTGAGGGGCCTAAAGGCGCTGCGCAACGTCTGTGGATGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGTAATTGATAAATctg  <  2:601606/150‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGG  >  NC_000913/2856109‑2856405

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