Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,385,208 T→A 5.2% V457E (GTG→GAG) ‡ yjeM → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,385,2080TA 5.2% 1198.6 / 35.4 324V457E (GTG→GAG) ‡yjeMputative transporter
Reads supporting (aligned to +/- strand):  ref base T (154/153);  new base A (7/10);  total (161/163)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 349 total aligned reads displayed.

CTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTATCG  >  NC_000913/4385060‑4385356
                                                                                                                                                     |                                                                                                                                                    
cTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTTGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      >  4:179462/1‑149 (MQ=255)
cTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      <  1:119552/149‑1 (MQ=255)
cTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      <  3:1573934/149‑1 (MQ=255)
 tCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                     >  1:1023763/1‑149 (MQ=255)
    cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                                  <  3:704360/149‑1 (MQ=255)
        ggCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                              <  4:1532994/149‑2 (MQ=255)
        ggCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                              <  3:472829/149‑2 (MQ=255)
         gCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaa                                                                                                                                               <  2:847829/147‑1 (MQ=255)
           tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAg                                                                                                                                              >  2:1153288/1‑146 (MQ=255)
           tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaatgat                                                                                                                                           <  2:330106/149‑5 (MQ=255)
                cGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAg                                                                                                                                              <  1:389180/141‑1 (MQ=255)
                 gTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACt                                                                                                                                     >  3:530142/1‑149 (MQ=255)
                  ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                    <  2:1127672/149‑1 (MQ=255)
                     tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                               >  1:1384259/1‑135 (MQ=255)
                     tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                   >  3:1487106/1‑147 (MQ=255)
                     tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                   <  4:1487106/147‑1 (MQ=255)
                     tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaatga                                                                                                                                            >  4:328186/1‑135 (MQ=255)
                      ttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGAc                                                                                                                                <  1:711196/149‑1 (MQ=255)
                         aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCAGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCTGGAGACTGGGACAGc                                                                                                                             >  2:547988/1‑149 (MQ=255)
                          aGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      <  4:555527/123‑1 (MQ=255)
                          aGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      >  3:555527/1‑123 (MQ=255)
                          aGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGca                                                                                                                            <  2:745782/149‑1 (MQ=255)
                             aCGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACAt                                                                                                                         >  1:184844/1‑149 (MQ=255)
                                tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgt                                                                                                                      >  1:931981/1‑149 (MQ=255)
                                tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgt                                                                                                                      >  3:1694562/1‑149 (MQ=255)
                                 cAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                     <  3:356522/149‑1 (MQ=255)
                                   ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTgg                                                                                                                                                    <  2:191546/116‑1 (MQ=255)
                                   ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTgg                                                                                                                                                    <  1:1135765/116‑1 (MQ=255)
                                   ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTgg                                                                                                                                                    >  1:191546/1‑116 (MQ=255)
                                      tctcGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatga                                                                                                                <  4:1017868/149‑1 (MQ=255)
                                         cGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                     <  4:1183074/141‑1 (MQ=255)
                                         cGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTg                                                                                                             <  2:184844/149‑1 (MQ=255)
                                         cGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTg                                                                                                             <  4:1066984/149‑1 (MQ=255)
                                          gACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                               <  4:611100/114‑1 (MQ=255)
                                            cAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTggcg                                                                                                          <  4:684682/149‑1 (MQ=255)
                                                 cGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCt                                                                                                     <  4:1390476/149‑1 (MQ=255)
                                                  gTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGATGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                      <  3:1277668/99‑1 (MQ=255)
                                                   tttGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                     >  2:222906/1‑99 (MQ=255)
                                                      gtgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgaga                                                                                                                                       <  2:1493092/110‑1 (MQ=255)
                                                      gtgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgaga                                                                                                                                       >  1:1493092/1‑110 (MQ=255)
                                                       tgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttc                                                                                               >  3:762040/1‑149 (MQ=255)
                                                        gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttct                                                                                              <  4:1576143/149‑1 (MQ=255)
                                                          ttttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttcttc                                                                                            >  4:1333738/1‑149 (MQ=255)
                                                          ttttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttcttc                                                                                            >  2:669709/1‑149 (MQ=255)
                                                           tttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTg                                                                                                             >  3:1350801/1‑131 (MQ=255)
                                                               aaaaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCt                                                                                       >  4:405544/1‑149 (MQ=255)
                                                               aaaaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCt                                                                                       >  3:1419222/1‑149 (MQ=255)
                                                                 aaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGt                                                                                     <  1:1331099/149‑1 (MQ=255)
                                                                 aaCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGt                                                                                     >  2:25469/1‑149 (MQ=255)
                                                                   cGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTa                                                                                   >  4:1649247/1‑149 (MQ=255)
                                                                       tttGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGa                                                                               <  1:1561644/149‑1 (MQ=255)
                                                                       tttGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGa                                                                               <  3:179462/149‑1 (MQ=255)
                                                                       cactgtcgCAATCATTGCCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCg                                                                                >  1:552202/8‑149 (MQ=255)
                                                                        ttGGCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGAt                                                                              <  1:672439/149‑1 (MQ=255)
                                                                         tGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATg                                                                             <  1:860638/149‑1 (MQ=255)
                                                                            cGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCg                                                                          <  3:916579/149‑1 (MQ=255)
                                                                             ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGa                                                                         <  2:1339434/149‑1 (MQ=255)
                                                                              gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                        >  2:1615083/1‑149 (MQ=255)
                                                                              gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                        >  2:725261/1‑149 (MQ=255)
                                                                                   gATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAcc                                                                   <  2:931981/149‑1 (MQ=255)
                                                                                                 gTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATgg                                                                            <  3:853574/127‑1 (MQ=255)
                                                                                                     tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                     <  3:223630/82‑1 (MQ=255)
                                                                                                     tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCaga                                                  >  2:1450676/1‑149 (MQ=255)
                                                                                                     tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTAAAGCCGGAGACTGGGACAGCACATTGTGGATGATTgg                                                                                                            >  4:1304522/1‑91 (MQ=255)
                                                                                                       gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGa                                                                               >  1:742019/1‑118 (MQ=255)
                                                                                                       gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGt                                                <  3:1847568/149‑1 (MQ=255)
                                                                                                          cTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTg                                               <  3:1326600/147‑1 (MQ=255)
                                                                                                          cTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGc                                             <  3:1075808/149‑1 (MQ=255)
                                                                                                                 tgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                       >  4:1875769/1‑132 (MQ=255)
                                                                                                                  gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGGGGCAAAAAATc                                     >  4:218442/1‑149 (MQ=255)
                                                                                                                     aTTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                        <  2:39603/127‑1 (MQ=255)
                                                                                                                     aTTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGc                                  >  2:861668/1‑149 (MQ=255)
                                                                                                                           gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCaaaaaa                                       <  2:1277372/138‑1 (MQ=255)
                                                                                                                           gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                            >  1:96344/1‑149 (MQ=255)
                                                                                                                             aCGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCg                          >  1:627314/1‑149 (MQ=255)
                                                                                                                               gTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGt                        >  1:1682075/1‑149 (MQ=255)
                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCaga                                                  <  4:1435564/122‑1 (MQ=255)
                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTg                                               >  3:1515057/1‑125 (MQ=255)
                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGTTg                       <  4:1677977/149‑1 (MQ=255)
                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                       >  4:916954/1‑149 (MQ=255)
                                                                                                                                tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                       >  3:59995/1‑149 (MQ=255)
                                                                                                                                 cTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTGTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTgg                      >  1:471726/1‑149 (MQ=255)
                                                                                                                                 cTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATg                                                                             <  3:1671859/94‑1 (MQ=255)
                                                                                                                                  ttCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCgg                         <  1:270275/145‑1 (MQ=255)
                                                                                                                                  ttCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGa                     >  3:1159388/1‑149 (MQ=255)
                                                                                                                                   tCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGaa                    <  2:610185/149‑1 (MQ=255)
                                                                                                                                     aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATa                  >  1:23374/1‑149 (MQ=255)
                                                                                                                                     aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATa                  >  3:816185/1‑149 (MQ=255)
                                                                                                                                      ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATaa                 <  2:627314/149‑1 (MQ=255)
                                                                                                                                       catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAc                                                                    >  4:537672/1‑97 (MQ=255)
                                                                                                                                        atcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAcc               <  4:1773886/149‑1 (MQ=255)
                                                                                                                                         tcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtc                                                                                                  <  3:681485/65‑1 (MQ=255)
                                                                                                                                           atTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa            <  1:1628196/149‑1 (MQ=255)
                                                                                                                                            ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt           >  4:710235/1‑149 (MQ=255)
                                                                                                                                            ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACAAGAACTATTGCAGCAGAGCGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt           >  4:651725/1‑149 (MQ=255)
                                                                                                                                              cAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAAc                                                              <  2:992115/96‑1 (MQ=255)
                                                                                                                                               aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc        >  4:835219/1‑149 (MQ=255)
                                                                                                                                               aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc        >  4:87932/1‑149 (MQ=255)
                                                                                                                                                  cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa     <  4:775575/149‑1 (MQ=255)
                                                                                                                                                  cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCAGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa     >  2:1651126/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTATCG  >  NC_000913/4385060‑4385356

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