Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,437 G→A 6.2% H109Y (CAC→TAC)  PP_1324 ← endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,509,4370GA 6.2% 227.1 / 3.4 80H109Y (CAC→TAC) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base A (3/2);  total (40/40)
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.87e-01

ATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAC  >  NC_002947/1509288‑1509585
                                                                                                                                                     |                                                                                                                                                    
aTTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      >  2:1137456/1‑150 (MQ=255)
aTTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      <  1:193365/150‑1 (MQ=255)
      tGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                <  2:2076843/150‑1 (MQ=255)
          cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                            <  1:1038465/150‑1 (MQ=255)
          cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                            >  2:1038465/1‑150 (MQ=255)
          cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTCTTCAGCCATTGAACCGGTCTGCCTGCATGCTGGCTAcccccc                                                                                                                                            >  1:174362/1‑147 (MQ=255)
            tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                >  2:566961/1‑144 (MQ=255)
            tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                <  1:566961/144‑1 (MQ=255)
             cTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAgg                                                                                                                                         <  1:2236109/150‑1 (MQ=255)
                gTGACCTTAAATGGGGCGGTGGCCTTTGCTACGACCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAAAGCAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCCTTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCg                                                                                                                                      >  1:1939534/1‑150 (MQ=255)
                 tGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGa                                                                                                                                     <  2:1903177/150‑1 (MQ=255)
                  gACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                    >  2:1818610/1‑150 (MQ=255)
                   aCCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                   >  2:2048629/1‑150 (MQ=255)
                   aCCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                   >  1:480479/1‑150 (MQ=255)
                      ttAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACg                                                                                                                                >  1:1277658/1‑150 (MQ=255)
                        aaaTGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                              <  2:2096848/150‑1 (MQ=255)
                            ggggCGGTGGCCTTTGCTACTCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCCAAg                                                                                                                          >  2:514603/1‑150 (MQ=255)
                             gggCGGTGGCCTTTGCTAGGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGc                                                                                                                         <  2:1590044/150‑1 (MQ=255)
                                    ggCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAggg                                                                                                                  <  2:480479/150‑1 (MQ=255)
                                     gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAGGCGGCAgggg                                                                                                                 >  2:806085/1‑150 (MQ=255)
                                     gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgggg                                                                                                                 <  2:1277658/150‑1 (MQ=255)
                                        tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      <  2:925404/110‑1 (MQ=255)
                                        tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      >  1:925404/1‑110 (MQ=255)
                                               cgcgCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTTAAGCGGTCTGCCTGCATGGTGGCTACGCTGCAGGGCGACCACGTCGGAGCGGCAGAGGTAGTTGcccc                                                                                                       >  2:529150/1‑150 (MQ=255)
                                                gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                            <  1:1537609/144‑1 (MQ=255)
                                                gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                            >  2:1537609/1‑144 (MQ=255)
                                                   caacaaCGACCAGCCCGCAATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTg                                                                                                   <  1:264413/150‑1 (MQ=255)
                                                    aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                  <  1:2048629/150‑1 (MQ=255)
                                                    aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                  <  2:1993838/150‑1 (MQ=255)
                                                    aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                  <  1:434255/150‑1 (MQ=255)
                                                     acaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTggg                                                                                                 >  2:1319995/1‑150 (MQ=255)
                                                      caacGACCAGCCCTCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                <  1:529150/150‑1 (MQ=255)
                                                       aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                               >  2:235813/1‑150 (MQ=255)
                                                          gACCAGCCCGCGATACATAGAGCAAAAAAAACGAATTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcctcc                                                                                            >  1:1436291/1‑150 (MQ=255)
                                                          gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                   <  1:1351468/111‑1 (MQ=255)
                                                          gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                   >  2:1351468/1‑111 (MQ=255)
                                                              aGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCt                                                                                        >  1:1612884/1‑150 (MQ=255)
                                                                cccGCGATACATAGAGCAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    <  2:1836497/119‑1 (MQ=255)
                                                                cccGCGATACATAGAGCAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    >  1:1836497/1‑119 (MQ=255)
                                                                 ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                <  1:803490/139‑1 (MQ=255)
                                                                 ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                >  2:803490/1‑139 (MQ=255)
                                                                 ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCa                                                                                     <  1:1512605/150‑1 (MQ=255)
                                                                        aCATAGAGCAAAAAAAACGAATTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAg                                                                              <  2:1436291/150‑1 (MQ=255)
                                                                        aCATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAg                                                                              >  1:1771330/1‑150 (MQ=255)
                                                                         cATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCGCAGTGGGCGTCCTTCTGAAGGaaaggg                                                                             >  2:916328/1‑146 (MQ=255)
                                                                         cATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                             <  2:1771330/150‑1 (MQ=255)
                                                                             gagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTGTGTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                         <  2:187071/150‑1 (MQ=255)
                                                                                cAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                          >  1:467683/1‑82 (MQ=38)
                                                                                cAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                          <  2:467683/82‑1 (MQ=38)
                                                                                 aaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcggc                                                                     >  2:1059041/1‑150 (MQ=255)
                                                                                 aaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGCTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcggc                                                                     <  1:1319995/150‑1 (MQ=255)
                                                                                   aaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                   <  1:128694/150‑1 (MQ=255)
                                                                                       aaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTCGCTGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCAc                                                               <  1:614726/150‑1 (MQ=255)
                                                                                        aCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCAcc                                                              <  2:1612884/150‑1 (MQ=255)
                                                                                          gAAGTGGGTTCAGCATTCGAAGGCGTTTTTTTTTTTTTTTTTTGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCaa                                                            >  2:976992/1‑150 (MQ=255)
                                                                                          gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCGGCAGGAACAGGGTGGCGGCAAGCACCac                                                            >  2:1775269/1‑149 (MQ=255)
                                                                                           aaGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAAc                                                           >  1:1869100/1‑150 (MQ=255)
                                                                                             gTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCg                                                         >  2:340275/1‑150 (MQ=255)
                                                                                                          tCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTgcgc                                            >  2:1974842/1‑150 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      >  1:852483/1‑42 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      <  2:852483/42‑1 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                <  1:1848163/96‑1 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                >  2:1848163/1‑96 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                         <  1:257859/119‑1 (MQ=255)
                                                                                                            gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                         >  2:257859/1‑119 (MQ=255)
                                                                                                             aaGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                             >  2:690624/1‑114 (MQ=255)
                                                                                                             aaGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                             <  1:690624/114‑1 (MQ=255)
                                                                                                                 cGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGa                                     <  1:1137456/150‑1 (MQ=255)
                                                                                                                  gTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                                    >  1:1298787/1‑150 (MQ=255)
                                                                                                                     tttCAGCCATTGAAGCGGTCGGCCATCATGATGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACt                                 <  1:514603/150‑1 (MQ=255)
                                                                                                                        cAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                              <  2:1869100/150‑1 (MQ=255)
                                                                                                                          gCCATTGAAGCTGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                            <  1:916328/150‑1 (MQ=255)
                                                                                                                          gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCATGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                            >  1:2071449/1‑150 (MQ=255)
                                                                                                                           ccATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAGGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                           <  1:806085/150‑1 (MQ=255)
                                                                                                                                gAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGAGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACAAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGCGAgcgc                      >  2:2230647/1‑150 (MQ=255)
                                                                                                                                  aGCGGTCTGCCTGCATGGTGGCTACCCTGCATGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCt                    <  2:2071449/150‑1 (MQ=255)
                                                                                                                                    cGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTcc                  <  1:1818610/150‑1 (MQ=255)
                                                                                                                                      gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGa                >  1:290503/1‑150 (MQ=255)
                                                                                                                                      gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGCAGATGCGTTGGTGAGCGCCTCCGa                >  2:540962/1‑150 (MQ=255)
                                                                                                                                       tCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGaa               <  1:976992/150‑1 (MQ=255)
                                                                                                                                           ccttcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGccc           <  1:588313/150‑1 (MQ=255)
                                                                                                                                            ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcg          >  2:1424517/1‑150 (MQ=255)
                                                                                                                                            ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                         <  2:1412188/71‑1 (MQ=37)
                                                                                                                                            ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                         >  1:1412188/1‑71 (MQ=37)
                                                                                                                                            ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGTGCGCCTCCGAAGcccc          >  2:1045913/1‑149 (MQ=255)
                                                                                                                                                aTGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTg      >  1:1284778/1‑150 (MQ=255)
                                                                                                                                                 tggtggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGccc                                          <  2:868753/113‑1 (MQ=255)
                                                                                                                                                 tggtggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGGGGCGGCAAGCACCAACCGTTTCTGCTTGCGccc                                          >  1:868753/1‑113 (MQ=255)
                                                                                                                                                    tggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAc  <  2:174362/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAC  >  NC_002947/1509288‑1509585

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