Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 305,808 T→C Q337R (CAG→CGG)  ecpD ← fimbrial adhesin EcpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913305,8080TC97.7% 309.1 / ‑2.4 87Q337R (CAG→CGG) ecpDfimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (43/42);  total (44/43)
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.37e-01

GAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGC  >  NC_000913/305672‑305941
                                                                                                                                        |                                                                                                                                     
gAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCt                                                                                                                                    >  1:533566/1‑140 (MQ=255)
   aTCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTaa                                                                                                                                 >  1:849122/1‑140 (MQ=255)
      aCGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAaaga                                                                                                                              >  1:341672/1‑140 (MQ=255)
               gATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCg                                                                                                                     >  1:204306/1‑140 (MQ=255)
                   aTCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaa                                                                                                                 <  2:760317/140‑1 (MQ=255)
                    tCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                >  1:642147/1‑140 (MQ=255)
                    tCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                >  2:483256/1‑140 (MQ=255)
                    tCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                <  2:259434/140‑1 (MQ=255)
                    tCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                <  1:602230/140‑1 (MQ=255)
                     cAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaaga                                                                                                               <  1:769696/140‑1 (MQ=255)
                      aGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAATGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                >  2:470435/1‑138 (MQ=255)
                      aGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAATGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGagaag                                                                                                                <  1:470435/138‑1 (MQ=255)
                      aGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAgaagaa                                                                                                              <  2:204306/140‑1 (MQ=255)
                      aGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAgaagaa                                                                                                              >  2:105458/1‑140 (MQ=255)
                      aGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAgaagaa                                                                                                              >  2:477538/1‑140 (MQ=255)
                       gCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGaaa                                                                                                             >  2:506361/1‑140 (MQ=255)
                        cGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGaaaa                                                                                                            >  1:840317/1‑140 (MQ=255)
                        cGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGaaaa                                                                                                            <  1:105458/140‑1 (MQ=255)
                         gtgtTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTc                                                                                                                      >  2:737296/1‑129 (MQ=255)
                         gtgtTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTc                                                                                                                      <  1:737296/129‑1 (MQ=255)
                          tgtTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACa                                                                                                          <  1:483256/140‑1 (MQ=255)
                          tgtTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACa                                                                                                          <  2:840317/140‑1 (MQ=255)
                           gtTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACAt                                                                                                         >  1:538529/1‑140 (MQ=255)
                               gAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGAt                                                                                                     <  1:506361/140‑1 (MQ=255)
                               gAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGAt                                                                                                     <  1:455498/140‑1 (MQ=255)
                                    aGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGt                                                                                                >  1:36357/1‑140 (MQ=255)
                                    aGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCACTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGt                                                                                                <  2:285682/140‑1 (MQ=255)
                                     gAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTc                                                                                               >  2:166174/1‑140 (MQ=255)
                                      aaaaTTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCt                                                                                              >  1:673233/1‑140 (MQ=255)
                                      aaaaTTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCt                                                                                              >  1:299022/1‑140 (MQ=255)
                                             aTACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGctct                                                                                       >  1:794204/1‑140 (MQ=255)
                                               aCCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTc                                                                                     <  1:454011/140‑1 (MQ=255)
                                                 cagccagACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCAt                                                                                   <  2:319159/140‑1 (MQ=255)
                                                   gccagACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCAtct                                                                                 >  2:70888/1‑140 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  1:51289/140‑1 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  2:36357/140‑1 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  2:341672/140‑1 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  1:70888/140‑1 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  2:82577/140‑1 (MQ=255)
                                                         cTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCACTGATACCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  2:800823/140‑1 (MQ=255)
                                                          tCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGaa                                                                          <  2:538529/140‑1 (MQ=255)
                                                           ccGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAg                                                                         <  2:614805/140‑1 (MQ=255)
                                                            cGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAgg                                                                        <  1:56444/140‑1 (MQ=255)
                                                            cGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAgg                                                                        >  2:463647/1‑140 (MQ=255)
                                                            cGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAgg                                                                        <  2:642147/140‑1 (MQ=255)
                                                                  tgtgGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGGAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAg                                                                         <  2:120541/133‑1 (MQ=255)
                                                                  tgtgGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGGAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAg                                                                         >  1:120541/1‑133 (MQ=255)
                                                                      gtgtTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATa                                                              >  2:134182/1‑140 (MQ=255)
                                                                        gtTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           >  2:95830/1‑125 (MQ=255)
                                                                        gtTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGa                                                                           <  1:95830/125‑1 (MQ=255)
                                                                        gtTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAgg                                                            <  2:794204/140‑1 (MQ=255)
                                                                         ttCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGgtg                                                                      <  2:669151/129‑1 (MQ=255)
                                                                         ttCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGgtg                                                                      >  1:669151/1‑129 (MQ=255)
                                                                            tGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGt                                                        >  2:66275/1‑140 (MQ=255)
                                                                             gAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTa                                                       <  2:123505/140‑1 (MQ=255)
                                                                              aaaGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTAt                                                      >  1:497033/1‑140 (MQ=255)
                                                                                 gCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGc                                                   >  1:460061/1‑140 (MQ=255)
                                                                                     aaattGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCAc                                               >  1:686921/4‑140 (MQ=255)
                                                                                              aGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCa                                      <  2:686921/140‑1 (MQ=255)
                                                                                                aTCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATg                                    <  2:108789/140‑1 (MQ=255)
                                                                                                aTCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATg                                    <  2:299022/140‑1 (MQ=255)
                                                                                                aTCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATg                                    <  2:747909/140‑1 (MQ=255)
                                                                                                  cTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGaa                                  >  2:316078/1‑140 (MQ=255)
                                                                                                     tGGTGTTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTg                               >  1:827352/1‑140 (MQ=255)
                                                                                                        tgtTGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCa                            >  1:476230/1‑140 (MQ=255)
                                                                                                          ttGATATCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCAtt                          >  2:846486/1‑140 (MQ=255)
                                                                                                             atatCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAg                       >  2:414509/1‑140 (MQ=255)
                                                                                                              tatCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGa                      <  1:66275/140‑1 (MQ=255)
                                                                                                              tatCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGa                      <  2:673233/140‑1 (MQ=255)
                                                                                                               atCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGAc                     >  1:832776/1‑140 (MQ=255)
                                                                                                                tCGCTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACt                    <  2:827352/140‑1 (MQ=255)
                                                                                                                   cTGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGaa                 <  2:128700/140‑1 (MQ=255)
                                                                                                                    tGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACt                                >  2:22264/1‑124 (MQ=255)
                                                                                                                    tGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACt                                <  1:22264/124‑1 (MQ=255)
                                                                                                                    tGATCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAAc                >  1:628440/1‑140 (MQ=255)
                                                                                                                       tCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAACCATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGc             <  2:567837/140‑1 (MQ=255)
                                                                                                                       tCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGc             >  1:567837/1‑140 (MQ=255)
                                                                                                                       tCCCGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGc             >  2:432919/1‑140 (MQ=255)
                                                                                                                        cccGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACt                    >  1:595811/1‑132 (MQ=255)
                                                                                                                        cccGAGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCa            >  1:272892/1‑140 (MQ=255)
                                                                                                                        cccGAGGTTCACCATCCGCTTAAAGAACTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACt                    <  2:595811/132‑1 (MQ=255)
                                                                                                                            aGGTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATc        >  1:251596/1‑140 (MQ=255)
                                                                                                                              gTTCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGt      >  1:857918/1‑140 (MQ=255)
                                                                                                                               ttCACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTa     <  1:316078/140‑1 (MQ=255)
                                                                                                                                 cACCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGCGTAGTAATAGGCGGAATTGCTTACCGGTTTCCATGAACTGCCATTCGGACTGCGATGCATATCGTAc  >  2:68828/1‑139 (MQ=255)
                                                                                                                                  aCCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGc  <  2:53306/140‑1 (MQ=255)
                                                                                                                                  aCCATCCGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGc  >  1:731223/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGC  >  NC_000913/305672‑305941

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