Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,767,114 G→T 5.4% T168T (ACG→ACT)  PP_3329 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,767,1140GT 5.4% 242.8 / 2.2 74T168T (ACG→ACT) PP_3329membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base T (2/2);  total (37/37)
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.06e-01

ATCATTGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGTGAAAAGCC  >  NC_002947/3766971‑3767238
                                                                                                                                               |                                                                                                                            
atcatTGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGAc                                                                                                                        >  2:598854/1‑150 (MQ=255)
 tcatTGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGAcc                                                                                                                       <  2:1160928/150‑1 (MQ=255)
     tGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCgtgt                                                                                                                   <  1:120343/150‑1 (MQ=255)
     tGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCgtgt                                                                                                                   <  1:265770/150‑1 (MQ=255)
       cACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCg                                                                                                                 >  1:669496/1‑150 (MQ=255)
                ccGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAg                                                                                                        >  2:1001564/1‑150 (MQ=255)
                      ccTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCa                                                                                                  >  2:1129373/1‑150 (MQ=255)
                        tACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCaaa                                                                                                <  2:846642/150‑1 (MQ=255)
                            gTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACTATCGACCgtgt                                                                                                                   <  1:890578/127‑1 (MQ=255)
                            gTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACTATCGACCgtgt                                                                                                                   >  2:890637/1‑127 (MQ=255)
                            gTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACTATCGACCgtgt                                                                                                                   <  1:890637/127‑1 (MQ=255)
                            gTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACTATCGACCgtgt                                                                                                                   >  2:890578/1‑127 (MQ=255)
                                acaAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGa                                                                                        >  2:1213516/1‑150 (MQ=255)
                                acaAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGa                                                                                        >  2:838725/1‑150 (MQ=255)
                                     gTCGTGTAGTGCTTACCCTCGTCACGATCGCCGTGCTCGTCCTGATCCCTTACCTCCTGTTGTTCTGGATTGATCACCTCTCCAAACTTGTTCTCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGcgccgc                                                                                   <  1:1471239/150‑1 (MQ=255)
                                      tCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCg                                                                                  >  1:1449605/1‑150 (MQ=255)
                                        gtgtAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCgtg                                                                                <  1:838725/150‑1 (MQ=255)
                                          gtAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGtt                                                                              >  2:108601/1‑150 (MQ=255)
                                           tAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTg                                                                             >  1:526716/1‑150 (MQ=255)
                                             gTGCTTTTCACCGTGACGATTGCATTGGTGGTGCTGATCATCTACACCTTGGTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGaa                                                                           >  2:493556/1‑150 (MQ=255)
                                             gTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGaa                                                                           >  1:1146873/1‑150 (MQ=255)
                                              tGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAac                                                                          >  2:289046/1‑150 (MQ=255)
                                              tGCTTTTCACCGTGACGAATGCAGTGGTGGTCCTGATCATCTACACCTTGCGGCTCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAac                                                                          >  2:1296260/1‑150 (MQ=255)
                                                    tCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTAGACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAACAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGt                                                                    >  1:1331347/1‑150 (MQ=255)
                                                    tCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGt                                                                    <  1:1001564/150‑1 (MQ=255)
                                                          tGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTgg                                                              <  1:1213516/150‑1 (MQ=255)
                                                                 tGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGc                                                       <  2:1374718/150‑1 (MQ=255)
                                                                         tggtCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGAGCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTAcc                                               >  2:985362/1‑150 (MQ=255)
                                                                            tCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTg                                                               <  2:1624560/131‑1 (MQ=255)
                                                                            tCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTg                                                               >  1:1624560/1‑131 (MQ=255)
                                                                            tCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGc                                            <  1:598854/150‑1 (MQ=255)
                                                                            tCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGc                                            <  2:1449605/150‑1 (MQ=255)
                                                                              cTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCt                                          >  1:79556/1‑150 (MQ=255)
                                                                                gATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGa                                        <  2:669496/150‑1 (MQ=255)
                                                                                 atcatcTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGAc                                       >  1:1239278/1‑150 (MQ=255)
                                                                                      cTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGTTCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCaa                                  <  2:863864/150‑1 (MQ=255)
                                                                                      cTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCaa                                  <  1:1129373/150‑1 (MQ=255)
                                                                                      cTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCaa                                  <  1:79542/150‑1 (MQ=255)
                                                                                      cTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGCGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCaa                                  <  1:379723/150‑1 (MQ=255)
                                                                                       tACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAg                                 >  1:452884/1‑150 (MQ=255)
                                                                                            cTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATc                                                                                          <  2:1548519/88‑1 (MQ=255)
                                                                                            cTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATc                                                                                          >  1:1548519/1‑88 (MQ=255)
                                                                                            cTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTcc                            >  1:553216/1‑150 (MQ=255)
                                                                                              tgctgcGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTg                          <  2:526716/150‑1 (MQ=255)
                                                                                                      tGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGcc                                           <  1:369618/125‑1 (MQ=255)
                                                                                                      tGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGcc                                           >  2:369618/1‑125 (MQ=255)
                                                                                                      tGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAAcc                  >  1:1647794/1‑150 (MQ=255)
                                                                                                        gATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGATTCCTGGGGTAACCCa                >  1:933840/1‑150 (MQ=255)
                                                                                                          ttGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATg              <  1:493556/150‑1 (MQ=255)
                                                                                                            gATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGaa            <  2:1146873/150‑1 (MQ=255)
                                                                                                             aTCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGaa            <  1:1314573/149‑1 (MQ=255)
                                                                                                             aTCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGaa            >  2:1314573/1‑149 (MQ=255)
                                                                                                              tCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGt          <  2:1239278/150‑1 (MQ=255)
                                                                                                              tCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGt          >  1:761245/1‑150 (MQ=255)
                                                                                                              tCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGt          <  1:985362/150‑1 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGtt                                                                              >  1:950571/1‑80 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGtt                                                                              <  2:950571/80‑1 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTg                          >  2:191224/1‑132 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTg                          <  1:191224/132‑1 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCAt               <  2:318958/143‑1 (MQ=255)
                                                                                                                aCCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCAt               >  1:318958/1‑143 (MQ=255)
                                                                                                                 ccTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTgggg                       >  1:948771/1‑134 (MQ=255)
                                                                                                                 ccTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTgggg                       <  2:948771/134‑1 (MQ=255)
                                                                                                                   tctcCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGTGaaaa     <  1:108601/150‑1 (MQ=255)
                                                                                                                   tctcCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGTGaaaa     <  1:1296260/150‑1 (MQ=255)
                                                                                                                    ctcCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTaa                    >  1:1446135/1‑134 (MQ=255)
                                                                                                                    ctcCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTaa                    <  2:1446135/134‑1 (MQ=255)
                                                                                                                      ccAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGTGAAAAgcc  >  1:1312760/1‑150 (MQ=255)
                                                                                                                         aaCTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCt                                          >  1:1272789/1‑107 (MQ=255)
                                                                                                                         aaCTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCt                                          <  2:1272789/107‑1 (MQ=255)
                                                                                                                               ggCGCTTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCgt                                                                                 >  2:1231797/1‑62 (MQ=255)
                                                                                                                               ggCGCTTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCgt                                                                                 <  1:1231797/62‑1 (MQ=255)
                                                                                                                                             aCGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGaa            >  2:812223/1‑117 (MQ=255)
                                                                                                                                             aCGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGaa            <  1:812223/117‑1 (MQ=255)
                                                                                                                                               gATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGaa                                                                           <  1:1430733/52‑1 (MQ=255)
                                                                                                                                               gATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGaa                                                                           >  2:1430733/1‑52 (MQ=255)
                                                                                                                                               |                                                                                                                            
ATCATTGCACTGAGCACCGGAGCCTACGGTACACAAGGTCGTGTAGTGCTTTTCACCGTGACGATTGCAGTGGTGGTCCTGATCATCTACACCTTGCTGCGCTGGATTGATCACCTCTCCAAACTAGGGCGCGTAGGCGAGACGATCGACCGTGTCGAAAAAGCAGTGATCAAAGCCATCGAGCAGCGTGTTGAACACCCGTTTCTGGGCGCTGCAGCCTACCCGCCTGACCTCAAGCTTCCTGGGGTAACCCATGAAGTGAAAAGCC  >  NC_002947/3766971‑3767238

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: ‑