Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,547,319 C→A 67.5% R742S (CGT→AGT)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,547,3190CA67.5% 147.5 / 132.6 135R742S (CGT→AGT) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base A (34/57);  total (56/79)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 135 total aligned reads displayed.

CGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAT  >  NC_000913/4547221‑4547411
                                                                                                  |                                                                                            
cGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGGGCGTACCGACTGGCg                                                                                             >  1:2992472/1‑100 (MQ=255)
 gCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGt                                                                                            <  2:845936/100‑1 (MQ=255)
  cTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTg                                                                                           <  1:1581949/100‑1 (MQ=255)
   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTgg                                                                                          <  2:2491283/100‑1 (MQ=255)
   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTg                                                                                           <  1:1165437/99‑1 (MQ=255)
      ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTa                                                                                       <  1:652516/100‑1 (MQ=255)
      ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                       <  2:2895189/100‑1 (MQ=255)
      ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                       <  2:2684330/100‑1 (MQ=255)
         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGc                                                                                    <  2:872480/100‑1 (MQ=255)
         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGc                                                                                    >  2:2902769/1‑100 (MQ=255)
         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGc                                                                                    <  2:3345363/100‑1 (MQ=255)
          gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGcc                                                                                   >  1:1482976/1‑100 (MQ=255)
           ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                  <  1:2512371/100‑1 (MQ=255)
           ccGTTAAACGATACGGTGGTGCTTCTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTCGTTATGCCg                                                                                  <  1:2804015/100‑1 (MQ=255)
            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                 >  1:342856/1‑100 (MQ=255)
              ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                               <  1:7272/100‑1 (MQ=255)
              ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                               <  1:2391936/100‑1 (MQ=255)
                aaaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtgctg                                                                             <  2:1025227/100‑1 (MQ=255)
                 aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                            <  2:279076/100‑1 (MQ=255)
                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGcc                                                                           >  2:1186769/1‑100 (MQ=255)
                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                           <  2:2556611/100‑1 (MQ=255)
                  aCGATACGGTGGTGCTTGTTAAAGCGCATGGCGCAAAAGATGCACCAGTCGAAAACCCGACGGGGGTGCGTCCCGACTGGCGTGGTTATGCCGtgctgca                                                                           <  1:1125843/100‑2 (MQ=255)
                   cGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCt                                                                          >  2:16648/1‑100 (MQ=255)
                     aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                        <  2:2339669/100‑1 (MQ=255)
                        cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                     <  2:1045817/100‑1 (MQ=255)
                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGcc                                                                    <  2:695792/100‑1 (MQ=255)
                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCATTATGcc                                                                    >  1:1058990/1‑100 (MQ=255)
                            ggtgCTTGTTAAAGCGCCTGGCGCAAATGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                 <  1:1295244/100‑1 (MQ=255)
                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCAcn                                                                 <  1:1183396/99‑2 (MQ=255)
                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                <  1:26321/100‑1 (MQ=255)
                              tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                               >  1:388343/1‑100 (MQ=255)
                              tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                               <  2:1805498/100‑1 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                              <  1:58093/100‑1 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                              <  2:2186337/100‑1 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                              <  1:2732940/100‑1 (MQ=255)
                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGCCTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                           >  1:138713/1‑100 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        <  1:2524508/100‑1 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        <  1:3306047/100‑1 (MQ=255)
                                      aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACGAGACGGGGGAGCGTACCGACTGACGGGGTTATGCCGTACTGCCTTATCACACTCAATATCg                                                         >  2:2394843/1‑98 (MQ=255)
                                      aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                       >  2:697970/1‑100 (MQ=255)
                                          cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                   >  1:2022794/1‑100 (MQ=255)
                                             cTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGTCGTGCTGCCTTATGCCACTGAATATCGGGAAAATag                                                <  2:2755617/100‑1 (MQ=255)
                                             cTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                 >  1:1238521/1‑99 (MQ=255)
                                              tGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATaga                                               <  1:2359046/100‑1 (MQ=255)
                                                  gcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAAGATCGGGAAAATAGAGTgg                                           >  1:2518584/1‑100 (MQ=255)
                                                   cAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                          >  1:630943/1‑100 (MQ=255)
                                                   cAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                          <  1:2170340/100‑1 (MQ=255)
                                                    aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcg                                         >  1:3119392/1‑100 (MQ=255)
                                                    aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                          >  1:1314946/1‑99 (MQ=255)
                                                       aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                      <  1:3314932/100‑1 (MQ=255)
                                                       aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                      <  2:1197671/100‑1 (MQ=255)
                                                       aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                      <  2:2646300/100‑1 (MQ=255)
                                                        gATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTgg                                     <  1:3060439/100‑1 (MQ=255)
                                                         aTGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGa                                    <  1:1399707/100‑1 (MQ=255)
                                                          tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                   >  2:3081273/1‑100 (MQ=255)
                                                          tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATCGAGCGGCGCTggg                                    >  1:2234726/1‑98 (MQ=255)
                                                           gCAAAAGTCGAAAGCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTCGAt                                   <  1:4731/99‑1 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                  <  2:401282/100‑1 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                  >  2:126186/1‑100 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                  <  2:342857/100‑1 (MQ=255)
                                                            gaaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTGTGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  1:3217196/2‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATat                                 >  1:1993724/1‑99 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 <  1:2911420/100‑1 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGACTATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  2:609039/1‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  2:2037792/1‑100 (MQ=255)
                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  2:2843062/100‑1 (MQ=255)
                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  2:3008509/100‑1 (MQ=255)
                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  2:3120893/100‑1 (MQ=255)
                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  2:1336839/100‑1 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  1:926421/99‑1 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                               <  2:657708/100‑1 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                               <  2:26257/100‑1 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                               >  1:3190743/1‑100 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTAAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  1:2303312/1‑98 (MQ=255)
                                                               aaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCaa                              <  1:1210337/100‑1 (MQ=255)
                                                               aaGGCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCaa                              <  1:289882/100‑1 (MQ=255)
                                                                 gTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATa                            >  2:2031135/1‑100 (MQ=255)
                                                                    gAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                         <  2:824910/100‑1 (MQ=255)
                                                                    gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                         >  1:1935518/1‑100 (MQ=255)
                                                                    gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                         <  2:1154991/100‑1 (MQ=255)
                                                                       aaCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                      <  2:3135187/100‑1 (MQ=255)
                                                                       aaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                      >  1:1615207/1‑100 (MQ=255)
                                                                        aCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                     <  1:1000371/100‑1 (MQ=255)
                                                                         ccAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                    >  2:2256891/1‑100 (MQ=255)
                                                                         ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCGCCTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                    >  2:806327/1‑100 (MQ=255)
                                                                         ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                      <  1:2826027/98‑1 (MQ=255)
                                                                             aCGGGGGTGCGTACCGACTGGAGTGGTTTTGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                     <  2:595717/95‑1 (MQ=255)
                                                                               gggggTGCGTACCGACTGGAGTGTTTATGCCGTGCTGCCTTATGCCACTGAACATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGAt                 <  1:2722310/97‑1 (MQ=255)
                                                                                  ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc           <  1:523610/100‑1 (MQ=255)
                                                                                   gTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc           <  1:1576480/99‑1 (MQ=255)
                                                                                    tGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGa         >  2:2685679/1‑100 (MQ=255)
                                                                                      cGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAt        >  1:314132/1‑99 (MQ=255)
                                                                                      cGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt       >  2:2934743/1‑100 (MQ=255)
                                                                                      cGTACCGACTGGAGTGGTTATGCAGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAt        >  1:283917/1‑99 (MQ=255)
                                                                                       gTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt      >  2:1593664/1‑100 (MQ=255)
                                                                                       gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt      <  1:723582/100‑1 (MQ=255)
                                                                                       gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt      <  2:1030320/100‑1 (MQ=255)
                                                                                        tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa     >  2:62295/1‑100 (MQ=255)
                                                                                          ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa   <  2:870427/100‑1 (MQ=255)
                                                                                          ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa   >  2:3035565/1‑100 (MQ=255)
                                                                                                  |                                                                                            
CGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAT  >  NC_000913/4547221‑4547411

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