Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 250,749 T→C G236G (GGT→GGC mbhA → putative motility protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913250,7490TC85.2% 669.0 / 104.4 284G236G (GGT→GGCmbhAputative motility protein
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base C (109/133);  total (132/152)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 285 total aligned reads displayed.

ATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATCTTTA  >  NC_000913/250658‑250841
                                                                                           |                                                                                            
   ccAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGc                                                                                          >  2:2114335/1‑93 (MQ=255)
   ccAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGc                                                                                          >  2:7307282/1‑93 (MQ=255)
   ccAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGc                                                                                          >  2:6744190/1‑93 (MQ=255)
     aaaaaTCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCAt                                                                                        <  1:7105281/93‑1 (MQ=255)
     aaaaaTCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCAt                                                                                        >  2:5528163/1‑93 (MQ=255)
       aaaTCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATgg                                                                                      >  1:347992/1‑93 (MQ=255)
         aTCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATgggg                                                                                    <  2:5413418/93‑1 (MQ=255)
          tCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGa                                                                                   >  2:2870042/1‑93 (MQ=255)
            cGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATa                                                                                 >  1:5604010/1‑93 (MQ=255)
                aaGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAgt                                                                             >  1:100140/1‑93 (MQ=255)
                aaGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAgt                                                                             >  2:182028/1‑93 (MQ=255)
                  gcgcGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAgtgg                                                                           <  1:3382168/93‑1 (MQ=255)
                   cgcgGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAgtggt                                                                          <  2:1067373/93‑1 (MQ=255)
                     cgGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGc                                                                        <  1:5257670/93‑1 (MQ=255)
                     cgGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGc                                                                        <  1:4597872/93‑1 (MQ=255)
                      gggCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCa                                                                       <  2:7035988/93‑1 (MQ=255)
                      gggCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCa                                                                       >  2:6357338/1‑93 (MQ=255)
                       ggCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAg                                                                      >  1:2647059/1‑93 (MQ=255)
                       ggCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAg                                                                      <  1:1915797/93‑1 (MQ=255)
                       ggCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAg                                                                      >  2:3795803/1‑93 (MQ=255)
                        gCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGc                                                                     <  2:613688/93‑1 (MQ=255)
                        gCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGc                                                                     <  1:7294252/93‑1 (MQ=255)
                        gCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGc                                                                     <  1:5838059/93‑1 (MQ=255)
                        gCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGc                                                                     <  2:2190248/93‑1 (MQ=255)
                          aaCCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCg                                                                   >  1:7488974/1‑93 (MQ=255)
                           aCCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGc                                                                  >  1:2135454/1‑93 (MQ=255)
                            ccGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCt                                                                 <  2:5907174/93‑1 (MQ=255)
                              ggCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTgg                                                               >  2:1158599/1‑93 (MQ=255)
                               gcgcATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGt                                                              >  1:3788343/1‑93 (MQ=255)
                               gcgcATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGt                                                              <  1:1979893/93‑1 (MQ=255)
                                cgcATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTg                                                             <  1:5993181/93‑1 (MQ=255)
                                 gcATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGc                                                            <  1:304418/93‑1 (MQ=255)
                                  cATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCa                                                           >  2:4544835/1‑93 (MQ=255)
                                    ttGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCaaa                                                         >  2:4907055/1‑93 (MQ=255)
                                    ttGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCaaa                                                         >  2:1319934/1‑93 (MQ=255)
                                    ttGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCaaa                                                         >  2:2313179/1‑93 (MQ=255)
                                    ttGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCaaa                                                         >  2:468460/1‑93 (MQ=255)
                                     tGAGATTATGGCGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCaaaa                                                        <  1:6501181/93‑1 (MQ=255)
                                      gagaTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAg                                                       <  2:3477459/93‑1 (MQ=255)
                                      gagaTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAg                                                       >  1:5167071/1‑93 (MQ=255)
                                        gaTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCt                                                     <  2:2968889/93‑1 (MQ=255)
                                        gaTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCt                                                     >  1:2077254/1‑93 (MQ=255)
                                        gaTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCt                                                     >  1:4304288/1‑93 (MQ=255)
                                         aTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTg                                                    >  1:6860975/1‑93 (MQ=255)
                                         aTTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTg                                                    <  1:5680941/93‑1 (MQ=255)
                                          ttATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTgg                                                   <  1:5363059/93‑1 (MQ=255)
                                          ttATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTgg                                                   <  2:3189585/93‑1 (MQ=255)
                                          ttATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTgg                                                   >  2:4368135/1‑93 (MQ=255)
                                          ttATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTgg                                                   >  1:2849655/1‑93 (MQ=255)
                                                tGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAagc                                             <  1:89614/93‑1 (MQ=255)
                                                tGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAagc                                             >  2:2678633/1‑93 (MQ=255)
                                                 gCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAagca                                            <  2:526182/93‑1 (MQ=255)
                                                   tGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAagcagc                                          >  1:4962531/1‑93 (MQ=255)
                                                   tGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAagcagc                                          >  2:1274648/1‑93 (MQ=255)
                                                    gACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAgcagca                                         <  2:5768514/93‑1 (MQ=255)
                                                    gACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAgcagca                                         <  1:2924504/93‑1 (MQ=255)
                                                    gACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAgcagca                                         >  2:6228142/1‑93 (MQ=255)
                                                     aCCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGcagcag                                        <  2:5042596/93‑1 (MQ=255)
                                                         aaaGTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGc                                    >  2:3749734/1‑93 (MQ=255)
                                                            gTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCttt                                 <  2:1899039/93‑1 (MQ=255)
                                                            gTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCttt                                 <  1:4834490/93‑1 (MQ=255)
                                                            gTGCGTCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCttt                                 <  1:2342469/93‑1 (MQ=255)
                                                                 tCCGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAg                            <  2:227120/93‑1 (MQ=255)
                                                                 tCCGATACGCTGTATAAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAg                            >  2:1615159/1‑93 (MQ=255)
                                                                  ccGATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGc                           <  2:5649107/93‑1 (MQ=255)
                                                                    gATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGa                         <  1:16971/93‑1 (MQ=255)
                                                                    gATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGa                         >  1:631528/1‑93 (MQ=255)
                                                                    gATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGa                         >  2:4451314/1‑93 (MQ=255)
                                                                    gATACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGa                         >  1:1260659/1‑93 (MQ=255)
                                                                      tACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAc                       >  2:86620/1‑93 (MQ=255)
                                                                      tACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAc                       <  2:271428/93‑1 (MQ=255)
                                                                      tACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAc                       <  1:3629545/93‑1 (MQ=255)
                                                                      tACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAc                       <  1:399724/93‑1 (MQ=255)
                                                                      tACGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAc                       >  2:4935001/1‑93 (MQ=255)
                                                                        cGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAcgc                     <  1:1581085/93‑1 (MQ=255)
                                                                        cGCTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAcgc                     <  2:1828350/93‑1 (MQ=255)
                                                                         gCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAAcgcg                    <  2:4226320/93‑1 (MQ=255)
                                                                          cTGTATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGt                   >  2:4821688/1‑93 (MQ=255)
                                                                             tATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaa                >  2:3948645/1‑93 (MQ=255)
                                                                             tATCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaa                <  2:324858/93‑1 (MQ=255)
                                                                              ttCAATGCTTGGGTCAGCATGGGGATAAAGTGGTGCAGCTGCTGGTTCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaaa               <  2:2598380/92‑1 (MQ=255)
                                                                              aTCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaaa               >  2:6289419/1‑93 (MQ=255)
                                                                              aTCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaaa               >  1:269129/1‑93 (MQ=255)
                                                                              aTCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTaaa               <  1:5484791/93‑1 (MQ=255)
                                                                               tCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAAt              >  1:4466472/1‑93 (MQ=255)
                                                                               tCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAAt              <  2:1996408/93‑1 (MQ=255)
                                                                               tCAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAAt              >  1:2482739/1‑93 (MQ=255)
                                                                                cAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATg             <  1:6079888/93‑1 (MQ=255)
                                                                                cAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATg             >  1:1632760/1‑93 (MQ=255)
                                                                                cAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATg             >  1:1319351/1‑93 (MQ=255)
                                                                                cAATACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATg             >  2:6459507/1‑93 (MQ=255)
                                                                                 aaTACTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGc            <  2:7198836/93‑1 (MQ=255)
                                                                                   tACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTg          >  2:6713512/1‑93 (MQ=255)
                                                                                    aCTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAACTGCTGa         <  2:4723240/93‑1 (MQ=255)
                                                                                    aCTTTGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGa         <  2:6592669/93‑1 (MQ=255)
                                                                                      tttGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAAt       >  1:366688/1‑93 (MQ=255)
                                                                                       ttGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATc      <  2:6196741/93‑1 (MQ=255)
                                                                                       ttGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATc      <  2:3623352/93‑1 (MQ=255)
                                                                                        tGGCCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATCt     <  2:6860975/93‑1 (MQ=255)
                                                                                          gTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATCttt   <  2:918155/93‑1 (MQ=255)
                                                                                           |                                                                                            
ATTCCAAAAATCCGCAAAGCGCGGGCAACCGGCGCATTGAGATTATGGTGCTGACCAAAAGTGCGTCCGATACGCTGTATCAATACTTTGGTCAGCATGGGGATAAAGTGGTGCAGCCGCTGGTGCAAAAGCTGGATAAGCAGCAGGTGCTTTCGCAGCGAACGCGTTAAATGCTGAATCTTTA  >  NC_000913/250658‑250841

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