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,2800TC99.6% 713.3 / ‑5.1 223T109A (ACC→GCC) b2335putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (98/124);  total (98/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 228 total aligned reads displayed.

GTTATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGC  >  NC_000913/2449191‑2449369
                                                                                         |                                                                                         
gTTATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGat                                                                                        >  1:1646937/1‑93 (MQ=255)
 ttATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGata                                                                                       >  2:2691423/1‑93 (MQ=255)
  tATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGatat                                                                                      <  1:1524531/93‑1 (MQ=255)
  tATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGatat                                                                                      >  2:4528964/1‑93 (MQ=255)
   aTTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATAtt                                                                                     >  2:650664/1‑93 (MQ=255)
   aTTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATAtt                                                                                     <  1:9169/93‑1 (MQ=255)
   aTTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATAtt                                                                                     >  1:2400919/1‑93 (MQ=255)
    ttGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTg                                                                                    >  1:5789471/1‑93 (MQ=255)
     tGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTgg                                                                                   <  1:1177648/93‑1 (MQ=255)
     tGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTgg                                                                                   <  2:759758/93‑1 (MQ=255)
      gATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGt                                                                                  >  1:3210326/1‑93 (MQ=255)
      gATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGt                                                                                  <  2:3511904/93‑1 (MQ=255)
      gATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGt                                                                                  >  2:5751901/1‑93 (MQ=255)
       atatTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTa                                                                                 <  2:6260670/93‑1 (MQ=255)
       atatTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTa                                                                                 >  1:4327862/1‑93 (MQ=255)
        tatTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTAc                                                                                >  1:1627720/1‑93 (MQ=255)
          ttAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTg                                                                              <  2:1742693/93‑1 (MQ=255)
          ttAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTg                                                                              <  2:3702513/93‑1 (MQ=255)
            aaaTGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAg                                                                            <  2:93197/93‑1 (MQ=255)
              aTGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTa                                                                          <  2:3366202/93‑1 (MQ=255)
              aTGGCGCCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTCGTACTGAgga                                                                          >  1:5727564/1‑91 (MQ=255)
                ggCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAcc                                                                        >  2:5261795/1‑93 (MQ=255)
                ggCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAcc                                                                        <  1:5054343/93‑1 (MQ=255)
                ggCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTAcc                                                                        >  1:4790737/1‑93 (MQ=255)
                  cGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGt                                                                      >  2:637603/1‑93 (MQ=255)
                   gTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGtt                                                                     >  2:2300774/1‑93 (MQ=255)
                    tCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGttt                                                                    >  2:1222541/1‑93 (MQ=255)
                     ccAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTc                                                                   <  2:4068459/93‑1 (MQ=255)
                     ccAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTc                                                                   <  1:1071462/93‑1 (MQ=255)
                     ccAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTc                                                                   >  1:3277659/1‑93 (MQ=255)
                       aGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGc                                                                 >  2:5683623/1‑93 (MQ=255)
                        gCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGcc                                                                >  1:2187144/1‑93 (MQ=255)
                          tATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGt                                                              <  2:2077581/93‑1 (MQ=255)
                          tATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGt                                                              <  2:4179130/93‑1 (MQ=255)
                            ttttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTa                                                            >  1:1299857/1‑93 (MQ=255)
                            ttttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTa                                                            >  1:6462333/1‑93 (MQ=255)
                              ttCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAAt                                                          >  1:5683059/1‑93 (MQ=255)
                               tCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATc                                                         >  1:5076578/1‑93 (MQ=255)
                                cACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTACTCg                                                        >  1:2044200/1‑93 (MQ=255)
                                  ccGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGtt                                                      <  2:1321051/93‑1 (MQ=255)
                                  ccGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGtt                                                      >  2:2763891/1‑93 (MQ=255)
                                  ccGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGtt                                                      >  1:450275/1‑93 (MQ=255)
                                   cGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTg                                                     <  1:2646353/93‑1 (MQ=255)
                                    gACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGttaa                                                    >  2:5455978/1‑91 (MQ=255)
                                    gACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGa                                                    >  1:2531526/1‑93 (MQ=255)
                                       ggCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGt                                                 >  1:1969525/1‑93 (MQ=255)
                                         cAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAt                                               <  2:3228772/93‑1 (MQ=255)
                                          aaaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTAtt                                              >  1:4487864/1‑93 (MQ=255)
                                           aaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTg                                             <  1:4023082/93‑1 (MQ=255)
                                           aaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTg                                             >  1:759025/1‑93 (MQ=255)
                                           aaCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTg                                             >  2:1114922/1‑93 (MQ=255)
                                            aCAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGc                                            >  1:123221/1‑93 (MQ=255)
                                              aaTGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGccc                                          <  1:1017910/93‑1 (MQ=255)
                                               aTGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCt                                         <  1:3858060/93‑1 (MQ=255)
                                               aTGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGCACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCt                                         >  1:2768942/1‑93 (MQ=255)
                                                tGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTCATCGTTGAGGTATTGCCCtt                                        >  1:1902079/1‑93 (MQ=255)
                                                tGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCtt                                        >  1:3610919/1‑93 (MQ=255)
                                                tGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCtt                                        >  2:1851388/1‑93 (MQ=255)
                                                 gAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTa                                       <  1:338379/93‑1 (MQ=255)
                                                  aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGTTATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  1:5857449/93‑1 (MQ=255)
                                                  aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  1:5891078/93‑1 (MQ=255)
                                                  aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  2:3095967/93‑1 (MQ=255)
                                                  aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  2:2850116/93‑1 (MQ=255)
                                                  aGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTaa                                      <  1:542541/93‑1 (MQ=255)
                                                    ttATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGt                                    >  2:5329999/1‑93 (MQ=255)
                                                     tATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGtt                                   <  2:5248869/93‑1 (MQ=255)
                                                      aTCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTg                                  <  2:676848/93‑1 (MQ=255)
                                                       tCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGc                                 <  2:198815/93‑1 (MQ=255)
                                                         ggCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCAt                               >  1:6236563/1‑93 (MQ=255)
                                                         ggCGCTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCAt                               <  1:4647772/93‑1 (MQ=255)
                                                           cgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATcc                             >  2:5813470/1‑93 (MQ=255)
                                                           cgcTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATcc                             <  2:6080184/93‑1 (MQ=255)
                                                             cTATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGCTTCGCCGTTAATCGCTGAGGCCTTGCCCTTAAGTTGCATCCGc                           >  2:1594728/1‑93 (MQ=255)
                                                              tATTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCa                          >  1:5178062/1‑93 (MQ=255)
                                                               aTTTTCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCaa                         <  1:3773578/93‑1 (MQ=255)
                                                                ttttCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAg                        >  2:1437965/1‑93 (MQ=255)
                                                                 tttCTATCCGAATGCCCAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGt                       <  2:1452619/93‑1 (MQ=255)
                                                                 tttCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGt                       <  1:3475223/93‑1 (MQ=255)
                                                                 tttCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGt                       <  2:5535262/93‑1 (MQ=255)
                                                                 tttCTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGt                       <  2:244912/93‑1 (MQ=255)
                                                                    cTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGt                    <  1:4182306/93‑1 (MQ=255)
                                                                    cTATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGt                    <  2:5166008/93‑1 (MQ=255)
                                                                     tATCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGtt                   <  1:412891/93‑1 (MQ=255)
                                                                      aTCCGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTa                  >  2:379432/1‑93 (MQ=255)
                                                                        ccGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGc                <  2:4270679/93‑1 (MQ=255)
                                                                         cGAATGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCa               >  2:5098143/1‑93 (MQ=255)
                                                                           aaTGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAg             <  2:5852014/93‑1 (MQ=255)
                                                                            aTGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGa            >  1:554311/1‑93 (MQ=255)
                                                                            aTGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGa            <  2:2130074/93‑1 (MQ=255)
                                                                            aTGCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGa            <  2:4211104/93‑1 (MQ=255)
                                                                              gCCCAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACt          <  1:6038690/93‑1 (MQ=255)
                                                                               cccAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTg         <  2:587933/93‑1 (MQ=255)
                                                                                 cAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCTTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:769996/93‑1 (MQ=255)
                                                                                 cAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       <  1:2855771/93‑1 (MQ=255)
                                                                                 cAGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGtt       >  2:1555987/1‑93 (MQ=255)
                                                                                  aGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGttt      <  2:3657421/93‑1 (MQ=255)
                                                                                  aGCCCGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGttt      <  2:4733160/93‑1 (MQ=255)
                                                                                    cccGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtgt    <  2:1076409/93‑1 (MQ=255)
                                                                                     ccGGCGATATTGGTACTGAGTACCGTTTCGCGGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTtgtg   <  1:6020276/93‑1 (MQ=255)
                                                                                      cGGCGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGc  <  1:2901406/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTTATTGATATTAAATGGCGTCCAGCTATTTTCACCGACGGCAAACAATGAGTTATCGGCGCTATTTTCTATCCGAATGCCCAGCCCGGTGATATTGGTACTGAGTACCGTTTCGCCGTTAATCGTTGAGGTATTGCCCTTAAGTTGCATCCGCAAGTCGTTAGCAAGACTGTTTGTGC  >  NC_000913/2449191‑2449369

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