Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,541,157 C→T R159C (CGT→TGT)  fimA → major type 1 subunit fimbrin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,541,1570CT96.2% 702.1 / 6.0 211R159C (CGT→TGT) fimAmajor type 1 subunit fimbrin
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (100/103);  total (105/106)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 222 total aligned reads displayed.

AGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG  >  NC_000913/4541067‑4541246
                                                                                          |                                                                                         
aGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCgtgt                                                                                         >  1:474409/1‑93 (MQ=255)
aGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCgtgt                                                                                         >  2:646577/1‑93 (MQ=255)
aGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCgtgt                                                                                         >  1:203899/1‑93 (MQ=255)
 gAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGtt                                                                                        >  2:798059/1‑93 (MQ=255)
 gAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGGTCAGAAACAACCCTGAATAACTGAACCAATACCATTCCGTTCCAGTCGTGtt                                                                                        >  1:2406360/1‑93 (MQ=255)
  aaCGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTa                                                                                       <  2:3090490/93‑1 (MQ=255)
    cGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTAtt                                                                                     <  1:1280537/93‑1 (MQ=255)
     gggTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTAttt                                                                                    <  2:3943592/93‑1 (MQ=255)
      ggTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTAtttt                                                                                   <  1:221163/93‑1 (MQ=255)
          ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCaa                                                                               >  2:2236385/1‑93 (MQ=255)
           tgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAAc                                                                              >  1:1222898/1‑93 (MQ=255)
            gcgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAAcc                                                                             <  2:3458581/93‑1 (MQ=255)
             cgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCg                                                                            >  2:2942094/1‑93 (MQ=255)
              gcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCgg                                                                           <  1:3223277/93‑1 (MQ=255)
              gcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCgg                                                                           <  1:1459904/93‑1 (MQ=255)
               cTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCggg                                                                          >  1:1175416/1‑93 (MQ=255)
                tGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCgggg                                                                         <  2:4535388/93‑1 (MQ=255)
                 gACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGc                                                                        <  2:4260893/93‑1 (MQ=255)
                 gACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGc                                                                        >  1:1845609/1‑93 (MQ=255)
                   cGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCg                                                                      <  2:1494104/93‑1 (MQ=255)
                     cTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                    >  1:2418038/1‑93 (MQ=255)
                      tGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATGCCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCaa                                                                   <  1:1356463/93‑1 (MQ=255)
                      tGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCaa                                                                   >  2:1351045/1‑93 (MQ=255)
                      tGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCaa                                                                   >  1:876852/1‑93 (MQ=255)
                       ggATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAc                                                                  >  2:2866211/1‑93 (MQ=255)
                       ggATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAc                                                                  >  2:4661178/1‑93 (MQ=255)
                        gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAcc                                                                 >  1:3405030/1‑93 (MQ=255)
                        gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAcc                                                                 <  1:210305/93‑1 (MQ=255)
                        gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                 <  2:2904080/93‑1 (MQ=255)
                          tGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAcccc                                                               <  1:508003/93‑1 (MQ=255)
                              gCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                                                           <  2:1274119/93‑1 (MQ=255)
                                gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGtgc                                                         <  1:991204/93‑1 (MQ=255)
                                    tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTgctgct                                                     >  1:1823014/1‑93 (MQ=255)
                                        gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCAGGGTGCTGTTAATg                                                 >  2:3412874/1‑93 (MQ=255)
                                         ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                <  1:2918980/93‑1 (MQ=255)
                                         ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                >  2:1079370/1‑93 (MQ=255)
                                           cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCgg                                              >  2:96585/1‑93 (MQ=255)
                                             gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGAt                                            >  2:564519/1‑93 (MQ=255)
                                               aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                          <  2:419935/93‑1 (MQ=255)
                                               aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                          <  1:1674909/93‑1 (MQ=255)
                                                acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                         >  2:211873/1‑93 (MQ=255)
                                                  aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                       >  1:168687/1‑93 (MQ=255)
                                                   acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                      <  1:2805865/93‑1 (MQ=255)
                                                   acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                      >  2:1604807/1‑93 (MQ=255)
                                                     ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                    >  1:3004859/1‑93 (MQ=255)
                                                      cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGTTACCGGGGCCGCCACCCCGGGTGCTGCTAATGCGGATGCGACCGTc                                   >  2:2441130/1‑93 (MQ=255)
                                                      cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtta                                   >  2:666817/1‑92 (MQ=255)
                                                      cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                   <  1:2951662/93‑1 (MQ=255)
                                                       tGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                  <  1:3195187/93‑1 (MQ=255)
                                                       tGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                  >  1:3731055/1‑93 (MQ=255)
                                                       tGAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGCGCTGCTAATGCGGATGCGACCTTCa                                  >  2:3688796/1‑93 (MQ=255)
                                                        gAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                 >  2:2485177/1‑93 (MQ=255)
                                                        gAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                 <  1:1511798/93‑1 (MQ=255)
                                                        gAATAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                 <  1:2529009/93‑1 (MQ=255)
                                                         aaTAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                <  2:1558963/93‑1 (MQ=255)
                                                         aaTAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                <  1:4605676/93‑1 (MQ=255)
                                                          aTAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                               <  1:4354056/93‑1 (MQ=255)
                                                           tAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                              >  1:1602326/1‑93 (MQ=255)
                                                           tAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                              <  2:1152996/93‑1 (MQ=255)
                                                           tAACGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCAGTGCCGCAACCCCGTGGGCTGCTAATGCGGATGCGAGCTTCCAGGt                              >  1:58253/1‑93 (MQ=255)
                                                             aCGGAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                            <  2:2255302/93‑1 (MQ=255)
                                                              cGGAACCAATACCATGCCGTTCCAGGCGTGGTAGTTTGCAACCGGGGCCGCAGCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                           <  2:4118170/93‑1 (MQ=255)
                                                               ggAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                          <  1:3521922/93‑1 (MQ=255)
                                                                gAACCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                         >  1:1404387/1‑93 (MQ=255)
                                                                 aaCCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                        <  2:477662/93‑1 (MQ=255)
                                                                 aaCCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                        <  1:1122288/93‑1 (MQ=255)
                                                                  aCCAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                       >  1:4192616/1‑93 (MQ=255)
                                                                   ccAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                      <  2:872026/93‑1 (MQ=255)
                                                                   ccAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                      >  2:4367019/1‑93 (MQ=255)
                                                                   ccAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                      >  1:4155655/1‑93 (MQ=255)
                                                                    cAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                     <  2:1844754/93‑1 (MQ=255)
                                                                    cAATACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                     >  2:4067551/1‑93 (MQ=255)
                                                                     aaTACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                    <  1:3561342/93‑1 (MQ=255)
                                                                     aaTACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                    >  1:4115385/1‑93 (MQ=255)
                                                                      aTACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                   >  1:1011903/1‑93 (MQ=255)
                                                                        aCCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                 <  1:1588415/93‑1 (MQ=255)
                                                                        aCCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                 <  2:504781/93‑1 (MQ=255)
                                                                         ccATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                >  2:2062065/1‑93 (MQ=255)
                                                                         ccATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGAATCAATAAc                >  1:3988968/1‑93 (MQ=255)
                                                                            ttCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCtt             <  1:4424432/93‑2 (MQ=255)
                                                                            ttCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa             <  1:4073810/93‑1 (MQ=255)
                                                                            ttCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa             >  1:4204871/1‑93 (MQ=255)
                                                                            ttCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa             >  1:3891982/1‑93 (MQ=255)
                                                                             tCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCACCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc            <  2:4004479/93‑1 (MQ=255)
                                                                             tCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc            <  2:515088/93‑1 (MQ=255)
                                                                             tCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc            <  2:4547578/93‑1 (MQ=255)
                                                                               cGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc          <  2:2022985/93‑1 (MQ=255)
                                                                                gTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa         >  1:4286857/1‑93 (MQ=255)
                                                                                gTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa         <  1:2649775/93‑1 (MQ=255)
                                                                                gTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa         >  1:1452369/1‑93 (MQ=255)
                                                                                  tCCAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg       <  1:556754/93‑1 (MQ=255)
                                                                                    cAGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt     <  2:2670810/93‑1 (MQ=255)
                                                                                     aGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc    >  1:2743315/1‑93 (MQ=255)
                                                                                     aGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc    >  2:2777228/1‑93 (MQ=255)
                                                                                     aGGCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc    >  2:878589/1‑93 (MQ=255)
                                                                                      ggCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa   >  2:2956038/1‑93 (MQ=255)
                                                                                      ggCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa   >  1:2473804/1‑93 (MQ=255)
                                                                                      ggCGTGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa   <  2:3841737/93‑1 (MQ=255)
                                                                                       gcgtGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg  >  1:3941908/1‑93 (MQ=255)
                                                                                       gcgtGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg  >  1:3069243/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG  >  NC_000913/4541067‑4541246

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