Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,728,215 T→A 7.0% K386I (AAA→ATA)  xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,2150TA 7.0% 530.1 / 24.2 157K386I (AAA→ATA) xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base T (72/74);  new base A (2/9);  total (74/83)
Fisher's exact test for biased strand distribution p-value = 6.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

Only 100 of 164 total aligned reads displayed.

CGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAC  >  NC_000913/3728067‑3728362
                                                                                                                                                    |                                                                                                                                                   
cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                                                                                                                     <  1:135458/149‑1 (MQ=255)
     tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                                                                                                                 <  3:305443/148‑1 (MQ=255)
          gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTCACACTTTGCGGTTTAATACCGCa                                                                                                                                           <  2:90494/149‑1 (MQ=255)
            cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaa                                                                                                                                         <  1:323096/149‑3 (MQ=255)
            cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaa                                                                                                                                         <  1:74162/149‑3 (MQ=255)
            cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaa                                                                                                                                         <  1:189310/149‑3 (MQ=255)
            cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaa                                                                                                                                          <  4:403040/148‑2 (MQ=255)
               ccTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                                                                                                                                       <  3:90582/148‑1 (MQ=255)
                           cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                                                                                                                            <  4:338294/131‑1 (MQ=255)
                            aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                                                                                                          <  4:198810/148‑1 (MQ=255)
                            aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                                                                                                         <  2:87150/149‑1 (MQ=255)
                              tGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                                                                                                        <  4:145751/148‑1 (MQ=255)
                              tGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                                                                                                        <  3:165740/148‑1 (MQ=255)
                               gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGTTCTCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGcc                                                                                                                      >  1:66711/1‑149 (MQ=255)
                               gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                                                                                                                       >  3:352902/1‑148 (MQ=255)
                                   gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                                                                                                                  <  2:266256/149‑1 (MQ=255)
                                    ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                                                                                                                  >  3:121547/1‑148 (MQ=255)
                                     cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                                                                                                         <  2:295311/140‑1 (MQ=255)
                                         caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACTCCAGTACTCACTACGCGCCCCGCCCCCACTCAATGTAACACTTTGCGGTTTAATAGCGCGGGCATGCACGACATCTATGCCATCTGCTAg                                                                                                             >  4:249637/1‑148 (MQ=255)
                                         caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATGCCGCAGGCATGCACGACATCCATGCCATCTGCCAg                                                                                                             >  4:264295/1‑148 (MQ=255)
                                            ccccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCa                                                                                                         >  2:64039/1‑149 (MQ=255)
                                             cccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                                                                                                                     >  1:209857/1‑104 (MQ=255)
                                             cccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                                                                                                                     <  2:209857/104‑1 (MQ=255)
                                                ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa                                                                                                      <  3:181311/148‑1 (MQ=255)
                                                 cccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc                                                                                                     <  3:40868/148‑1 (MQ=255)
                                                   cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa                                                                                                  >  2:105652/1‑149 (MQ=255)
                                                    gTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa                                                                                                  <  3:95357/148‑1 (MQ=255)
                                                             tCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCa                                                                                         >  4:254608/1‑148 (MQ=255)
                                                                aGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc                                                                                           >  4:13129/1‑143 (MQ=255)
                                                                aGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc                                                                                           <  3:13129/143‑1 (MQ=255)
                                                                aGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCa                                                                                      >  4:73443/1‑148 (MQ=255)
                                                                        aCCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgc                                                                              >  3:274077/1‑148 (MQ=255)
                                                                          cGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGcc                                                                           >  1:254625/1‑149 (MQ=255)
                                                                           gCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                                                                                                                                      <  1:272021/89‑1 (MQ=255)
                                                                            cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAGTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa                                                                                                  >  4:402302/1‑124 (MQ=255)
                                                                            cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAg                                                                         >  2:818/1‑149 (MQ=255)
                                                                            cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAg                                                                         >  1:121043/1‑149 (MQ=255)
                                                                             tGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGt                                                                        >  1:100506/1‑149 (MQ=255)
                                                                                tatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTc                                                                      <  3:168769/148‑1 (MQ=255)
                                                                                  tCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAg                                                                         <  4:287960/143‑1 (MQ=255)
                                                                                    cGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                  <  3:767/148‑1 (MQ=255)
                                                                                     gCCAGCATCTGACTCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTgg                                                                 >  3:29056/1‑148 (MQ=255)
                                                                                     gCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTggg                                                                <  2:132936/149‑1 (MQ=255)
                                                                                           aTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATg                                                          >  1:119295/1‑149 (MQ=255)
                                                                                                aCGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtga                                                      <  4:274077/148‑1 (MQ=255)
                                                                                                  gCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgatg                                                    <  3:348732/148‑1 (MQ=255)
                                                                                                   ccAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgatga                                                   >  4:286013/1‑148 (MQ=255)
                                                                                                   ccAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAg                                                  <  2:121043/149‑1 (MQ=255)
                                                                                                    cAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc                                                                                               <  1:114007/103‑1 (MQ=255)
                                                                                                    cAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc                                                                                               >  2:114007/1‑103 (MQ=255)
                                                                                                    cAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAg                                                  <  4:209153/148‑1 (MQ=255)
                                                                                                          tCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAAcc                                           >  2:235233/1‑149 (MQ=255)
                                                                                                          tCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAAcc                                           >  2:236982/1‑149 (MQ=255)
                                                                                                           cACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCa                                                                     <  4:18995/122‑1 (MQ=255)
                                                                                                           cACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCa                                                                     >  3:18995/1‑122 (MQ=255)
                                                                                                           cACTACGCGCCCCGCCCCCAATCAACGTAACACTTGGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTGCCAGCACTGCTCGCGCCAGTTCATTGGGGCCGTGTTGATGAGTCAAAcc                                           >  4:341873/1‑148 (MQ=255)
                                                                                                            aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCa                                          <  4:121547/148‑1 (MQ=255)
                                                                                                            aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGGTCATTGGGGCCATGTTGATGAGGCAAACCaa                                         >  2:195653/1‑149 (MQ=255)
                                                                                                             cTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACTCCTTCCAGCACTGCTAGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaa                                         <  3:168317/148‑1 (MQ=255)
                                                                                                             cTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaaa                                        >  2:62474/1‑149 (MQ=255)
                                                                                                             cTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaa                                         >  3:226449/1‑148 (MQ=255)
                                                                                                              tACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAg                                       >  2:252753/1‑149 (MQ=255)
                                                                                                               aCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAg                                       >  3:347031/1‑148 (MQ=255)
                                                                                                                  cgccccgcccCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGaaa                                    >  4:354703/1‑148 (MQ=255)
                                                                                                                    ccccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAc                                  <  3:375483/148‑1 (MQ=255)
                                                                                                                      ccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAccc                                >  4:169641/1‑148 (MQ=255)
                                                                                                                      ccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAccc                                >  4:133646/1‑148 (MQ=255)
                                                                                                                      ccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAccc                                <  3:233274/148‑1 (MQ=255)
                                                                                                                       cgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtg                                                       <  4:152682/124‑1 (MQ=255)
                                                                                                                        gccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAccccc                              >  4:301643/1‑148 (MQ=255)
                                                                                                                            ccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCt                                                                                 <  3:113091/93‑1 (MQ=255)
                                                                                                                            ccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCt                                                                                 >  4:113091/1‑93 (MQ=255)
                                                                                                                            ccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCg                          >  4:258643/1‑148 (MQ=255)
                                                                                                                            ccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGc                         <  1:121789/149‑1 (MQ=255)
                                                                                                                            ccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGc                         >  2:31348/1‑149 (MQ=255)
                                                                                                                              aaTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGcc                        <  3:302972/148‑1 (MQ=255)
                                                                                                                              aaTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCt                       <  1:64039/149‑1 (MQ=255)
                                                                                                                                tCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgg                     >  1:321242/1‑149 (MQ=255)
                                                                                                                                 cAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgg                     <  3:361286/148‑1 (MQ=255)
                                                                                                                                 cAACGTAACACTTCGCGGTTTAATACCGCAGGCGTGCACGACATCCATGCCATCTGCCAGCGCATAACCCACCCCCTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                    <  2:66711/149‑1 (MQ=255)
                                                                                                                                       aaCACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaa                                         >  3:196114/1‑122 (MQ=255)
                                                                                                                                       aaCACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaa                                         <  4:196114/122‑1 (MQ=255)
                                                                                                                                       aaCACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGattat              >  2:269153/1‑149 (MQ=255)
                                                                                                                                        acacTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGattat              <  3:308025/148‑1 (MQ=255)
                                                                                                                                          agtttGCGGTTTAATACCGCAGGCATGCACCACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtg            <  3:329885/146‑1 (MQ=255)
                                                                                                                                          acTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt           >  1:219801/1‑149 (MQ=255)
                                                                                                                                          acTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt           >  2:277649/1‑149 (MQ=255)
                                                                                                                                          acTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt           >  1:333561/1‑149 (MQ=255)
                                                                                                                                          acTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt           <  1:236982/149‑1 (MQ=255)
                                                                                                                                           cTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCACGCGCCACTTCTTTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtg          >  2:291118/1‑149 (MQ=255)
                                                                                                                                            tttGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgt         >  1:237688/1‑149 (MQ=255)
                                                                                                                                             ttGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgtg        <  1:252753/149‑1 (MQ=255)
                                                                                                                                             ttGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgt         >  3:297675/1‑148 (MQ=255)
                                                                                                                                             ctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                    >  2:17111/2‑137 (MQ=255)
                                                                                                                                             ctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                    <  1:17111/136‑1 (MQ=255)
                                                                                                                                              tGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGatt                <  3:240396/140‑1 (MQ=255)
                                                                                                                                               gCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTgg       <  3:254608/148‑1 (MQ=255)
                                                                                                                                                  gTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                  <  1:211621/86‑1 (MQ=255)
                                                                                                                                                  gTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                  >  2:211621/1‑86 (MQ=255)
                                                                                                                                                  gTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTa   <  2:332036/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAC  >  NC_000913/3728067‑3728362

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: