Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,529,396 T→C pseudogene (1507/2049 nt) rhsE → RhsE protein in rhs element

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,529,3960TC98.5% 486.9 / ‑3.6 135pseudogene (1507/2049 nt)rhsERhsE protein in rhs element
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (64/69);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCC  >  NC_000913/1529260‑1529530
                                                                                                                                        |                                                                                                                                      
tgaATGGGGTAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGcata                                                                                                                                     >  2:652825/1‑140 (MQ=255)
tgaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGcata                                                                                                                                     >  1:333175/1‑140 (MQ=255)
 gaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTAAAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCatat                                                                                                                                    <  2:63147/140‑1 (MQ=255)
 gaATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCatat                                                                                                                                    >  2:521747/1‑140 (MQ=255)
   aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCa                                                                                                                                  >  1:79144/1‑140 (MQ=255)
   aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCa                                                                                                                                  >  2:1171045/1‑140 (MQ=255)
   aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCa                                                                                                                                  >  1:667769/1‑140 (MQ=255)
   aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCa                                                                                                                                  >  1:1315798/1‑140 (MQ=255)
     ggggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCAcc                                                                                                                                <  1:1129559/140‑1 (MQ=255)
     ggggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCAcc                                                                                                                                >  2:912436/1‑140 (MQ=255)
       ggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCatat                                                                                                                                    >  2:1141354/1‑134 (MQ=255)
       ggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCAcccc                                                                                                                              >  1:1174575/1‑140 (MQ=255)
       ggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGACGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCAcccc                                                                                                                              >  2:1100518/1‑140 (MQ=255)
         cAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCgg                                                                                                                            >  2:1288710/1‑140 (MQ=255)
         cAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCgg                                                                                                                            >  1:262966/1‑140 (MQ=255)
          aaCCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGa                                                                                                                           <  2:7524/140‑1 (MQ=255)
          aaCCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGa                                                                                                                           >  1:987479/1‑140 (MQ=255)
               gCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGa                                                                                                                      <  2:333175/140‑1 (MQ=255)
                cTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGAt                                                                                                                     <  2:534828/140‑1 (MQ=255)
                 ttAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCatat                                                                                                                                    <  2:1348748/124‑1 (MQ=255)
                 ttAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCatat                                                                                                                                    >  1:1348748/1‑124 (MQ=255)
                 ttAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGAtt                                                                                                                    >  1:1066487/1‑140 (MQ=255)
                 ttAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGAtt                                                                                                                    >  1:536734/1‑140 (MQ=255)
                       aggagAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTg                                                                                                              >  1:677594/1‑140 (MQ=255)
                         gagaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTgag                                                                                                            >  1:274928/1‑140 (MQ=255)
                          agaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTgaga                                                                                                           >  2:146850/1‑140 (MQ=255)
                          agaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTgaga                                                                                                           <  1:1219250/140‑1 (MQ=255)
                          agaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTgaga                                                                                                           >  2:660804/1‑140 (MQ=255)
                           gaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGagag                                                                                                          <  2:778450/140‑1 (MQ=255)
                            aaCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGa                                                                                                                           <  1:615705/122‑1 (MQ=255)
                            aaCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGa                                                                                                                           >  2:615705/1‑122 (MQ=255)
                            aaCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAgg                                                                                                         <  1:146850/140‑1 (MQ=255)
                            aaCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAgg                                                                                                         >  1:992330/1‑140 (MQ=255)
                             aCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGt                                                                                                        <  1:660804/140‑1 (MQ=255)
                                cGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                                                                                                     <  1:106452/140‑1 (MQ=255)
                                cGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                                                                                                     <  2:912089/140‑1 (MQ=255)
                                  cATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGtggatg                                                                                                   <  1:652825/140‑1 (MQ=255)
                                  cATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGtggatg                                                                                                   <  2:1362502/140‑1 (MQ=255)
                                      cccTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAat                                                                                               <  1:842449/139‑1 (MQ=255)
                                      aCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGtt                                                                                                               <  1:376570/124‑1 (MQ=255)
                                      aCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGtt                                                                                                               >  2:376570/1‑124 (MQ=255)
                                      aCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGagag                                                                                                          <  2:216553/129‑1 (MQ=255)
                                      aCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGagag                                                                                                          >  1:216553/1‑129 (MQ=255)
                                       ccTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAata                                                                                              <  1:912436/140‑1 (MQ=255)
                                        cTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAatat                                                                                             <  2:1066487/140‑1 (MQ=255)
                                        cTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAatat                                                                                             <  2:954507/140‑1 (MQ=255)
                                        cTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAatat                                                                                             >  1:735191/1‑140 (MQ=255)
                                           ccccAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                                                                                          <  2:1315798/138‑1 (MQ=255)
                                               aGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAg                                                                                      >  2:1240738/1‑140 (MQ=255)
                                                gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGt                                                                                     >  2:397769/1‑140 (MQ=255)
                                                gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGt                                                                                     <  1:299795/140‑1 (MQ=255)
                                                gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGt                                                                                     >  2:79955/1‑140 (MQ=255)
                                                gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGt                                                                                     >  2:639631/1‑140 (MQ=255)
                                                 ccGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTa                                                                                    >  2:594593/1‑140 (MQ=255)
                                                 ccGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTa                                                                                    <  2:274928/140‑1 (MQ=255)
                                                     aCCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCg                                                                                >  1:1096581/1‑140 (MQ=255)
                                                      ccGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGTTTGGGTTGCGAGGTGGTTGGAATTTGTatccgtatccgt                                                                               >  2:538408/1‑140 (MQ=17)
                                                       cGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGtt                                                                              <  2:192888/140‑1 (MQ=255)
                                                          cTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGaa                                                                           <  2:677594/140‑1 (MQ=255)
                                                           tGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTAt                                                                                   <  1:887423/131‑1 (MQ=255)
                                                           tGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAAt                                                                          <  2:1204731/140‑1 (MQ=255)
                                                           tGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAAt                                                                          >  2:493698/1‑140 (MQ=255)
                                                              cAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATccc                                                                       >  2:729293/1‑140 (MQ=255)
                                                               aGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCa                                                                      <  1:1288710/140‑1 (MQ=255)
                                                                gggCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCAt                                                                     <  1:802870/140‑1 (MQ=255)
                                                                  gcagcagcaTGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                                                                   <  1:158395/140‑1 (MQ=255)
                                                                  gcagcagcaTGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                                                                   <  1:1240738/140‑1 (MQ=255)
                                                                  gcagcagcaTGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                                                                   <  1:431613/140‑1 (MQ=255)
                                                                    agcagcaTGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACaa                                                                 >  1:1132992/1‑140 (MQ=255)
                                                                           tgatgaGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                                          <  1:729997/140‑1 (MQ=255)
                                                                           tgatgaGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                                          <  1:521747/140‑1 (MQ=255)
                                                                           tgatgaGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                                          <  2:421933/140‑1 (MQ=255)
                                                                                aggagTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                                          <  1:732035/135‑1 (MQ=255)
                                                                                aggagTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                                          >  2:732035/1‑135 (MQ=255)
                                                                                      cGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATgggg                                               <  1:639631/140‑1 (MQ=255)
                                                                                           cTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                                          <  1:1171045/140‑1 (MQ=255)
                                                                                           cTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                                          <  2:987479/140‑1 (MQ=255)
                                                                                            tGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGAt                                         >  1:310622/1‑140 (MQ=255)
                                                                                            tGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGAt                                         >  1:26036/1‑140 (MQ=255)
                                                                                             gTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATg                                        <  2:1132992/140‑1 (MQ=255)
                                                                                              tACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGc                                       >  1:837099/1‑140 (MQ=255)
                                                                                              tACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGc                                       >  2:206215/1‑140 (MQ=255)
                                                                                                   taACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTg                                  >  2:707353/1‑140 (MQ=255)
                                                                                                     aCCGTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTgag                                >  1:67995/1‑140 (MQ=255)
                                                                                                        gTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGAtt                                   <  1:810784/134‑1 (MQ=255)
                                                                                                        gTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGAtt                                   >  2:810784/1‑134 (MQ=255)
                                                                                                        gTCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAat                             >  2:262982/1‑140 (MQ=255)
                                                                                                         tCACCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAata                            <  2:1096581/140‑1 (MQ=255)
                                                                                                           aCCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGAt                                         <  1:919269/125‑1 (MQ=255)
                                                                                                           aCCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGAt                                         >  2:919269/1‑125 (MQ=255)
                                                                                                           aCCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGa                                     <  2:300491/129‑1 (MQ=255)
                                                                                                           aCCGGCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATg                          >  2:1237050/1‑140 (MQ=255)
                                                                                                              ggCACTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACa                       <  1:594593/140‑1 (MQ=255)
                                                                                                                  cTACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAg                   <  1:896851/140‑1 (MQ=255)
                                                                                                                   tACGATCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAgg                  <  2:1174575/140‑1 (MQ=255)
                                                                                                                       aTCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGa              <  2:330933/140‑1 (MQ=255)
                                                                                                                       aTCCGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGa              >  1:17923/1‑140 (MQ=255)
                                                                                                                          cGTTGCAGGGGCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTa           >  1:306133/1‑140 (MQ=255)
                                                                                                                                 ggggCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATg    >  1:910499/1‑140 (MQ=255)
                                                                                                                                  gggCGGCATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGc   >  2:512295/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAATTTATGCC  >  NC_000913/1529260‑1529530

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