Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,874,676 T→C 5.4% E178G (GAA→GGA)  rsuA ← pseudouridylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,874,6760TC 5.4% 220.5 / 3.6 75E178G (GAA→GGA) rsuApseudouridylate synthase
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base C (2/2);  total (34/41)
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.86e-01

ATCATTTTCAGTTAATAAACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGATTTGCTGAAATAACTTCATA  >  CP000730/1874537‑1874815
                                                                                                                                           |                                                                                                                                           
aTCATTTTCAGTTAATAAACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTc                                                                                                                                            <  2:699896/141‑1 (MQ=255)
 tCATTTTCAGTTAATAAACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCa                                                                                                                                           >  2:7291/1‑141 (MQ=255)
  cATTTTCAGTTAATAAACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCAt                                                                                                                                          <  2:735876/141‑1 (MQ=255)
            taataaACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTAc                                                                                                                                <  1:692447/141‑1 (MQ=255)
            taataaACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTAc                                                                                                                                <  1:433775/141‑1 (MQ=255)
             aataaACGATATTCTCCTGAATCTAAATTGCTGTCAAATTCACAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTACTTCATAAATCGTTCCg                                                                                                                               >  2:133090/1‑141 (MQ=255)
              ataaACCATATTCTCCTCAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACgt                                                                                                                              <  2:160988/141‑1 (MQ=255)
               taaACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACgtg                                                                                                                             <  1:481998/141‑1 (MQ=255)
               taaACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACgtg                                                                                                                             <  2:494001/141‑1 (MQ=255)
                aaaCGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACgtgt                                                                                                                            <  1:585522/141‑1 (MQ=255)
                aaaCGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACgtgt                                                                                                                            >  1:27204/1‑141 (MQ=255)
                 aaCGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGtgtg                                                                                                                           >  1:523260/1‑141 (MQ=255)
                    gATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACg                                                                                                                        <  1:98543/141‑1 (MQ=255)
                       atTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGttt                                                                                                                     <  2:600252/141‑1 (MQ=255)
                       atTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTGCATAACTTGTTACGTGTGACGttt                                                                                                                     >  2:352005/1‑141 (MQ=255)
                       atTCTCCTGAATCTAAATTGCTGTCCAATTCGAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGttt                                                                                                                     <  1:443330/141‑1 (MQ=255)
                       atTCCCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGttt                                                                                                                     <  1:472664/141‑1 (MQ=255)
                         tctcCTGTATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATGTTGACTGCTAGGAAACATGCGCTTATCTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGa                                                                                                                   <  2:46423/141‑1 (MQ=255)
                          ctcCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGAt                                                                                                                  >  2:583795/1‑141 (MQ=255)
                            ccTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTTTTACGTGTGACGTTTGATTa                                                                                                                >  1:727291/1‑141 (MQ=255)
                            ccTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCCTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTa                                                                                                                >  1:513040/1‑141 (MQ=255)
                                aaTCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATCATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCaa                                                                                                            <  1:459959/141‑1 (MQ=255)
                                  tCTATATTGCTTTCCAATTCTACATCAGCATTTTTTATACGTCTTCAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAata                                                                                                          <  1:168001/141‑1 (MQ=255)
                                  tCTAAATTGCTGTCGAATTCTAAATCAGCAATTTTGATACATCTTAAATGTTATACCCCATTTTGAATGCTATGATACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACTTGTGACGTTTGATTATCAata                                                                                                          <  2:219630/141‑1 (MQ=255)
                                  acTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAata                                                                                                          <  2:419582/140‑1 (MQ=255)
                                      aagtGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTAAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCATTCATAAACTATTACGTGTGACCTTTGATTATCacattaca                                                                                                      >  2:338352/4‑134 (MQ=255)
                                      aaTTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATTTAAg                                                                                                      <  1:550307/141‑1 (MQ=255)
                                      aaTTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAg                                                                                                      >  2:13424/1‑141 (MQ=255)
                                      aaTTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACATCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATCGTTACGTGTGACGTTTGATTCTCAATATAAg                                                                                                      >  2:82358/1‑141 (MQ=255)
                                           cTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTCTAACTTGGTGTTAGCTTCCTTCATATCTTGTTGCTTGTGACGTTTGATTATCAATATAAGTTAat                                                                                                 <  2:170746/141‑1 (MQ=255)
                                           cTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTCTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAat                                                                                                 <  2:734103/141‑1 (MQ=255)
                                                 aactcaaaaTCAGCAATTTTCATACTTCTTAAATGTAATACCTCACTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCa                                                                                           <  2:108622/135‑1 (MQ=255)
                                                 aaTTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATATCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCa                                                                                           <  2:412249/141‑1 (MQ=255)
                                                 aaTTCCAAATCAGCAATTTTGATACGTCTTAAACGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATCATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCa                                                                                           <  2:667120/141‑1 (MQ=255)
                                                     cTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTCCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCt                                                                                       >  2:107078/1‑141 (MQ=255)
                                                     cTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTCCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCt                                                                                       >  2:664256/1‑141 (MQ=255)
                                                          tCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAAcc                                                                                  >  1:362615/1‑141 (MQ=255)
                                                            aGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCtt                                                                                <  1:404415/141‑1 (MQ=255)
                                                              taaTTTTGATACGTCTCAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGc                                                                              <  2:703856/140‑1 (MQ=255)
                                                              aaaTTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGc                                                                              <  2:52366/140‑1 (MQ=255)
                                                                  tttGATACGTCTTAAATGTAATAACTCGTTTTGAATGCTATGAAACATTCGTTTAACTCGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTAGCAATATACGTTAATATTGCAGGCTTAACCTTGCCAtc                                                                          <  2:603910/141‑1 (MQ=255)
                                                                  tttGATACATCTTAAATGTATAACCGCTGGGTGAATGCTATGATACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGATTAACCTTGCCAtc                                                                          <  2:603985/141‑1 (MQ=255)
                                                                        aCGTTTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCAtcagtcag                                                                    <  2:63793/141‑1 (MQ=255)
                                                                            ctcaaaTGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGtt                                                                <  2:94766/138‑1 (MQ=255)
                                                                                  tGTAATACCGTATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACAccc                                                          <  1:10452/141‑1 (MQ=255)
                                                                                    tAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCctc                                                        >  2:405614/1‑141 (MQ=255)
                                                                                    tAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCctc                                                        >  2:250065/1‑141 (MQ=255)
                                                                                    tAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCACTGTTACACCctc                                                        >  2:263386/1‑141 (MQ=255)
                                                                                      aTACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCtt                                                      >  2:132636/1‑141 (MQ=255)
                                                                                       tACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTCCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCttt                                                     <  1:664256/141‑1 (MQ=255)
                                                                                       tACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTCCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCttt                                                     <  1:107078/141‑1 (MQ=255)
                                                                                        aCCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTa                                                    >  2:704034/1‑141 (MQ=255)
                                                                                         ccTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTCCACCCTCTTTaa                                                   >  2:245336/1‑141 (MQ=255)
                                                                                         ccTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTaa                                                   <  1:7291/141‑1 (MQ=255)
                                                                                         ccTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTaa                                                   >  2:767660/1‑141 (MQ=255)
                                                                                         ccTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTaa                                                   >  2:207562/1‑141 (MQ=255)
                                                                                            cATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTaaaaa                                                >  2:682410/1‑140 (MQ=255)
                                                                                            cATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAAGCTTAACCTTGCCCTCAGTCAGTGTTACACCCTCTTTAAAAg                                                >  1:238348/1‑141 (MQ=255)
                                                                                                 tGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGa                                           <  2:236969/141‑1 (MQ=255)
                                                                                                   aaTGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAATTTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAAt                                         >  1:427953/1‑141 (MQ=255)
                                                                                                   aaTGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAAt                                         >  1:499255/1‑141 (MQ=255)
                                                                                                    aTGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATg                                        >  2:546552/1‑141 (MQ=255)
                                                                                                    aTGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATg                                        >  2:731140/1‑141 (MQ=255)
                                                                                                       cTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcg                                     >  2:197646/1‑141 (MQ=255)
                                                                                                       cTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcg                                     <  1:690927/141‑1 (MQ=255)
                                                                                                       cTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTCGCGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcg                                     <  1:424297/141‑1 (MQ=255)
                                                                                                        tATGAAACGTTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCATTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcgt                                    <  2:513040/141‑1 (MQ=255)
                                                                                                        tATGAAACATTCGTTTAACTTGATTATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcgt                                    <  2:727291/141‑1 (MQ=255)
                                                                                                         aTGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATgtcgtc                                   <  2:34273/141‑1 (MQ=255)
                                                                                                                   cGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATa                         >  1:189295/1‑141 (MQ=255)
                                                                                                                      ttAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGa                      >  2:98317/1‑141 (MQ=255)
                                                                                                                      ttAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGa                      >  2:776329/1‑141 (MQ=255)
                                                                                                                      ttAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGa                      >  1:416286/1‑141 (MQ=255)
                                                                                                                          cTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGATTTg                  <  2:603395/141‑1 (MQ=255)
                                                                                                                           ttGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGACATCAGTCAGTGTTACACCCTCTTTAAAAGCTGGAATGTCGTCTTCAGTGATAGGATTTGc                 >  1:65206/1‑141 (MQ=255)
                                                                                                                                          ttCTTAAATTGTGACGTGTGACGTTTGATTATTAATATAAGTTAATATTGCAGGCCTAACCTTGCCATCAGTGAGTGTTACACCCTCTTTAATAGCTTGAATGTCGTCTTCAGTGATAGGATTTGCTGATATAACTTCatt  <  2:432421/141‑2 (MQ=255)
                                                                                                                                           |                                                                                                                                           
ATCATTTTCAGTTAATAAACGATATTCTCCTGAATCTAAATTGCTGTCCAATTCTAAATCAGCAATTTTGATACGTCTTAAATGTAATACCTCATTTTGAATGCTATGAAACATTCGTTTAACTTGATGATATTTTCCTTCATAAATTGTTACGTGTGACGTTTGATTATCAATATAAGTTAATATTGCAGGCTTAACCTTGCCATCAGTCAGTGTTACACCCTCTTTAAAAGCTTGAATGTCGTCTTCAGTGATAGGATTTGCTGAAATAACTTCATA  >  CP000730/1874537‑1874815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: