Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,029 Δ1 bp 6.1% coding (557/1557 nt) yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,0290A. 6.1% 482.5 / 22.4 131coding (557/1557 nt)yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (64/59);  new base . (4/4);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 141 total aligned reads displayed.

GTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGC  >  NC_000913/1900884‑1901177
                                                                                                                                                 |                                                                                                                                                     
gTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccacca                                                                                                                                                    >  2:437379/1‑149 (MQ=255)
     gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                               >  2:563061/1‑149 (MQ=255)
     gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                               >  1:406987/1‑149 (MQ=255)
      aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                              <  1:576583/149‑1 (MQ=255)
       cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTggg                                                                                                                                            <  1:563061/149‑1 (MQ=255)
        gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgggg                                                                                                                                           >  2:968963/1‑149 (MQ=255)
        gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgggg                                                                                                                                           <  2:406987/149‑1 (MQ=255)
                   ttAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                >  2:1301852/1‑149 (MQ=255)
                   ttAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCt                                                                                                                                       <  3:454589/142‑1 (MQ=38)
                   ttAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCt                                                                                                                                       <  1:1279591/142‑1 (MQ=37)
                     aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGt                                                                                                                                    >  2:137474/1‑143 (MQ=255)
                      aaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCATCACCGTGGGGTGCTGGTTAACGaa                                                                                                                             >  3:1097824/1‑149 (MQ=255)
                       acacTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            >  1:860599/1‑149 (MQ=255)
                         acTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcg                                                                                                                          >  3:300325/1‑149 (MQ=255)
                          cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgc                                                                                                                         <  3:120853/149‑1 (MQ=255)
                          cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgc                                                                                                                         <  1:835993/149‑1 (MQ=255)
                           ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                >  2:655986/1‑141 (MQ=255)
                           ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                <  1:655986/141‑1 (MQ=255)
                           ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgcg                                                                                                                        >  4:1430367/1‑149 (MQ=255)
                            tCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGt                                                                                                                       <  4:300325/149‑1 (MQ=255)
                              gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCa                                                                                                                     >  2:902052/1‑149 (MQ=255)
                               atgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                >  4:7870/1‑137 (MQ=255)
                               atgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                <  3:7870/137‑1 (MQ=255)
                                 gatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGt                                                                                                                                    <  4:984655/131‑1 (MQ=255)
                                 gatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGt                                                                                                                                    >  3:984655/1‑131 (MQ=255)
                                 gatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCg                                                                                                                  >  1:1106389/1‑149 (MQ=255)
                                    gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                               >  1:1178003/1‑149 (MQ=255)
                                    gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                               <  1:968963/149‑1 (MQ=255)
                                     atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACTACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  4:607190/149‑1 (MQ=255)
                                     atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            >  3:913199/1‑135 (MQ=255)
                                     atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            <  4:913199/135‑1 (MQ=255)
                                     atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  1:437379/149‑1 (MQ=255)
                                     atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  3:1430367/149‑1 (MQ=255)
                                        gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           <  1:902052/149‑1 (MQ=255)
                                        gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctgtc                                                                                                                                  <  2:775117/126‑6 (MQ=37)
                                        gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctgtc                                                                                                                                  >  1:775117/1‑121 (MQ=37)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                              >  4:1073230/1‑110 (MQ=255)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                      <  1:339656/117‑1 (MQ=255)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                      >  2:339656/1‑117 (MQ=255)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            <  1:127025/127‑1 (MQ=255)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            >  2:127025/1‑127 (MQ=255)
                                             gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                      <  3:933470/149‑1 (MQ=255)
                                              ccAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTaa                                                                                                                                 >  2:683182/1‑121 (MQ=255)
                                                aaTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTgg                                                                                                                                     <  2:132922/115‑1 (MQ=255)
                                                aaTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTgg                                                                                                                                     >  1:132922/1‑115 (MQ=255)
                                                  tCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATc                                                                                                                           >  4:498283/1‑123 (MQ=255)
                                                  tCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCACACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATc                                                                                                                           <  3:498283/123‑1 (MQ=255)
                                                       gCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCATCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAc                                                                                            <  4:1097824/149‑1 (MQ=255)
                                                        cATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                           <  2:1106389/149‑1 (MQ=255)
                                                         aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCg                                                                                          >  3:233120/1‑149 (MQ=255)
                                                         aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtc                                                                                                                                      <  2:210652/105‑2 (MQ=25)
                                                         aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtc                                                                                                                                      >  1:210652/1‑104 (MQ=37)
                                                                   ccTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAAt                                                                                                                            >  1:335471/1‑105 (MQ=255)
                                                                     tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 <  1:1309385/146‑1 (MQ=255)
                                                                       tCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 <  3:815415/144‑1 (MQ=255)
                                                                       tCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 >  4:815415/1‑144 (MQ=255)
                                                                        cGGAATGTGGAAACCGAAACCTGCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTa                                                                           <  4:214616/149‑1 (MQ=255)
                                                                         ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           <  4:374993/116‑1 (MQ=255)
                                                                         ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           >  3:374993/1‑116 (MQ=255)
                                                                                 gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           >  2:910956/1‑108 (MQ=255)
                                                                                 gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCca                                                                  >  3:203152/1‑149 (MQ=255)
                                                                                            cTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTg                                                                                                             >  4:794627/1‑95 (MQ=255)
                                                                                            cTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTg                                                                                                             <  3:794627/95‑1 (MQ=255)
                                                                                            cTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGA‑‑ACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAg                                                     >  2:1123716/1‑149 (MQ=255)
                                                                                             ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGAt                                                                                                          >  1:639775/1‑97 (MQ=255)
                                                                                             ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                      >  1:388149/1‑149 (MQ=255)
                                                                                                   caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                      <  2:1277156/63‑1 (MQ=255)
                                                                                                   caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGt                                                >  4:446188/1‑149 (MQ=255)
                                                                                                    accaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                               >  2:871638/1‑149 (MQ=255)
                                                                                                     ccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGAt                                                                                                          <  3:318467/89‑1 (MQ=255)
                                                                                                      caGACTCAGACCAATCATTAACAGGAAGCTCTGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                >  4:900146/1‑66 (MQ=255)
                                                                                                      caGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                <  3:900146/66‑1 (MQ=255)
                                                                                                       aGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAAc                                            <  1:1094777/149‑1 (MQ=255)
                                                                                                        gACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAAcc                                           <  4:233120/149‑1 (MQ=255)
                                                                                                           tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                      >  1:716483/1‑55 (MQ=255)
                                                                                                           tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                               <  1:201937/142‑1 (MQ=255)
                                                                                                              gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCt                                     >  2:1283873/1‑149 (MQ=255)
                                                                                                                  aaTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTg                                 >  4:989101/1‑149 (MQ=255)
                                                                                                                      aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATc                                                                                     >  1:1060975/1‑93 (MQ=255)
                                                                                                                      aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATc                                                                                     <  2:1060975/93‑1 (MQ=255)
                                                                                                                       ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 >  3:322820/1‑96 (MQ=255)
                                                                                                                       ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 <  4:322820/96‑1 (MQ=255)
                                                                                                                        tAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcat                                                                 <  4:31040/111‑1 (MQ=255)
                                                                                                                         aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATc                                                                                     <  2:757858/90‑1 (MQ=255)
                                                                                                                         aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAt                          <  4:203152/149‑1 (MQ=255)
                                                                                                                         aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAt                          >  3:1450286/1‑149 (MQ=255)
                                                                                                                          aCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGc                                                       >  2:731299/1‑119 (MQ=255)
                                                                                                                          aCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGc                                                       <  1:731299/119‑1 (MQ=255)
                                                                                                                              gAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCg                     >  2:139898/1‑149 (MQ=255)
                                                                                                                                   tCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAt                                                   <  4:1220334/114‑1 (MQ=255)
                                                                                                                                   tCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAt                                                   >  3:1220334/1‑114 (MQ=255)
                                                                                                                                   tCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc                >  3:1148967/1‑149 (MQ=255)
                                                                                                                                     agacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                                <  2:241554/83‑1 (MQ=255)
                                                                                                                                     agacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                                >  1:241554/1‑83 (MQ=255)
                                                                                                                                      gacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                               <  1:902770/115‑1 (MQ=255)
                                                                                                                                      gacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCTCATCCAAc             >  2:1187843/1‑149 (MQ=255)
                                                                                                                                                caccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                      >  1:9716/1‑98 (MQ=255)
                                                                                                                                                caccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                      <  2:9716/98‑1 (MQ=255)
                                                                                                                                                caccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTAc                                   <  2:604198/117‑1 (MQ=255)
                                                                                                                                                caccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACgg   <  3:994401/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGC  >  NC_000913/1900884‑1901177

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