Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,724 A→C 6.6% T236P (ACC→CCC)  lsrC → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,7240AC 6.6% 296.1 / 10.1 88T236P (ACC→CCC) lsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base C (3/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTC  >  NC_000913/1603576‑1603871
                                                                                                                                                    |                                                                                                                                                   
ttATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGa                                                                                                                                                     <  1:245992/149‑1 (MQ=255)
       aCGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTAc                                                                                                                                              >  3:599186/1‑149 (MQ=255)
           gCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCggg                                                                                                                                          >  1:155600/1‑149 (MQ=255)
            cGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGc                                                                                                                                         <  2:79914/149‑1 (MQ=255)
             gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCt                                                                                                                                        >  2:225888/1‑149 (MQ=255)
             gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                        >  2:419323/1‑149 (MQ=255)
             gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                        >  2:419684/1‑149 (MQ=255)
              aTAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCTTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCACATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTg                                                                                                                                       <  4:654176/149‑1 (MQ=255)
              aTAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTg                                                                                                                                       <  1:586025/149‑1 (MQ=255)
               tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                      >  3:520206/1‑149 (MQ=255)
                aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                     <  4:599186/149‑1 (MQ=255)
                aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                     >  4:246804/1‑149 (MQ=255)
                aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                     <  2:394324/149‑1 (MQ=255)
                  tttACAGGGCGCTCGTCATCTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttg                                                                                                                                   <  3:138657/149‑4 (MQ=255)
                  tttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGgaga                                                                                                                                   >  2:597109/1‑149 (MQ=255)
                    tACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttggg                                                                                                                                 <  1:225888/149‑6 (MQ=255)
                       aGGGCGCTCGTCAACTGGGCGTTCGTATTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                              <  2:194218/149‑1 (MQ=255)
                       aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                              <  2:582856/149‑1 (MQ=255)
                       aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                              >  3:357922/1‑149 (MQ=255)
                       aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                              >  3:273294/1‑149 (MQ=255)
                        gggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                             <  4:554292/149‑1 (MQ=255)
                         ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                  >  3:329876/1‑143 (MQ=255)
                         ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                  <  4:329876/143‑1 (MQ=255)
                          gcgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCa                                                                                                                           >  4:352543/1‑149 (MQ=255)
                           cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCaa                                                                                                                          <  1:419684/149‑1 (MQ=255)
                                    aCTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCAGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCt                                                                                                                 <  3:298523/149‑1 (MQ=255)
                                     cTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTg                                                                                                                >  4:69274/1‑149 (MQ=255)
                                        ggCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGt                                                                                                             <  4:520206/149‑1 (MQ=255)
                                           gTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCt                                                                                                          >  2:148675/1‑149 (MQ=255)
                                           gTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCt                                                                                                          <  4:305071/149‑1 (MQ=255)
                                           gTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCt                                                                                                          <  1:597109/149‑1 (MQ=255)
                                             tCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTgg                                                                                                        <  1:605977/149‑1 (MQ=255)
                                               gTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggc                                                                                                      <  3:10132/149‑1 (MQ=255)
                                                tACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggcg                                                                                                     >  3:658555/1‑149 (MQ=255)
                                                tACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggcg                                                                                                     <  1:206776/149‑1 (MQ=255)
                                                   tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTTCGTGCTGGGCGGCa                                                                                                  >  1:692814/1‑149 (MQ=255)
                                                      aGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCGAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                               <  4:658555/149‑1 (MQ=255)
                                                        ccATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGt                                                                                             <  3:352543/149‑1 (MQ=255)
                                                         cATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTTGTTTTATCCCCAACCAGACCGGCTCCgg                                                                                                                                           >  1:8257/1‑102 (MQ=255)
                                                         cATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCgg                                                                                                                                           <  2:8257/102‑1 (MQ=255)
                                                          aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGttt                                                                                           >  3:82453/1‑149 (MQ=255)
                                                          aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                            >  2:129581/1‑148 (MQ=255)
                                                          aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                            <  1:129581/148‑1 (MQ=255)
                                                              gCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTc                                                                                       <  4:566942/149‑1 (MQ=255)
                                                                aTTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCgg                                                                                     >  1:620199/1‑149 (MQ=255)
                                                                 ttGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                    <  2:511111/149‑1 (MQ=255)
                                                                 ttGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                    <  1:419323/149‑1 (MQ=255)
                                                                      gCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTc                                                                               >  3:589784/1‑149 (MQ=255)
                                                                        aTTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAAATGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCTGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCg                                                                             >  2:185373/1‑149 (MQ=255)
                                                                          tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                           <  2:672452/149‑1 (MQ=255)
                                                                          tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                           <  1:456594/149‑1 (MQ=255)
                                                                          tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                           <  1:179312/149‑1 (MQ=255)
                                                                            tCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGc                                                                         <  2:155600/149‑1 (MQ=255)
                                                                              gTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTg                                                                                          >  4:247238/1‑130 (MQ=255)
                                                                              gTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTg                                                                                          <  3:247238/130‑1 (MQ=255)
                                                                               ttGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTTCTCGGTGGTTCCGGTGCGAt                                                                      >  4:497577/1‑149 (MQ=255)
                                                                                   aCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                               >  1:241309/1‑120 (MQ=255)
                                                                                   aCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                               <  2:241309/120‑1 (MQ=255)
                                                                                    cGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                 <  2:585797/149‑1 (MQ=255)
                                                                                          cATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCgg                                                           >  3:78190/1‑149 (MQ=255)
                                                                                             ggcggcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTAc                                                        >  1:300362/1‑149 (MQ=255)
                                                                                             ggcggcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTAc                                                        >  3:658742/1‑149 (MQ=255)
                                                                                                       gCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATg                                              >  2:591984/1‑149 (MQ=255)
                                                                                                        cGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATgg                                             >  4:642762/1‑149 (MQ=255)
                                                                                                            aaTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTcc                                         >  2:378838/1‑149 (MQ=255)
                                                                                                              ttGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCa                                                                    <  3:497156/120‑1 (MQ=255)
                                                                                                              ttGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCa                                                                    >  4:497156/1‑120 (MQ=255)
                                                                                                              ttGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGGGCGATCCTTGGTGCGGTACTCGGCGCATGGTTCCTg                                       >  4:414317/1‑149 (MQ=255)
                                                                                                                gtgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGAc                                     <  1:148675/149‑1 (MQ=255)
                                                                                                                  gtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGc                                   <  4:73257/149‑1 (MQ=255)
                                                                                                                   tttGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCa                                  >  2:128521/1‑149 (MQ=255)
                                                                                                                     tGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGa                                <  2:4931/149‑1 (MQ=255)
                                                                                                                        ttCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCg                             >  2:274334/1‑149 (MQ=255)
                                                                                                                           gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATa                          <  1:378838/149‑1 (MQ=255)
                                                                                                                                ttGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTa                     >  1:42944/1‑149 (MQ=255)
                                                                                                                                    ttttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTgg                 >  2:322050/1‑149 (MQ=255)
                                                                                                                                     tttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGt                <  3:62326/149‑1 (MQ=255)
                                                                                                                                           cccAACCAGACCGGTACCGGGCTGGAGATGAAAACAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGtt          >  1:237025/1‑149 (MQ=255)
                                                                                                                                            ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  2:669765/1‑149 (MQ=255)
                                                                                                                                            ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  4:486258/1‑149 (MQ=255)
                                                                                                                                            ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  4:106935/1‑149 (MQ=255)
                                                                                                                                            ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  2:386309/1‑149 (MQ=255)
                                                                                                                                            ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:574701/1‑149 (MQ=255)
                                                                                                                                             cAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgc        <  3:246804/149‑1 (MQ=255)
                                                                                                                                              ctccAGCCCGGTATGGGTATGGGGATGAAAGCAATTGCAGCCTGCTTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGAGGCAGATCGATAGCGTACTGGTGCTGTTgcg       >  4:82579/3‑149 (MQ=255)
                                                                                                                                                ccAGCCCGGTAGCGGTCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCa     >  3:543888/1‑149 (MQ=255)
                                                                                                                                                ccAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCa     >  3:148423/1‑149 (MQ=255)
                                                                                                                                                   gACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTc  <  3:69274/149‑1 (MQ=255)
                                                                                                                                                   gACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTc  <  1:86871/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTC  >  NC_000913/1603576‑1603871

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