Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,449,280 T→C T109A (ACC→GCC)  b2335 ← putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2800TC98.6% 449.3 / ‑2.6 144T109A (ACC→GCC) b2335putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (72/70);  total (72/72)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 148 total aligned reads displayed.

TTATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGC  >  NC_000913/2449192‑2449369
                                                                                        |                                                                                         
ttATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGata                                                                                       >  2:2829488/1‑93 (MQ=255)
ttATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGata                                                                                       >  2:3128187/1‑93 (MQ=255)
ttATTGATATTAAATGGAGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGata                                                                                       <  2:3954867/93‑1 (MQ=255)
  aTTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATAtt                                                                                     <  2:694182/93‑1 (MQ=255)
    tGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTgg                                                                                   >  2:2463222/1‑93 (MQ=255)
    tGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTgg                                                                                   <  2:1698063/93‑1 (MQ=255)
    tGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCACAGCCCGGCGATATTgg                                                                                   >  1:2008798/1‑93 (MQ=255)
      atatTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTGTCCGAATGCCCAGCCCGGCGATATTGGtc                                                                                 <  1:4207114/93‑2 (MQ=255)
      atatTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTa                                                                                 <  2:3440247/93‑1 (MQ=255)
         ttAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTg                                                                              >  1:1020067/1‑93 (MQ=255)
           taaTGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAg                                                                            <  2:1316765/92‑1 (MQ=255)
           aaaTGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAg                                                                            <  2:1451032/93‑1 (MQ=255)
              tGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAc                                                                         >  2:822196/1‑93 (MQ=255)
              tGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAc                                                                         <  2:2775187/93‑1 (MQ=255)
               ggCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAcc                                                                        >  1:2760673/1‑93 (MQ=255)
                 cGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGt                                                                      <  1:1579124/93‑1 (MQ=255)
                 cGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGt                                                                      <  2:3186732/93‑1 (MQ=255)
                  gTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGtt                                                                     >  2:263373/1‑93 (MQ=255)
                   tCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGttt                                                                    >  1:2536111/1‑93 (MQ=255)
                   tCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGttt                                                                    >  2:1279949/1‑93 (MQ=255)
                   gcccgCTATTTTCACCGACGGCAAACAATGAGTTAGCGGCGCTATTTTCTATCCGAAGGCCCAGCCCGGCGATATTGGTACTGAGTACCGttt                                                                    <  1:2300335/89‑1 (MQ=255)
                    ccAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTc                                                                   <  2:4712685/93‑1 (MQ=255)
                    ccAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTc                                                                   >  2:3414204/1‑93 (MQ=255)
                     cAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCg                                                                  >  1:321098/1‑93 (MQ=255)
                     cAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCAGCGATATTGGTACTGAGTACCGTTTCg                                                                  >  1:3286772/1‑93 (MQ=255)
                         tATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGt                                                              >  1:2566026/1‑93 (MQ=255)
                          aTTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCgct                                                             >  1:707323/1‑91 (MQ=255)
                          aTTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGtt                                                             <  2:2149531/93‑1 (MQ=255)
                          aTTTTAACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGtt                                                             >  2:3701333/1‑93 (MQ=255)
                           ttttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTa                                                            >  1:3592183/1‑93 (MQ=255)
                           ttttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTa                                                            >  1:1758621/1‑93 (MQ=255)
                           ttttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTa                                                            >  1:4141400/1‑93 (MQ=255)
                            tttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTTCCGTTTCGCCGTTaa                                                           <  2:3884984/93‑1 (MQ=255)
                             ttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAAt                                                          <  1:547758/93‑1 (MQ=255)
                                cccGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGt                                                       <  2:3316335/92‑1 (MQ=255)
                                aCCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGt                                                       >  2:3267845/1‑93 (MQ=255)
                                aCCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGt                                                       >  2:3267843/1‑93 (MQ=255)
                                aCCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGt                                                       <  2:2863279/93‑1 (MQ=255)
                                aCCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGt                                                       >  1:3169745/1‑93 (MQ=255)
                                 ccGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGtt                                                      <  2:1234652/93‑1 (MQ=255)
                                   gACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGa                                                    <  2:2017082/93‑1 (MQ=255)
                                   gACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGa                                                    <  2:4231813/93‑1 (MQ=255)
                                     cGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAgg                                                  >  2:1824940/1‑93 (MQ=255)
                                      ggCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCACCGTTAATCGTTGAGGt                                                 >  2:495592/1‑93 (MQ=255)
                                        cAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCGGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAt                                               <  2:1455714/93‑1 (MQ=255)
                                        cAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAt                                               <  2:1340177/93‑1 (MQ=255)
                                        cAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAt                                               <  2:81638/93‑1 (MQ=255)
                                         aaaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAtt                                              >  1:1332731/1‑93 (MQ=255)
                                          aaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTg                                             <  2:4373585/93‑1 (MQ=255)
                                            cAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGcc                                           <  1:1621306/93‑1 (MQ=255)
                                            cAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGcc                                           >  1:2614671/1‑93 (MQ=255)
                                             aaTGAGTTGTCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGccc                                          <  2:4665684/93‑1 (MQ=255)
                                             aaTGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGccc                                          >  1:614046/1‑93 (MQ=255)
                                              aTGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCt                                         <  1:3948709/93‑1 (MQ=255)
                                               tGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCtt                                        <  1:3700061/93‑1 (MQ=255)
                                                gAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTa                                       >  1:3942827/1‑93 (MQ=255)
                                                 aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  1:4270987/93‑1 (MQ=255)
                                                  gTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAg                                     >  1:1179394/1‑93 (MQ=255)
                                                    tATCGGCGCTATTTTCTATCCGAATGCCCCGCCCGGCGATATTGGTACTGAGCACCGGTTCGCCGTTAATCGTTGAGGTATTGTCCTTAAgct                                   >  1:1751734/1‑91 (MQ=255)
                                                    tATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGtt                                   >  1:3413210/1‑93 (MQ=255)
                                                    tATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGtt                                   >  1:4297722/1‑93 (MQ=255)
                                                    tATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGAATTGTTCTTAAGtt                                   >  2:4390320/1‑93 (MQ=255)
                                                    tATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATAGTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGATGTATTGCCCTTAAGtt                                   >  1:3666137/1‑93 (MQ=255)
                                                       cGGCGCTATTTTCTATCCGACTGCCCTGCCCGGCGATATTGGTACTGCGTCCCGTTTCGCCGTTCCTCGTTGAGGTATTGCCCTTAAGTTGCa                                <  1:4200929/93‑1 (MQ=255)
                                                         gcgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATc                              <  1:3239240/93‑1 (MQ=255)
                                                         gcgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATc                              <  2:268881/93‑1 (MQ=255)
                                                         gcgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATc                              >  1:1635252/1‑93 (MQ=255)
                                                          cgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATcc                             <  2:3561159/93‑1 (MQ=255)
                                                             tATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCa                          >  2:2226872/1‑93 (MQ=255)
                                                             tATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCa                          <  2:683122/93‑1 (MQ=255)
                                                              aTTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCaa                         >  1:1399427/1‑93 (MQ=255)
                                                                   cTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGt                    <  2:584477/93‑1 (MQ=255)
                                                                     aTCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTa                  <  2:2706595/93‑1 (MQ=255)
                                                                     aTCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTa                  >  2:1871311/1‑93 (MQ=255)
                                                                     aTCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTa                  >  2:4184163/1‑93 (MQ=255)
                                                                     aTCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTa                  <  2:3720273/93‑1 (MQ=255)
                                                                       ccGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGc                <  1:3797710/93‑1 (MQ=255)
                                                                        cGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCa               >  1:742482/1‑93 (MQ=255)
                                                                        cGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCa               >  2:3994898/1‑93 (MQ=255)
                                                                        cGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCa               >  1:1923293/1‑93 (MQ=255)
                                                                         gAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCaa              >  1:4347916/1‑93 (MQ=255)
                                                                         gAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCaa              <  1:2632033/93‑1 (MQ=255)
                                                                         gAATGCACAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCaa              <  2:2864240/93‑1 (MQ=255)
                                                                          aaTGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAg             >  1:930554/1‑93 (MQ=255)
                                                                            tGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCGAGAc           <  1:1889607/93‑1 (MQ=255)
                                                                            tGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGAc           >  1:4666394/1‑93 (MQ=255)
                                                                              cccAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTg         <  1:1974586/93‑1 (MQ=255)
                                                                               ccAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGt        <  1:1223220/93‑1 (MQ=255)
                                                                               ccAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGt        <  2:1313114/93‑1 (MQ=255)
                                                                                cAGCCCGTCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:633221/93‑1 (MQ=255)
                                                                                cAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:2080530/93‑1 (MQ=255)
                                                                                cAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:3380573/93‑1 (MQ=255)
                                                                                cAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:4169686/93‑1 (MQ=255)
                                                                                 aGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGttt      <  2:907165/93‑1 (MQ=255)
                                                                                  gCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtg     >  1:2820621/1‑93 (MQ=255)
                                                                                   cccGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtgt    <  1:847440/93‑1 (MQ=255)
                                                                                   cccGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtgt    <  1:3062209/93‑1 (MQ=255)
                                                                                    ccGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtgtg   >  1:155406/1‑93 (MQ=255)
                                                                                     cGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGc  <  1:4380619/93‑1 (MQ=255)
                                                                                     cGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGc  >  2:877445/1‑93 (MQ=255)
                                                                                        |                                                                                         
TTATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGC  >  NC_000913/2449192‑2449369

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