Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,930,602 A→G 5.2% D294D (GAT→GAC) ‡ peaA ← quinohemoprotein amine dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,930,6020AG 5.2% 750.1 / 18.1 230D294D (GAT→GAC) ‡peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (108/110);  new base G (6/6);  total (114/116)
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.26e-01

Only 100 of 232 total aligned reads displayed.

GCCTTGATGAACGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCTCATCACCCCAC  >  NC_002947/3930453‑3930750
                                                                                                                                                     |                                                                                                                                                     
  cTTGATGAACGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATc                                                                                                                                                     >  1:1009663/1‑150 (MQ=255)
          aCGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCGATGCCAACCt                                                                                                                                             <  1:4725434/150‑1 (MQ=255)
           cGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCtt                                                                                                                                            <  1:3596008/150‑1 (MQ=255)
           cGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCtt                                                                                                                                            <  2:1895465/150‑1 (MQ=255)
            gCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCGATGCCAACCtgc                                                                                                                                           <  1:3788338/150‑3 (MQ=255)
             ccGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTca                                                                                                                                          <  1:1263701/150‑1 (MQ=255)
             ccGGTTGTACCGCCAGCAACGGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTca                                                                                                                                          <  1:2864803/150‑1 (MQ=255)
              cGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTcac                                                                                                                                         >  2:4588278/1‑150 (MQ=255)
                     aCCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGccc                                                                                                                                 <  2:1800305/150‑1 (MQ=255)
                      ccGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGcccc                                                                                                                                >  2:327809/1‑150 (MQ=255)
                       cGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCGATGCCAAcc                                                                                                                                              >  1:1425942/1‑136 (MQ=255)
                       cGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCGATGCCAAcc                                                                                                                                              <  2:1425942/136‑1 (MQ=255)
                       cGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCg                                                                                                                               >  2:5593025/1‑150 (MQ=255)
                       cGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCg                                                                                                                               >  2:2804659/1‑150 (MQ=255)
                       cGCCAGCAACTGCGCCTTGCCTTCTTTTAAGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCg                                                                                                                               >  2:2804919/1‑150 (MQ=255)
                         ccAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGcc                                                                                                                             <  1:3189567/150‑1 (MQ=255)
                          cAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCa                                                                                                                            >  2:3749542/1‑150 (MQ=255)
                          cAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCa                                                                                                                            <  2:3786498/150‑1 (MQ=255)
                          cAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCa                                                                                                                            >  1:4343743/1‑150 (MQ=255)
                             cAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTc                                                                                                                         <  1:5505865/150‑1 (MQ=255)
                                 tGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCAGCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATaa                                                                                                                     >  1:2066099/1‑150 (MQ=255)
                                   cgcCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAAcc                                                                                                                   >  2:4063762/1‑150 (MQ=255)
                                     ccttgccttCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGt                                                                                                                 <  1:259501/150‑1 (MQ=255)
                                     ccttgccttCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGt                                                                                                                 >  1:2548710/1‑150 (MQ=255)
                                      cttgccttCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGtt                                                                                                                <  2:5131038/150‑1 (MQ=255)
                                       atgcctaCATTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTCCCCCGCCATTCATAACCGTTg                                                                                                               <  2:2648267/149‑1 (MQ=255)
                                          ccttATTTTTCGGCGCTGAAATCCAGGACGCGCTCATCGTGATCGGCTTCGAACATTCGGCCTTTCATCTCGCCATCCAGAGCGTCGAACACCTGGCGCAGGTTGGCATAGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTAc                                                                                                            <  2:2321159/150‑1 (MQ=255)
                                           cttcttTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACa                                                                                                           <  2:5619948/150‑1 (MQ=255)
                                            ttcttTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACAGTGCCCCGCCCTTCATAACCGTTGTACAg                                                                                                          >  2:3413986/1‑150 (MQ=255)
                                               ttttACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGAt                                                                                                       >  1:4910131/1‑150 (MQ=255)
                                                  tACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGc                                                                                                    <  1:562768/150‑1 (MQ=255)
                                                     ggCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGa                                                                                                 <  1:5556288/150‑1 (MQ=255)
                                                          tGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGc                                                                                            <  2:4001917/150‑1 (MQ=255)
                                                           gAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCt                                                                                           >  2:3460243/1‑150 (MQ=255)
                                                            aaaTCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTg                                                                                          >  2:5356666/1‑150 (MQ=255)
                                                             aaTCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGc                                                                                         <  1:3552754/150‑1 (MQ=255)
                                                             aaTCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCAGCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGc                                                                                         <  1:351341/150‑1 (MQ=255)
                                                               tCCAGGCCGCGCTCGTCGTGATCGGCTGCGAACATACGGCCTGTCCTCGCGCCATCCGGCGCGGCGAACCCCTGGCGCAGGTTGGCATCGCCTCCCTGCACAGTGCCCCGCCCTGCATAGCCGGTGTCCAGGCTGGCCGAGCCGCTGCCg                                                                                       >  2:4927526/1‑150 (MQ=255)
                                                               tCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCg                                                                                       >  2:2096169/1‑150 (MQ=255)
                                                                 cAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGtt                                                                                     >  1:466091/1‑150 (MQ=255)
                                                                    gCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGaa                                                                                  <  2:2540077/150‑1 (MQ=255)
                                                                     ccGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAAc                                                                                 <  1:3486944/150‑1 (MQ=255)
                                                                        cgcTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGgt                                                                              <  2:493728/150‑1 (MQ=255)
                                                                        cgcTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGgt                                                                              <  1:2113668/150‑1 (MQ=255)
                                                                             aTCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTAcc                                                                         <  1:1418611/150‑1 (MQ=255)
                                                                                   aTCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGCGTACCGTCGGc                                                                   <  1:1114401/150‑1 (MQ=255)
                                                                                       gCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTaa                                                               >  1:5158074/1‑150 (MQ=255)
                                                                                       gCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTaa                                                               <  2:4230667/150‑1 (MQ=255)
                                                                                       gCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTaa                                                               >  1:4432075/1‑150 (MQ=255)
                                                                                         ttCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAgc                                                             >  1:3688492/1‑150 (MQ=255)
                                                                                           cGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTc                                                                                                                         >  2:2733231/1‑88 (MQ=255)
                                                                                            gAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGcc                                                          >  2:2224029/1‑150 (MQ=255)
                                                                                             aaCATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGccc                                                         >  1:4489753/1‑150 (MQ=255)
                                                                                              aCATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCt                                                        <  2:360071/150‑1 (MQ=255)
                                                                                              aCATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCt                                                        <  1:3997980/150‑1 (MQ=255)
                                                                                              aCATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGAGGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCt                                                        <  2:360155/150‑1 (MQ=255)
                                                                                               cATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCtt                                                       <  2:4489753/150‑1 (MQ=255)
                                                                                               cATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCtt                                                       <  1:3963838/150‑1 (MQ=255)
                                                                                               cATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCtt                                                       <  1:697390/150‑1 (MQ=255)
                                                                                                aTACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTg                                                      >  1:2973512/1‑150 (MQ=255)
                                                                                                aTACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTg                                                      <  1:2282749/150‑1 (MQ=255)
                                                                                                aTACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCGTg                                                      >  1:3054912/1‑150 (MQ=255)
                                                                                                  aCGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGAc                                                    <  2:2476271/150‑1 (MQ=255)
                                                                                                   cGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGAcc                                                   <  2:1378152/150‑1 (MQ=255)
                                                                                                    ggCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCt                                                  <  2:1304413/150‑1 (MQ=255)
                                                                                                    ggCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCGTGAccg                                                  >  1:4067160/1‑149 (MQ=255)
                                                                                                     gCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTc                                                 >  2:4872044/1‑150 (MQ=255)
                                                                                                     gCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTc                                                 <  2:5322217/150‑1 (MQ=255)
                                                                                                     gCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTc                                                 >  2:3073949/1‑150 (MQ=255)
                                                                                                     gCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTc                                                 <  1:1446071/150‑1 (MQ=255)
                                                                                                           cATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCtt                                           <  2:1493655/150‑1 (MQ=255)
                                                                                                            aTCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTg                                          >  1:3025678/1‑150 (MQ=255)
                                                                                                             tctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGa                                         <  2:4775844/150‑1 (MQ=255)
                                                                                                             tctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGa                                         <  1:4581137/150‑1 (MQ=255)
                                                                                                             tctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGa                                         >  1:1762668/1‑150 (MQ=255)
                                                                                                              ctcGCCATCCAGTGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGaa                                        >  2:208327/1‑150 (MQ=255)
                                                                                                              ctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGaa                                        <  1:537566/150‑1 (MQ=255)
                                                                                                              ctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGaa                                        <  2:4376708/150‑1 (MQ=255)
                                                                                                              ctcGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGaa                                        <  2:1067767/150‑1 (MQ=255)
                                                                                                                   cATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTAt                                   <  2:462990/150‑1 (MQ=255)
                                                                                                                   cATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCAACGCCAAACTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGAGTAACGTAGGCGTAACAGCCCTTGACCTCGAACTAGAATGTAt                                   >  1:613515/1‑150 (MQ=255)
                                                                                                                    aTCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATc                                  >  2:3552628/1‑150 (MQ=255)
                                                                                                                      ccAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGc                                <  2:5222052/150‑1 (MQ=255)
                                                                                                                            cggcgAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCtctc                          <  2:2173432/150‑1 (MQ=255)
                                                                                                                            cggcgAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCtctc                          <  2:300671/150‑1 (MQ=255)
                                                                                                                                   cacCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGt                   <  2:658812/150‑1 (MQ=255)
                                                                                                                                    acCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCGTGCAGGTCTCGCCCTCTCCGGCCGt                   >  2:5460856/1‑150 (MQ=255)
                                                                                                                                      cTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGAc                >  2:1148172/1‑150 (MQ=255)
                                                                                                                                      cTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGAc                >  1:1367755/1‑150 (MQ=255)
                                                                                                                                      cTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGAAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGAc                >  2:1042763/1‑150 (MQ=255)
                                                                                                                                        ggCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGc              <  2:1669213/150‑1 (MQ=255)
                                                                                                                                        ggCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGc              >  1:2644800/1‑150 (MQ=255)
                                                                                                                                         gcgcAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCt             <  2:101776/150‑1 (MQ=255)
                                                                                                                                           gcAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCtca           >  2:5085364/1‑150 (MQ=255)
                                                                                                                                            aaGGTTGGCGATGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCtcat          >  1:2931724/2‑150 (MQ=255)
                                                                                                                                            aaGGTTGGCGATGCCAACCGTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCtcat          >  2:2716186/2‑150 (MQ=255)
                                                                                                                                             aGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCtcatc         >  1:1937743/1‑150 (MQ=255)
                                                                                                                                               gTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCTCATCAc       <  1:574203/150‑1 (MQ=255)
                                                                                                                                               gTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCTCATCAc       >  1:2850845/1‑150 (MQ=255)
                                                                                                                                                  ggCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCTCATCAcccc    <  1:1726805/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTTGATGAACGCCGGTTGTACCGCCAGCAACTGCGCCTTGCCTTCTTTTACGGCGCTGAAATCCAGGCCGCGCTCATCGTGATCGGCTTCGAACATACGGCCTTTCATCTCGCCATCCAGCGCGGCGAACACCTGGCGCAGGTTGGCATCGCCAACCTTCACATTGCCCCGCCATTCATAACCGTTGTACAGGATGGCCGAGCCGCTGCCGTTGAACGGTGTACCGTCGGCGTAAGCGCCCTTGACCTCGACCTTGAAGGTATCGCCCTCTCCGGCCGTGACGCTCATCACCCCAC  >  NC_002947/3930453‑3930750

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