Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,578,008 G→C 5.6% R80R (CGC→CGG) ‡ ttgC ← efflux pump outer membrane protein TtgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,578,0080GC 5.6% 282.5 / 4.2 87R80R (CGC→CGG) ‡ttgCefflux pump outer membrane protein TtgC
Reads supporting (aligned to +/- strand):  ref base G (36/46);  new base C (3/2);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TGGTGATGCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCGGCCACTGCGCAG  >  NC_002947/1577863‑1578153
                                                                                                                                                 |                                                                                                                                                 
tggTGATGCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGAcgcg                                                                                                                                                 >  4:50004/1‑148 (MQ=255)
 ggTGATGCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGAcgcg                                                                                                                                                 >  3:178530/1‑147 (MQ=255)
 ggTGATGCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGAcgcg                                                                                                                                                 >  1:173282/1‑147 (MQ=255)
       gCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTcgc                                                                                                                                         >  8:119345/1‑149 (MQ=255)
        ccAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTTGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACCTGCGCGTcgcg                                                                                                                                        <  6:245945/149‑1 (MQ=255)
        ccAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTcgcg                                                                                                                                        >  1:101978/1‑149 (MQ=255)
          aGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGt                                                                                                                                      <  2:288581/149‑1 (MQ=255)
             ttCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTAt                                                                                                                                   >  8:248738/1‑149 (MQ=255)
             ttCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACACCGCGACCTGCGCGTCGCggcgt                                                                                                                                    <  8:79902/148‑4 (MQ=255)
              tCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCgg                                                                                                                                       <  2:81733/144‑1 (MQ=255)
              tCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCgg                                                                                                                                       >  1:81733/1‑144 (MQ=255)
                 cccGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGAc                                                                                                                               <  8:97114/149‑1 (MQ=255)
                 cccGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGAc                                                                                                                               <  3:50004/149‑1 (MQ=255)
                        gCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCACGTCGCGGTTATTGACCAGCGAg                                                                                                                        >  7:272610/1‑149 (MQ=255)
                          gACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGt                                                                                                                      <  8:214947/149‑1 (MQ=255)
                           aCATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTcgcg                                                                                                                                        <  7:9002/130‑1 (MQ=255)
                           aCATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTcgcg                                                                                                                                        >  8:9002/1‑130 (MQ=255)
                           aCATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTc                                                                                                                     >  5:227149/1‑149 (MQ=255)
                             aTGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTg                                                                                                                   <  4:123888/149‑1 (MQ=255)
                              tGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCATGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  <  7:119345/149‑1 (MQ=255)
                              tGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  <  5:169110/149‑1 (MQ=255)
                              tGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  <  6:227149/149‑1 (MQ=255)
                              tGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTg                                                                                                                   >  8:301972/1‑148 (MQ=255)
                              tGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGGCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  <  8:203850/149‑1 (MQ=255)
                                   gCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCa                                                                                                              <  2:101978/148‑1 (MQ=255)
                                   gCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAg                                                                                                             >  2:169652/1‑149 (MQ=255)
                                    ccGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGc                                                                                                            <  8:200952/149‑1 (MQ=255)
                                     cGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAg                                                                                                                        >  6:105874/1‑136 (MQ=255)
                                     cGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAg                                                                                                                        <  5:105874/136‑1 (MQ=255)
                                      gggACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTg                                                                                                          >  6:178385/1‑149 (MQ=255)
                                       ggACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGt                                                                                                         >  7:91930/1‑149 (MQ=255)
                                       ggACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGt                                                                                                         <  2:171696/149‑1 (MQ=255)
                                         acacGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTg                                                                                                       <  7:230548/149‑1 (MQ=255)
                                         acacGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTg                                                                                                       <  5:178385/149‑1 (MQ=255)
                                         acacGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCCGGAACAGGTCGGCGCGCTGAATGAGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGATATTGACCAGCGAGGTCTGGATAAGCTGTTg                                                                                                       >  6:158118/1‑149 (MQ=255)
                                          cacGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGAcc                                                                                                                              >  6:141415/1‑125 (MQ=255)
                                          cacGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGAcc                                                                                                                              <  5:141415/125‑1 (MQ=255)
                                           acGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCa                                                                                                     <  7:118285/149‑1 (MQ=255)
                                                 gcggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  >  2:131463/1‑130 (MQ=255)
                                                 gcggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTgg                                                                                                                  <  1:131463/130‑1 (MQ=255)
                                                 gcggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAg                                                                                                             >  6:15156/1‑135 (MQ=255)
                                                 gcggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAg                                                                                                             <  5:15156/135‑1 (MQ=255)
                                                  cggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCgg                                                                                              <  6:286579/149‑1 (MQ=255)
                                                   ggcTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGAGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCggg                                                                                             >  4:100502/1‑149 (MQ=255)
                                                    gctgccgctgccGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGt                                                                                            <  3:220131/149‑1 (MQ=255)
                                                          ggtgccgGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTTTTGCGCACGGTAGGCGTCGAGGTTCTACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGGTGTTGCAGCGCCGGGTCGTGga                                                                                      <  5:158118/147‑1 (MQ=255)
                                                              ccgGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGTATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGgaa                                                                                     <  2:259871/146‑1 (MQ=255)
                                                              ccgGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGTATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGgaa                                                                                     >  1:259871/1‑146 (MQ=255)
                                                              ccgGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGgaagaa                                                                                  <  1:185238/149‑1 (MQ=255)
                                                                ggTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACt                                                                                <  1:11055/149‑1 (MQ=255)
                                                                   ggCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGgc                                                                             <  4:178530/149‑1 (MQ=255)
                                                                    gCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGc                                                                                                      >  5:57193/1‑123 (MQ=255)
                                                                    gCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGc                                                                                                      <  6:57193/123‑1 (MQ=255)
                                                                    gCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGgcg                                                                            <  2:94433/149‑1 (MQ=255)
                                                                      cGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAg                                                                                                    <  5:20430/123‑1 (MQ=255)
                                                                      cGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAg                                                                                                    >  6:20430/1‑123 (MQ=255)
                                                                          aCCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGcc                                                                      <  4:244120/149‑1 (MQ=255)
                                                                            cgccgGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCt                                                                    >  6:101030/1‑149 (MQ=255)
                                                                             gccgGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTg                                                                   >  3:264768/1‑149 (MQ=255)
                                                                              ccgGGAACACGTCGGCGCGCTGAATGCGGTATTGCGCACGGGAGGCGTCGAGGTTCAACGCCGCGACGCGCTGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTgg                                                                  >  6:200075/1‑148 (MQ=255)
                                                                               cgGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGtt                                                                 <  7:301972/149‑1 (MQ=255)
                                                                               cgGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGtt                                                                 <  7:207423/149‑1 (MQ=255)
                                                                                    aCAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTcag                                                              >  1:153362/1‑147 (MQ=255)
                                                                                     ctggTCTGCGCGCTGAATGCGGTTTTGCGCACGGTAGGCTTCGAGGTTCAACGCCGCGACGGGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTcagcag                                                           <  8:272610/147‑1 (MQ=255)
                                                                                       ggTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCg                                                         <  7:248738/149‑1 (MQ=255)
                                                                                        gTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCgg                                                        >  8:27093/1‑149 (MQ=255)
                                                                                        gTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCg                                                         >  1:234932/1‑148 (MQ=255)
                                                                                        gTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCg                                                         <  2:234932/148‑1 (MQ=255)
                                                                                         tCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCgg                                                        >  2:130754/1‑148 (MQ=255)
                                                                                           ggCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCaa                                                     <  6:290388/149‑1 (MQ=255)
                                                                                                gcTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGgcgc                                                                           >  4:98565/1‑122 (MQ=255)
                                                                                                gcTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGgcgc                                                                           <  3:98565/122‑1 (MQ=255)
                                                                                                  tGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATccgc                                              >  4:49770/1‑149 (MQ=255)
                                                                                                    aaTGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATccgcc                                             >  7:103429/1‑148 (MQ=255)
                                                                                                    aaTGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATccgc                                              <  8:91930/147‑1 (MQ=255)
                                                                                                     aTGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGa                                           <  7:27093/149‑1 (MQ=255)
                                                                                                        cGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGAtt                                         <  8:103429/148‑1 (MQ=255)
                                                                                                        cGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGAtt                                         <  1:169652/148‑1 (MQ=255)
                                                                                                         ggTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGc                                       >  6:43245/1‑149 (MQ=255)
                                                                                                               gcgcACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCgg                                  <  5:43245/148‑1 (MQ=255)
                                                                                                                cgcACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGc                                 >  8:66094/1‑148 (MQ=255)
                                                                                                                 gcACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGa                               <  4:264768/149‑1 (MQ=255)
                                                                                                                   aCGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAg                              >  3:237516/1‑148 (MQ=255)
                                                                                                                      gTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAg                          >  3:214731/1‑149 (MQ=255)
                                                                                                                                aGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTg                <  1:126988/149‑1 (MQ=255)
                                                                                                                                   ccaataaccgcGACCTGCGCGTCGCGGTTATAGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGGGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCgg             >  1:35170/8‑149 (MQ=255)
                                                                                                                                      ataaccgcGACCTGCGCGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCGGCCa          >  8:54169/5‑149 (MQ=255)
                                                                                                                                         gccgcGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCGGCCACTg       >  3:214476/1‑149 (MQ=255)
                                                                                                                                          ccgcGACCTGCGCGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTAGGCGAGTAAGCAGGGCCTTGCGGCCACTgc      >  6:148999/1‑149 (MQ=255)
                                                                                                                                              gACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCGGCCACTGCGCAg  >  6:314126/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGGTGATGCCAGATTCGCCCGTTTGCGACATGTTCGCCGGGACACGCTGGCGGCTGCCGCTGCCGGTGGCCGAAACCGCCGGGAACAGGTCGGCGCGCTGAATGCGGTATTGCGCACGGTAGGCGTCGAGGTTCAACGCCGCGACGCGCAGGTCGCGGTTATTGACCAGCGAGGTCTGGATCAGCTGTTGCAGCGCCGGGTCGTGGAAGAACTGGCGCCAGCCCTGTTCAGCAGCGGCAACATCCGCCGATTGCGTCGGCGAGTAAGCAGGGCCTTGCGGCCACTGCGCAG  >  NC_002947/1577863‑1578153

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