Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,960,777 A→G S4P (TCC→CCC) 
L134L (CTT→CTC
ppdC ←
b2824 ←
prepilin peptidase dependent protein C
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,7770AG82.5% 619.8 / 134.7 291S4P (TCC→CCC) 
L134L (CTT→CTC
ppdC
b2824
prepilin peptidase dependent protein C
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (129/111);  total (155/136)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 295 total aligned reads displayed.

TGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTT  >  NC_000913/2960686‑2960867
                                                                                           |                                                                                          
tGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGaa                                                                                           >  2:2399871/1‑93 (MQ=255)
 ggAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAg                                                                                          <  1:7447071/93‑1 (MQ=255)
 ggAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAg                                                                                          >  1:2632142/1‑93 (MQ=255)
     aCCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGa                                                                                      >  2:297426/1‑93 (MQ=255)
     aCCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGa                                                                                      >  1:864703/1‑93 (MQ=255)
        cGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACAt                                                                                   >  1:5292656/1‑93 (MQ=255)
        cGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACAt                                                                                   >  1:3465097/1‑93 (MQ=255)
        cGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACAt                                                                                   >  2:1539712/1‑93 (MQ=255)
         gATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATa                                                                                  >  2:6830202/1‑93 (MQ=255)
         gATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATa                                                                                  >  2:1116500/1‑93 (MQ=255)
             aCGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGc                                                                              >  2:3622557/1‑93 (MQ=255)
              cGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGcc                                                                             <  2:6480818/93‑1 (MQ=255)
               gCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGccc                                                                            >  2:254121/1‑93 (MQ=255)
                  gTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCctct                                                                         >  1:1466165/1‑93 (MQ=255)
                  gTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCctct                                                                         >  2:1188385/1‑93 (MQ=255)
                   tGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCctctc                                                                        <  1:4783808/93‑1 (MQ=255)
                   tGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCctctc                                                                        <  1:1426996/93‑1 (MQ=255)
                    gACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCctctct                                                                       <  2:192545/93‑1 (MQ=255)
                     aCAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCtt                                                                      >  2:217102/1‑93 (MQ=255)
                       aaTCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTc                                                                    <  2:2225195/93‑1 (MQ=255)
                       aaTCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTc                                                                    <  1:5615526/93‑1 (MQ=255)
                       aaTCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTc                                                                    <  1:6557896/93‑1 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCa                                                                   <  2:505369/93‑1 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCa                                                                   >  1:1143334/1‑93 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCa                                                                   <  1:3902839/93‑1 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCa                                                                   <  1:1314794/93‑1 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCa                                                                   <  1:4506449/93‑1 (MQ=255)
                        atcatcACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCa                                                                   <  1:1703569/93‑1 (MQ=255)
                         tcatcaCCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAg                                                                  >  1:212384/1‑93 (MQ=255)
                            tcaCCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCgg                                                               >  2:3947454/1‑93 (MQ=255)
                             caCCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGa                                                              >  2:4530126/1‑93 (MQ=255)
                               ccATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACa                                                            >  1:2750985/1‑93 (MQ=255)
                                 aTCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACaaa                                                          >  2:5789370/1‑93 (MQ=255)
                                  tCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACaaaa                                                         >  1:3019503/1‑93 (MQ=255)
                                  tCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACaaaa                                                         >  1:3126239/1‑93 (MQ=255)
                                     acacCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATc                                                      >  1:5444451/1‑93 (MQ=255)
                                      cacCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCg                                                     >  2:128648/1‑93 (MQ=255)
                                        ccATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCt                                                   <  2:1751638/93‑1 (MQ=255)
                                        ccATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCt                                                   <  2:4829064/93‑1 (MQ=255)
                                         cATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTc                                                  >  1:4983318/1‑93 (MQ=255)
                                          aTCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTcc                                                 >  2:729564/1‑93 (MQ=255)
                                          aTCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTcc                                                 <  1:6781985/93‑1 (MQ=255)
                                          aTCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCCGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTcc                                                 <  1:7305676/93‑1 (MQ=255)
                                           tCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCa                                                >  2:3853899/1‑93 (MQ=255)
                                           tCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCa                                                <  2:6705078/93‑1 (MQ=255)
                                           tCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCa                                                >  2:3766230/1‑93 (MQ=255)
                                           tCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCa                                                >  2:2728832/1‑93 (MQ=255)
                                           tCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGACATCTCTTTCAGCGGACAAAAATCGCTCCa                                                >  1:4850829/1‑93 (MQ=255)
                                             gCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGc                                              <  2:6615142/93‑1 (MQ=255)
                                             gCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGc                                              <  2:3385749/93‑1 (MQ=255)
                                               cAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcg                                            >  2:284129/1‑93 (MQ=255)
                                               cAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcg                                            >  2:7018360/1‑93 (MQ=255)
                                                aaCATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcgc                                           >  2:6289364/1‑93 (MQ=255)
                                                aaCATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcgc                                           >  2:6956229/1‑93 (MQ=255)
                                                 aCATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcgcg                                          <  2:4011784/93‑1 (MQ=255)
                                                 aCATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCcgcg                                          <  1:6328927/93‑1 (MQ=255)
                                                  cATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGt                                         >  1:6033893/1‑93 (MQ=255)
                                                   aTTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTg                                        >  1:517209/1‑93 (MQ=255)
                                                   aTTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTg                                        <  1:4890345/93‑1 (MQ=255)
                                                    ttACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTgg                                       >  1:2943434/1‑93 (MQ=255)
                                                    ttACCGCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTgg                                       >  1:4701044/1‑93 (MQ=255)
                                                     tACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGc                                      >  2:3263208/1‑93 (MQ=255)
                                                        cTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaa                                   >  1:2184976/1‑93 (MQ=255)
                                                         tCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaaa                                  <  2:226912/93‑1 (MQ=255)
                                                         tCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaaa                                  >  1:3692143/1‑93 (MQ=255)
                                                            ggCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaaaaaa                               <  1:5688531/93‑1 (MQ=255)
                                                            ggCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaaaaaa                               <  1:7217121/93‑1 (MQ=255)
                                                            ggCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGaaaaaa                               >  1:1230432/1‑93 (MQ=255)
                                                             gCAGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAAc                              <  2:5294901/93‑1 (MQ=255)
                                                               aGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAGCaa                            >  2:4139264/1‑93 (MQ=255)
                                                               aGGCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACaa                            >  1:2671251/1‑93 (MQ=255)
                                                                ggCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAAt                           >  2:5710366/1‑93 (MQ=255)
                                                                ggCTAAAGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAAt                           <  2:6322732/93‑1 (MQ=255)
                                                                  cTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGt                         <  1:2019834/93‑1 (MQ=255)
                                                                     aaGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTc                      <  1:4346823/93‑1 (MQ=255)
                                                                     aaGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTc                      >  1:6238644/1‑93 (MQ=255)
                                                                     aaGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTc                      <  2:1042152/93‑1 (MQ=255)
                                                                      aGCCTTGTTGATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTcc                     >  1:1572030/1‑93 (MQ=255)
                                                                      aGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTcc                     >  2:4663300/1‑93 (MQ=255)
                                                                          ttgttgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCa                 >  1:1007905/1‑93 (MQ=255)
                                                                           tgttgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCaa                <  1:4072011/93‑1 (MQ=255)
                                                                           tgttgATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCaa                <  1:774238/93‑1 (MQ=255)
                                                                            gttgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAAt               <  2:143709/93‑1 (MQ=255)
                                                                            gttgATTCTTCAGGAGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAAt               >  2:6439485/1‑93 (MQ=255)
                                                                              tgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGa             <  2:4242345/93‑1 (MQ=255)
                                                                              tgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGa             <  1:2346251/93‑1 (MQ=255)
                                                                              tgATTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGa             <  2:6089608/93‑1 (MQ=255)
                                                                               gATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGAc            <  1:5794629/93‑1 (MQ=255)
                                                                                aTTCTTCAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACt           <  2:1939039/93‑1 (MQ=255)
                                                                                 ttcttcAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACtt          <  2:776343/93‑1 (MQ=255)
                                                                                  tcttcAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTc         <  1:2253751/93‑1 (MQ=255)
                                                                                  tcttcAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTc         <  1:6890106/93‑1 (MQ=255)
                                                                                  tcttcAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTc         <  1:3371093/93‑1 (MQ=255)
                                                                                  tcttcAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTc         >  2:4211193/1‑93 (MQ=255)
                                                                                  tcttcAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTc         >  1:3431350/1‑93 (MQ=255)
                                                                                    ttcAGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGc       <  2:463540/93‑1 (MQ=255)
                                                                                     tcAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGcc      >  1:3633088/1‑93 (MQ=255)
                                                                                       aGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTg    >  1:6938451/1‑93 (MQ=255)
                                                                                       aGGGGAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTg    >  1:5001968/1‑93 (MQ=255)
                                                                                        ggggAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGt   <  1:4158471/93‑1 (MQ=255)
                                                                                           |                                                                                          
TGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGCCAACATTACCTCCGGCAGGCTAAAGCCTTGTTGATTCTTCAGGGAAGCTGACATAACGCCCTCTCTTTCAGCGGACAAAAATCGCTCCAGCCGCGTGGCGAAAAAACAATGTTTCCATCAATGACTTCGCCTGTT  >  NC_000913/2960686‑2960867

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