Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,736 T→A 10.4% K59* (AAA→TAA) ‡ PP_4851 ← phosphate starvation‑inducible protein PsiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,518,7360TA10.4% 183.0 / 9.5 67K59* (AAA→TAA) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base A (3/4);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCA  >  NC_002947/5518590‑5518876
                                                                                                                                                    |                                                                                                                                            
tCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTACATATTTTCTAGCTGCGATTGTGTCACATCCTATTAAAGGCAGGTCATAATCAACGCTTTGc                                                                                                                                            >  8:349096/1‑149 (MQ=255)
  aTTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCCTTGGTGGTGGCGTCAGCATTGCAGGTCTACATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGt                                                                                                                                          <  2:224757/149‑1 (MQ=255)
   tttCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGt                                                                                                                                          <  1:324949/148‑1 (MQ=255)
     tcttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGGTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCa                                                                                                                                       >  4:2540/1‑149 (MQ=255)
       ttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATc                                                                                                                                     >  3:18744/1‑149 (MQ=255)
        attctaGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCg                                                                                                                                    <  7:166798/144‑1 (MQ=255)
         cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTTTGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATAAACGCAAAGCGTTCAtga                                                                                                                                    <  1:126965/148‑3 (MQ=255)
         cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAtga                                                                                                                                    <  5:357076/148‑3 (MQ=255)
              ggCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTg                                                                                                                              <  4:18744/149‑1 (MQ=255)
              ggCAATAGCTAATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGGATTGCAGGTCTTCATTTTTTATTGCTGGGTTTGGGGCGCATCGTTTTTAAGGCAGGTGCTCATGAACGCTTTGCGTtcgtcg                                                                                                                                    >  8:166798/1‑143 (MQ=255)
                cAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTGTTCATTTTTTCTTGCGGGGTTTGGGGGCCATCCTTTTGGAGGCAGGTGCTCATGAAAGCTTTGCGTTCGTCGCCCTTgcg                                                                                                                            >  7:363534/1‑147 (MQ=255)
                cAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAg                                                                                                                            <  4:120465/149‑1 (MQ=255)
                     gCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGAGTTCATCGCCCTTGAGGGCCt                                                                                                                       >  3:252333/1‑149 (MQ=255)
                      cTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCtt                                                                                                                      >  3:77005/1‑149 (MQ=255)
                      cTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCtt                                                                                                                      <  4:77005/149‑1 (MQ=255)
                       tCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGGGTTCATCGCACTTGAGGGGCtt                                                                                                                      >  1:68482/1‑149 (MQ=255)
                           gAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCa                                                                                                                                       >  5:93422/1‑127 (MQ=255)
                           gAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCa                                                                                                                                       <  6:93422/127‑1 (MQ=255)
                                       cTCATCCCCCTTCAGGGCTTTGGTGGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAtga                                                                                                                                    >  1:196655/1‑117 (MQ=255)
                                       cTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAtga                                                                                                                                    <  2:196655/118‑3 (MQ=255)
                                         cATCCCCCTTCCGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTGTGCTTGCTGTGTTTGCGGCACATCCTTATTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTg                                                                                                   <  6:75545/149‑1 (MQ=255)
                                          aTCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTg                                                                                                   >  5:129538/1‑148 (MQ=255)
                                             ccccTTCAGGGCTTTGGTGGTTGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAg                                                                                                <  1:372295/148‑1 (MQ=255)
                                             ccccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGTGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAg                                                                                                <  4:369332/148‑1 (MQ=255)
                                              cccTTCAGGGAGTTGGTGGTGGCGGCAGCATTGCAGTTCTTCATTTTTTCTTGCTGGGGTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGt                                                                                              <  1:243685/149‑1 (MQ=255)
                                                   cAGGGCTTTGTTGGTGGCCGCAGCATGGCAGCTCTTCATTTTTTCTTGCTCCGTTTGGGGCACATCCTTTTTCAGCCTGGTTCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCa                                                                                         <  2:68482/149‑1 (MQ=255)
                                                                tggttccGTCAGCATTGCAGGTCTTCATTTTTTCTTTCAGGGTTTGGGTCACCTACTTTATCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCctg                                                                               <  4:252333/142‑1 (MQ=255)
                                                                             aTTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTgg                                                                 <  7:135630/149‑1 (MQ=255)
                                                                                      cTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAc                                                        >  8:301063/1‑149 (MQ=255)
                                                                                       ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAcc                                                       >  2:46944/1‑149 (MQ=255)
                                                                                          aTTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGc                                                                                                  >  6:310112/1‑103 (MQ=255)
                                                                                           ttttttCTTGCGTGGTTTGGGTCACTTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGGGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                   <  7:349096/149‑1 (MQ=255)
                                                                                                 cTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACt                                              <  3:2540/148‑1 (MQ=255)
                                                                                                  ttGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCATTTCCTGCTGCGCGGTATTGTCAGCGAAACCTTTGGCACTg                                             >  3:202698/1‑148 (MQ=255)
                                                                                                   ggcttggcttgTGGCCCATCCTTTTTCATGCAGGAGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaa                                           <  8:363534/148‑1 (MQ=255)
                                                                                                              gggCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgc                                  >  2:194683/1‑147 (MQ=255)
                                                                                                                 cacaTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGcc                                                                                                                        <  6:188281/58‑1 (MQ=255)
                                                                                                                 cacaTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGcc                                                                                                                        >  5:188281/1‑58 (MQ=255)
                                                                                                                  acaTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgc                            >  7:86748/1‑149 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGcc                                                                                                                        >  1:149377/1‑55 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGcc                                                                                                                        <  2:149377/55‑1 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgcca                          <  6:129538/149‑1 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgc                            >  1:83637/1‑147 (MQ=255)
                                                                                                                     tCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                         >  2:163171/1‑149 (MQ=255)
                                                                                                                      ccTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                         >  7:73860/1‑148 (MQ=255)
                                                                                                                       ctatttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTAATCTTTTCCTGCTCCGCGGTATTGGCACCGAAAACATGGGCAAGTAAAAACACCGCCACCGCCAGt                        >  2:377129/4‑148 (MQ=255)
                                                                                                                           ttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGAAGCGAAACCTTGGGCACTGAACAACACCGCCCGCGCCAGTAACgg                   >  3:212251/1‑149 (MQ=255)
                                                                                                                             caggcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGa                 >  4:144577/1‑149 (MQ=255)
                                                                                                                               ggcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGAtt               <  2:83637/149‑1 (MQ=255)
                                                                                                                                gcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACggg                  >  8:60712/1‑145 (MQ=255)
                                                                                                                                gcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACggg                  <  7:60712/145‑1 (MQ=255)
                                                                                                                                  aggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgc            <  8:86748/149‑1 (MQ=255)
                                                                                                                                   ggggCTCATGAACTCCTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgca           <  7:334595/146‑1 (MQ=255)
                                                                                                                                   ggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgca           <  3:144577/149‑1 (MQ=255)
                                                                                                                                      gCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAt                                                                    >  5:268180/1‑89 (MQ=255)
                                                                                                                                      gCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAt                                                                    <  6:268180/89‑1 (MQ=255)
                                                                                                                                        tCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACt      >  3:83964/1‑149 (MQ=255)
                                                                                                                                         gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTTGCAGCGACACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt     >  1:290838/2‑149 (MQ=255)
                                                                                                                                         gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGa                 >  5:180095/2‑137 (MQ=255)
                                                                                                                                         gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGa                 <  6:180095/136‑1 (MQ=255)
                                                                                                                                         cATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt     >  5:79477/1‑149 (MQ=255)
                                                                                                                                          aTGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt    >  4:46042/1‑149 (MQ=255)
                                                                                                                                          aTGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt    >  8:197158/1‑149 (MQ=255)
                                                                                                                                          aTGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt    >  6:98726/1‑149 (MQ=255)
                                                                                                                                            gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACt      <  5:163893/145‑1 (MQ=255)
                                                                                                                                            gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACt      >  6:163893/1‑145 (MQ=255)
                                                                                                                                             aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCa  <  7:301063/148‑1 (MQ=255)
                                                                                                                                                  tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagc                              <  6:92391/115‑1 (MQ=255)
                                                                                                                                                  tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagc                              >  5:92391/1‑115 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCA  >  NC_002947/5518590‑5518876

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