Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 306,943 T→C Y797C (TAC→TGC)  ecpC ← putative fimbrial usher protein EcpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913306,9430TC99.7% 1496.6 / ‑4.8 351Y797C (TAC→TGC) ecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (162/188);  total (163/188)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 351 total aligned reads displayed.

CCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCG  >  NC_000913/306845‑307041
                                                                                                  |                                                                                                  
 ccACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGc                                                                                                 <  2:565246/101‑1 (MQ=255)
 ccACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGc                                                                                                 <  2:565430/101‑1 (MQ=255)
  cACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCg                                                                                                <  1:564952/101‑1 (MQ=255)
  cACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCg                                                                                                <  2:565101/101‑1 (MQ=255)
    cGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGaa                                                                                              <  2:565159/101‑1 (MQ=255)
     gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGaaa                                                                                             >  2:565863/1‑101 (MQ=255)
     gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGaaa                                                                                             <  1:565268/101‑1 (MQ=255)
     gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGaaa                                                                                             <  2:565467/101‑1 (MQ=255)
     gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGaaa                                                                                             >  1:565862/1‑101 (MQ=255)
       tGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAAt                                                                                           <  1:565374/101‑1 (MQ=255)
        gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATc                                                                                          >  2:565866/1‑101 (MQ=255)
        gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATc                                                                                          <  2:565227/101‑1 (MQ=255)
        gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATc                                                                                          <  2:565321/101‑1 (MQ=255)
        gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATc                                                                                          <  1:565323/101‑1 (MQ=255)
         agagGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCg                                                                                         <  1:565157/101‑1 (MQ=255)
           agGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGAt                                                                                       <  2:565562/101‑1 (MQ=255)
           agGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGAt                                                                                       <  1:565304/101‑1 (MQ=255)
            ggCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATa                                                                                      <  2:565536/101‑1 (MQ=255)
             gCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAg                                                                                     >  1:565870/1‑101 (MQ=255)
               cGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTg                                                                                   <  1:565358/101‑1 (MQ=255)
                 cTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTggg                                                                                 <  2:565378/101‑1 (MQ=255)
                  tGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTgggg                                                                                >  1:565876/1‑101 (MQ=255)
                   gCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGt                                                                               <  1:565479/101‑1 (MQ=255)
                    cAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTa                                                                              <  1:565332/101‑1 (MQ=255)
                     aGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTggg                                                                                 >  2:565881/1‑97 (MQ=255)
                     aGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTAt                                                                             <  1:565294/101‑1 (MQ=255)
                      gACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTAtt                                                                            >  1:565882/1‑101 (MQ=255)
                       aCCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTAttt                                                                           <  2:565295/101‑1 (MQ=255)
                        ccACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTc                                                                          <  2:565301/101‑1 (MQ=255)
                          acaTCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCtt                                                                        <  2:565394/101‑1 (MQ=255)
                              cACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATcc                                                                    >  2:565888/1‑101 (MQ=255)
                              cACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGGATTTCTTATcc                                                                    >  1:565889/1‑101 (MQ=255)
                               aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCa                                                                   >  1:565893/1‑101 (MQ=255)
                               aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCa                                                                   <  1:565417/101‑1 (MQ=255)
                               aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCa                                                                   <  1:565618/101‑1 (MQ=255)
                                ccGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCAc                                                                  <  2:565492/101‑1 (MQ=255)
                                    cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTcc                                                              <  1:565342/101‑1 (MQ=255)
                                    cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTcc                                                              <  1:565548/101‑1 (MQ=255)
                                    cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTcc                                                              >  2:565899/1‑101 (MQ=255)
                                      cAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCAt                                                            <  2:565513/101‑1 (MQ=255)
                                       aGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATg                                                           <  2:565284/101‑1 (MQ=255)
                                       aGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATg                                                           >  1:565902/1‑101 (MQ=255)
                                        gACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGa                                                          <  2:565008/101‑1 (MQ=255)
                                         aCGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGAc                                                         >  1:565905/1‑101 (MQ=255)
                                         aCGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGAc                                                         >  1:565909/1‑101 (MQ=255)
                                          cGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACa                                                        >  2:565911/1‑101 (MQ=255)
                                           ggCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACaa                                                       >  1:565913/1‑101 (MQ=255)
                                            gCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACaaa                                                      <  2:565344/101‑1 (MQ=255)
                                               ccGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTc                                                   <  2:565388/101‑1 (MQ=255)
                                               ccGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTc                                                   >  1:565918/1‑101 (MQ=255)
                                               ccGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTc                                                   <  1:565436/101‑1 (MQ=255)
                                                 gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGc                                                 >  2:565921/1‑101 (MQ=255)
                                                 gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGc                                                 <  2:565311/101‑1 (MQ=255)
                                                 gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGc                                                 >  2:565924/1‑101 (MQ=255)
                                                  cgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGcc                                                <  1:565506/101‑1 (MQ=255)
                                                  cgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGcc                                                <  2:565744/101‑1 (MQ=255)
                                                     gcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGtt                                             >  1:565928/1‑101 (MQ=255)
                                                      ccTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCTCTGCAGc                                                                                                 >  2:565931/1‑48 (MQ=255)
                                                        tGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTc                                          <  2:565276/101‑1 (MQ=255)
                                                        tGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTc                                          >  1:565936/1‑101 (MQ=255)
                                                          gTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCAt                                        <  2:565403/101‑1 (MQ=255)
                                                          gTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCAt                                        >  1:565939/1‑101 (MQ=255)
                                                           ttGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATc                                       <  2:565444/101‑1 (MQ=255)
                                                            tGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCg                                      >  2:565944/1‑101 (MQ=255)
                                                            tGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCg                                      <  2:565460/101‑1 (MQ=255)
                                                             gAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCgg                                     >  1:565946/1‑101 (MQ=255)
                                                              aGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggt                                    >  1:565950/1‑101 (MQ=255)
                                                               gTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggtt                                   >  2:565953/1‑101 (MQ=255)
                                                                ttCCAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttc                                  <  2:565349/101‑1 (MQ=255)
                                                                   cAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcggg                               <  1:565177/101‑1 (MQ=255)
                                                                   cAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcggg                               >  1:565958/1‑101 (MQ=255)
                                                                   cAGAGCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCg                                      >  2:565960/1‑94 (MQ=255)
                                                                    agagCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggt                              <  2:565662/101‑1 (MQ=255)
                                                                      agCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggttc                            >  1:565967/1‑101 (MQ=255)
                                                                      agCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggttc                            >  2:565965/1‑101 (MQ=255)
                                                                       gCCACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggttcg                           <  2:565524/101‑1 (MQ=255)
                                                                         cACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGc                         >  1:565972/1‑101 (MQ=255)
                                                                         cACTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGc                         >  2:565974/1‑101 (MQ=255)
                                                                          acctCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCAt                                                            <  2:565039/62‑1 (MQ=255)
                                                                          aCTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGcc                        <  1:565464/101‑1 (MQ=255)
                                                                          aCTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGcc                        <  1:565278/101‑1 (MQ=255)
                                                                          aCTTCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGcc                        >  1:565977/1‑101 (MQ=255)
                                                                            ttCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGa                      <  1:565016/101‑1 (MQ=255)
                                                                            ttCGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGa                      >  2:565983/1‑101 (MQ=255)
                                                                              cGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGata                    >  2:565986/1‑101 (MQ=255)
                                                                              cGCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGGATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGata                    >  1:565987/1‑101 (MQ=255)
                                                                               gCAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat                   <  1:565245/101‑1 (MQ=255)
                                                                                cAGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATg                  >  1:565992/1‑101 (MQ=255)
                                                                                 aGGTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATgg                 >  1:565996/1‑101 (MQ=255)
                                                                                  ggTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGgt                >  2:565999/1‑101 (MQ=255)
                                                                                  ggTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGgt                <  1:565685/101‑1 (MQ=255)
                                                                                  ggTTTTATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGgt                <  1:565452/101‑1 (MQ=255)
                                                                                      ttATTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGgttgtt            >  1:566001/1‑101 (MQ=255)
                                                                                        aTTGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTaa          >  1:566007/1‑101 (MQ=255)
                                                                                         ttGCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAAt         >  2:566009/1‑101 (MQ=255)
                                                                                           gCCACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATcc       <  2:565498/101‑1 (MQ=255)
                                                                                            ccACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCg      >  1:566012/1‑101 (MQ=255)
                                                                                             cACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGt     >  2:566019/1‑101 (MQ=255)
                                                                                             cACTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGt     >  1:566018/1‑101 (MQ=255)
                                                                                               cTGCAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGc   <  1:565117/101‑1 (MQ=255)
                                                                                                  |                                                                                                  
CCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCG  >  NC_000913/306845‑307041

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