Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,906,511 A→T 7.6% L111Q (CTG→CAG)  PP_4315 ← epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,906,5110AT 7.6% 233.9 / 4.1 75L111Q (CTG→CAG) PP_4315epimerase
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base T (3/3);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTGAAC  >  NC_002947/4906365‑4906660
                                                                                                                                                  |                                                                                                                                                     
tttCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGc                                                                                                                                                    >  2:508624/1‑150 (MQ=255)
    gAGCAGCATGGCGTCATCGATCACCCCGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtg                                                                                                                                                <  2:1244077/150‑2 (MQ=255)
    gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATc                                                                                                                                                >  1:224151/1‑150 (MQ=255)
    gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTAGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATc                                                                                                                                                >  2:286869/1‑150 (MQ=255)
     agcagcATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                               <  2:224151/150‑1 (MQ=255)
       cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTTCGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtggac                                                                                                                                             <  2:1033777/150‑5 (MQ=255)
           aTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                              <  2:905048/145‑4 (MQ=255)
           aTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                              >  1:905048/1‑142 (MQ=255)
              gCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                      <  1:508624/150‑1 (MQ=255)
                 tCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTg                                                                                                                                   >  2:417586/1‑150 (MQ=255)
                      gATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgcg                                                                                                                              >  1:245745/1‑150 (MQ=255)
                          accacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                      >  1:1367499/1‑138 (MQ=255)
                          accacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                      <  2:1367499/138‑1 (MQ=255)
                           ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGACACGGCTGGCGCGCAACt                                                                                                                         >  2:1170202/1‑150 (MQ=255)
                            cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                        >  2:1104497/1‑150 (MQ=255)
                                  aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                               >  2:140900/1‑121 (MQ=255)
                                  aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                               <  1:140900/121‑1 (MQ=255)
                                   tAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTg                                                                                                                 <  2:1336611/150‑1 (MQ=255)
                                    aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGACGTGCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                                <  1:1170202/150‑1 (MQ=255)
                                    aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                               >  1:1318441/1‑135 (MQ=255)
                                    aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                               <  2:1318441/135‑1 (MQ=255)
                                      tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                             <  2:522990/135‑1 (MQ=255)
                                      tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                             >  1:522990/1‑135 (MQ=255)
                                      tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                              <  1:1104497/150‑1 (MQ=255)
                                            cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                        >  2:776197/1‑150 (MQ=255)
                                            cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                        <  1:108300/150‑1 (MQ=255)
                                             tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                             <  2:1411684/128‑1 (MQ=255)
                                             tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                             >  1:1411684/1‑128 (MQ=255)
                                                 gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCa                                                                                                                                                  >  2:310495/1‑103 (MQ=255)
                                                 gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCa                                                                                                                                                  <  1:310495/103‑1 (MQ=255)
                                                      ccGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCa                                                                                              <  1:776197/150‑1 (MQ=255)
                                                       cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCAc                                                                                             >  1:688340/1‑150 (MQ=255)
                                                         ccTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGTAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACtt                                                                                           >  1:987505/1‑150 (MQ=255)
                                                          cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                          <  1:1045140/134‑1 (MQ=255)
                                                          cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                          >  2:1045140/1‑134 (MQ=255)
                                                           ttCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgc                                                                                         <  1:402114/150‑1 (MQ=255)
                                                              gCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                      >  1:296504/1‑150 (MQ=255)
                                                                     tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                               >  2:359043/1‑150 (MQ=255)
                                                                     tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                               >  2:627877/1‑150 (MQ=255)
                                                                     tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGATTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                               >  1:1389947/1‑150 (MQ=255)
                                                                         gCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAg                                                                           >  1:511534/1‑150 (MQ=255)
                                                                            cAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTg                                                                        >  2:844879/1‑150 (MQ=255)
                                                                                 cttgcctaaGTCCCGCTGTGCTGCCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCACCCCCACGTCTGACGCGCAACTCGCCACTCCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGa                                                                   <  2:559934/143‑1 (MQ=255)
                                                                                   tGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                                 >  1:775828/1‑150 (MQ=255)
                                                                                   tGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                                 >  1:1281273/1‑150 (MQ=255)
                                                                                    gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACt                                                                                            >  2:403941/1‑122 (MQ=255)
                                                                                    gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACt                                                                                            <  1:403941/122‑1 (MQ=255)
                                                                                    gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACa                                                                >  1:1121031/1‑150 (MQ=255)
                                                                                            cccGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTaa                                                        <  2:775828/150‑1 (MQ=255)
                                                                                              cGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGc                                                      >  2:524468/1‑150 (MQ=255)
                                                                                                cTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTAAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCg                                                    >  2:75274/1‑150 (MQ=255)
                                                                                                 tGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGa                                                   <  2:245745/150‑1 (MQ=255)
                                                                                                   gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAgg                                                 <  1:417586/150‑1 (MQ=255)
                                                                                                    cATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGc                                                <  1:844879/150‑1 (MQ=255)
                                                                                                     aTGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGcc                                               >  1:222485/1‑150 (MQ=255)
                                                                                                      tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCa                                                             >  1:326277/1‑135 (MQ=255)
                                                                                                      tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCa                                                             <  2:326277/135‑1 (MQ=255)
                                                                                                            ccGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                             <  2:1409225/145‑1 (MQ=255)
                                                                                                            ccGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                             >  1:1409225/1‑145 (MQ=255)
                                                                                                              gCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCa                                      <  2:688340/150‑1 (MQ=255)
                                                                                                               cTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCAt                                     <  1:524468/150‑1 (MQ=255)
                                                                                                                tACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATg                                    >  1:1208176/1‑150 (MQ=255)
                                                                                                                tACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATg                                    <  1:855439/150‑1 (MQ=255)
                                                                                                                     gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGc                               >  2:545424/1‑150 (MQ=255)
                                                                                                                      ttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGca                              <  2:1281273/150‑1 (MQ=255)
                                                                                                                              cGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGa                      >  1:1300169/1‑150 (MQ=255)
                                                                                                                                   aaGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTcc                 >  1:1120865/1‑150 (MQ=255)
                                                                                                                                    aGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCa                <  1:545424/150‑1 (MQ=255)
                                                                                                                                     gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc               >  1:1278147/5‑150 (MQ=255)
                                                                                                                                        cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt            <  2:1389947/150‑1 (MQ=255)
                                                                                                                                        cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt            >  2:1311322/1‑150 (MQ=255)
                                                                                                                                         aTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg           >  1:140368/1‑150 (MQ=255)
                                                                                                                                         aTGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg           <  2:296504/150‑1 (MQ=255)
                                                                                                                                          tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg          <  1:359043/150‑1 (MQ=255)
                                                                                                                                          tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg          >  1:640878/1‑150 (MQ=255)
                                                                                                                                              cagcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGt      <  1:404336/150‑1 (MQ=255)
                                                                                                                                                  agGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTGAAc  <  2:1121031/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGTGAAC  >  NC_002947/4906365‑4906660

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