Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,091 T→A 9.4% intergenic (+21/+28) pepT → / ← roxA peptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0910TA 9.4% 263.4 / 16.2 86intergenic (+21/+28)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base T (44/33);  new base A (2/6);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGA  >  NC_000913/1186943‑1187234
                                                                                                                                                    |                                                                                                                                               
cGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                                                                 <  4:467478/149‑1 (MQ=255)
  tGCCCGAACCTGTTCACTGGCGGTTACAACTACCATGGTAAGCATGAGTTGGTGACTCTGGCAGGTATGGAAAAAGCGGTGCAGGTGATCCTCCGTATTGCCGACTTAACGGCGCAACGGAAGTAAGCGAAACGGGATGCGGCATGTGa                                                                                                                                               <  3:741956/149‑1 (MQ=255)
     ccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                            >  2:431189/1‑149 (MQ=255)
     ccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                            >  4:116922/1‑149 (MQ=255)
     ccGAACCTGTTCACTGGCGCTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                            <  2:466593/149‑1 (MQ=255)
          ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                       <  1:179526/149‑1 (MQ=255)
          ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                       <  3:665536/149‑1 (MQ=255)
          ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCAt                                                                                                                                       <  4:52442/149‑1 (MQ=255)
           cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                                                                      >  4:626845/1‑149 (MQ=255)
           cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATc                                                                                                                                      >  1:448918/1‑149 (MQ=255)
           cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                                                                      <  2:442117/149‑1 (MQ=255)
            tGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                                                                     <  2:192155/149‑1 (MQ=255)
                  cTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTTTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCTTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGGGAAAAGGGATGCGGCATCTGATGCCGCATCCGGCTTa                                                                                                                               >  2:766097/1‑149 (MQ=255)
                   tGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaa                                                                                                                              >  3:288417/1‑149 (MQ=255)
                      cGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATc                                                                                                                           >  3:417910/1‑149 (MQ=255)
                         ttACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                        <  2:542820/133‑1 (MQ=255)
                         ttACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCa                                                                                                                                        >  1:542820/1‑133 (MQ=255)
                             aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACtt                                                                                                                    >  1:526572/1‑149 (MQ=255)
                             aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACtt                                                                                                                    >  3:521221/1‑149 (MQ=255)
                             aaCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACtt                                                                                                                    >  3:751265/1‑149 (MQ=255)
                               cTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAc                                                                                                                  <  3:116922/149‑1 (MQ=255)
                                tATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                                                                 >  2:892687/1‑149 (MQ=255)
                                  tCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCt                                                                                                               >  4:304634/1‑149 (MQ=255)
                                           gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAc                                                                                                      <  1:892687/149‑1 (MQ=255)
                                             aTGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCa                                                                                                    >  1:197634/1‑149 (MQ=255)
                                               gAGTTTGTGACTCTGGAAGGTATGAAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAt                                                                                                  >  4:81103/1‑149 (MQ=255)
                                                  tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                                                               >  1:24088/1‑149 (MQ=255)
                                                      tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                                                           >  1:349088/1‑101 (MQ=255)
                                                      tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                                                           <  2:349088/101‑1 (MQ=255)
                                                         ctctGGAAGGTATAGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTAtt                                                                                        <  4:751265/149‑1 (MQ=255)
                                                          tctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                             <  2:762687/127‑1 (MQ=255)
                                                          tctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                             >  1:762687/1‑127 (MQ=255)
                                                            tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                                                                     <  4:39355/149‑1 (MQ=255)
                                                            tGGAAGGTATGAAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAc                                                                                     <  2:492521/149‑1 (MQ=255)
                                                              gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAc                                                                                                      >  4:109205/1‑130 (MQ=255)
                                                              gAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAc                                                                                                      <  3:109205/130‑1 (MQ=255)
                                                                 ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAAACCCGCTATTGACCAgcg                                                                                >  1:530682/1‑149 (MQ=255)
                                                                     tGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGcgccgc                                                                            >  4:196039/1‑149 (MQ=255)
                                                                      ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCcgcg                                                                           >  4:749877/1‑149 (MQ=255)
                                                                        aaaaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAg                                                                         <  1:431189/149‑1 (MQ=255)
                                                                            aGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATc                                                                     <  3:691129/149‑1 (MQ=255)
                                                                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGa                                                                 >  4:180726/1‑149 (MQ=255)
                                                                                 tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                <  4:63123/149‑1 (MQ=255)
                                                                                  gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                             <  4:763157/103‑1 (MQ=255)
                                                                                  gCAGGTGATCGGCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                             >  3:763157/1‑103 (MQ=255)
                                                                                   cAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTCCCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                              <  2:24088/149‑1 (MQ=255)
                                                                                   cAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                              <  4:253241/149‑1 (MQ=255)
                                                                                   cAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                              >  3:109836/1‑149 (MQ=255)
                                                                                     ggTGATCGTCCGTATTGCAGAGTTAACGGAGCAACGGAAGTAAGCGAAAAGGGATGCAGCATGTGATGCCGCATCCGACTTAATTCCAAGCTAACCCCTCGAAGAACCAATTCACGCTATTGACCAGCGCCGCGAGAATCGCGTGCAAt                                                            <  2:66673/149‑1 (MQ=255)
                                                                                      gTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt                                                            >  4:648216/1‑148 (MQ=255)
                                                                                      gTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt                                                            <  3:648216/148‑1 (MQ=255)
                                                                                      gTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg                                                           >  3:60284/1‑149 (MQ=255)
                                                                                       tGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGa                                                          <  2:448918/149‑1 (MQ=255)
                                                                                        gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc                                                         <  1:11784/149‑1 (MQ=255)
                                                                                         aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg                                                        >  1:777119/1‑149 (MQ=255)
                                                                                           cGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGa                                                      <  1:395154/149‑1 (MQ=255)
                                                                                            gTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCaa                                                                                                   <  1:320827/103‑1 (MQ=255)
                                                                                            gTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCaa                                                                                                   >  2:320827/1‑103 (MQ=255)
                                                                                             tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc                                                    >  4:569265/1‑149 (MQ=255)
                                                                                             tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc                                                    >  2:818881/1‑149 (MQ=255)
                                                                                                 tATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc                                                >  1:152526/1‑149 (MQ=255)
                                                                                                  aTTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCACTACCCGCTATTCACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCa                                               <  4:469695/149‑1 (MQ=255)
                                                                                                      ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcgc                                           >  1:279110/1‑149 (MQ=255)
                                                                                                        gAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCAt                                         >  4:483748/1‑149 (MQ=255)
                                                                                                             aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCa                                    >  3:660011/1‑149 (MQ=255)
                                                                                                              aCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaa                                   >  3:561652/1‑149 (MQ=255)
                                                                                                               cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaaa                                  <  4:106174/149‑1 (MQ=255)
                                                                                                                ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc                                                         <  1:553627/125‑1 (MQ=255)
                                                                                                                ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc                                                         >  2:553627/1‑125 (MQ=255)
                                                                                                                   gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtc                              <  4:191018/149‑1 (MQ=255)
                                                                                                                   gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtc                              >  1:878772/1‑149 (MQ=255)
                                                                                                                   gcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtc                              <  3:304634/149‑1 (MQ=255)
                                                                                                                      aCGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTcc                           >  1:861491/1‑149 (MQ=255)
                                                                                                                       cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCg                          >  4:773402/1‑149 (MQ=255)
                                                                                                                       cGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCg                          >  3:152629/1‑149 (MQ=255)
                                                                                                                         gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCa                        <  1:612411/149‑1 (MQ=255)
                                                                                                                         gAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCa                        <  3:626845/149‑1 (MQ=255)
                                                                                                                             tAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCa                    >  1:56139/1‑149 (MQ=255)
                                                                                                                              aaGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAg                   >  2:298511/1‑149 (MQ=255)
                                                                                                                                   aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCaa              <  2:279110/149‑1 (MQ=255)
                                                                                                                                      aGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt           <  3:389770/149‑1 (MQ=255)
                                                                                                                                         gATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGc        >  1:499036/1‑149 (MQ=255)
                                                                                                                                          aTGCGGCATGTGATGCCGCATCCTGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCCCCCAGCGCATCGCCAACATTATCCGCAGTCAGCGCAATGTTGCt       <  2:884519/149‑1 (MQ=255)
                                                                                                                                          aTGCGGCATCCCATGCCGCATGCGACTTAAATCCAAACTTACTCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt       >  2:743980/1‑149 (MQ=255)
                                                                                                                                              ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCg   >  3:359364/1‑149 (MQ=255)
                                                                                                                                               gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGa  <  3:566656/149‑1 (MQ=255)
                                                                                                                                               gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGa  >  4:781758/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGA  >  NC_000913/1186943‑1187234

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