Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,220,867 G→A 5.6% Q463Q (CAG→CAA)  PP_3699 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,220,8670GA 5.6% 213.1 / 2.6 71Q463Q (CAG→CAA) PP_3699hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base A (2/2);  total (35/36)
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.53e-01

GAGGTGGTCAGCCTCACCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTGGGC  >  NC_002947/4220719‑4221015
                                                                                                                                                    |                                                                                                                                                    
gAGGTGGTCAGCCTCACCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTTGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGt                                                                                                                                                     <  2:413183/150‑1 (MQ=255)
             tCACCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGAt                                                                                                                                        >  2:152852/1‑150 (MQ=255)
             tCACCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAATACCCTGGGAAGAt                                                                                                                                        >  1:127352/1‑150 (MQ=255)
              gaCCAAGATGACGTCGTGAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATg                                                                                                                                       <  2:512481/149‑1 (MQ=255)
               aCCAAGATGACGTCGTTAAACACGGTCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATgg                                                                                                                                      <  2:1241883/150‑1 (MQ=255)
               aCCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATgg                                                                                                                                      >  1:86978/1‑150 (MQ=255)
                 cAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCa                                                                                                                                    <  2:1283109/150‑1 (MQ=255)
                 cAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAATACCCTGGGAAGATGGCa                                                                                                                                    <  2:127352/150‑1 (MQ=255)
                      tGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGGGAAGACGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATgat                                                                                                                               >  1:1157383/1‑150 (MQ=255)
                         cgtcgtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATgatcga                                                                                                                            <  2:918278/150‑1 (MQ=255)
                          gtcgtAAAACACGACCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTATGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACTGTAGTCAATCTACTTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCCATgatcgat                                                                                                                           >  2:636915/1‑150 (MQ=255)
                           tcgtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATTGCAATgatcgatc                                                                                                                          <  1:874004/150‑1 (MQ=255)
                            cgtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGatcgatcg                                                                                                                         <  1:152852/150‑1 (MQ=255)
                             gtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGc                                                                                                                        >  1:255146/1‑150 (MQ=255)
                              ttAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCa                                                                                                                       >  2:912788/1‑150 (MQ=255)
                                         cAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCCGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAActggggctgg                                                                                                            >  1:307572/1‑150 (MQ=255)
                                          aaGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAActggggctggg                                                                                                           >  1:1402094/1‑150 (MQ=255)
                                               cgcATTTATGAGCGCAATAGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCTTTAGCACGGTAGTCAATCTACCTGCTGATCGCCACAGCCAGTACCCTGGGAAGATGGCAATGATCGATCTCAACTGGGGCTGGGGAAGc                                                                                                      <  2:560250/150‑1 (MQ=255)
                                                 cATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGccc                                                                                                    <  1:1049035/150‑1 (MQ=255)
                                                   tttATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCgg                                                                                                  >  1:735821/1‑150 (MQ=255)
                                                    ttATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGt                                                                                                 <  1:912788/150‑1 (MQ=255)
                                                       tGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACCGGGGCTGGGGAAGCCCGGTtga                                                                                              <  2:1157383/150‑1 (MQ=255)
                                                        gACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTtgat                                                                                             <  1:1037517/150‑1 (MQ=255)
                                                        gACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTtgat                                                                                             <  1:837095/150‑1 (MQ=255)
                                                          ccGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTtgatga                                                                                           >  1:1023388/1‑150 (MQ=255)
                                                            gCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATgaga                                                                                         <  2:1402094/150‑1 (MQ=255)
                                                            gCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATgaga                                                                                         <  1:423466/150‑1 (MQ=255)
                                                              aaTGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGAtt                                                                                       >  1:638989/1‑150 (MQ=255)
                                                                       cTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGa                                                                                                                                         <  1:294486/91‑1 (MQ=255)
                                                                       cTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGa                                                                                                                                         >  2:294486/1‑91 (MQ=255)
                                                                              cTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGaa                                                                                                        <  2:939444/117‑1 (MQ=255)
                                                                              cTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGaa                                                                                                        >  1:939444/1‑117 (MQ=255)
                                                                                   ttGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATgg                                                                                                                                      >  1:1139734/1‑82 (MQ=255)
                                                                                   ttGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATgg                                                                                                                                      <  2:1139734/82‑1 (MQ=255)
                                                                                         tgGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTcc                                                            <  2:735821/150‑1 (MQ=255)
                                                                                         tgGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAATACCCTGGGAAGATGGCAATgatc                                                                                                                              <  1:615275/84‑1 (MQ=255)
                                                                                         tgGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAATACCCTGGGAAGATGGCAATgatc                                                                                                                              >  2:615275/1‑84 (MQ=255)
                                                                                         tgGAATGCAAGAGTGGTTGTTGGCATGGTAGGCATTGTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGTAGAGCCTAGCCGCTcc                                                            <  2:307572/150‑1 (MQ=255)
                                                                                             aTGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcggcg                                                        >  2:1413309/1‑150 (MQ=255)
                                                                                              tGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcggcgg                                                       >  1:863592/1‑150 (MQ=255)
                                                                                               gCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGa                                                      <  2:638989/150‑1 (MQ=255)
                                                                                                      tGGTAGTTACCACGGTCATCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGtt                                               >  2:205856/1‑150 (MQ=255)
                                                                                                       ggTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTa                                              <  2:255146/150‑1 (MQ=255)
                                                                                                        gTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTAc                                             >  1:917368/1‑150 (MQ=255)
                                                                                                        gTCGTTAGCACGGTAGGCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTAc                                             >  2:950122/1‑150 (MQ=255)
                                                                                                         tCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAAt                                                                                                                                  <  2:830135/64‑1 (MQ=255)
                                                                                                         tCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAAt                                                                                                                                  >  1:830135/1‑64 (MQ=255)
                                                                                                         tCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTcc                                                                                     >  1:1419815/1‑109 (MQ=255)
                                                                                                         tCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTcc                                                                                     <  2:1419815/109‑1 (MQ=255)
                                                                                                         tCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACg                                            >  2:384655/1‑150 (MQ=255)
                                                                                                            ttAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTtgat                                                                                             >  1:1136831/1‑98 (MQ=255)
                                                                                                            ttAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTtgat                                                                                             <  2:1136831/98‑1 (MQ=255)
                                                                                                            ttAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTa                                         >  2:933341/1‑150 (MQ=255)
                                                                                                                 aCGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGagag                                                                        >  2:339253/1‑114 (MQ=255)
                                                                                                                 aCGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGagag                                                                        <  1:339253/114‑1 (MQ=255)
                                                                                                                  cGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATTGATCGCAACGGGGGCTGGGGACGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGAcc                                   <  2:86978/150‑1 (MQ=255)
                                                                                                                            tCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGc                                           <  2:418896/132‑1 (MQ=255)
                                                                                                                            tCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGc                                           >  1:418896/1‑132 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCg                                                                             <  1:547123/94‑1 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCg                                                                             >  2:547123/1‑94 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAAtt                        >  1:491720/1‑147 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAAtt                        <  2:491720/147‑1 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTa                     <  2:917368/150‑1 (MQ=255)
                                                                                                                                ccTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTACAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTa                     <  2:1023388/150‑1 (MQ=255)
                                                                                                                                 ctgctgATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAg                    >  2:1113249/1‑150 (MQ=255)
                                                                                                                                       aTCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTa                     >  1:1166103/1‑143 (MQ=255)
                                                                                                                                       aTCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTa                     <  2:1166103/143‑1 (MQ=255)
                                                                                                                                              cAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCAcc       >  2:100679/1‑150 (MQ=255)
                                                                                                                                                   aGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAggtg                                                  >  1:361623/1‑102 (MQ=255)
                                                                                                                                                   aGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAggtg                                                  <  2:361623/102‑1 (MQ=255)
                                                                                                                                                   aGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTGGGc  >  1:1165069/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGGTGGTCAGCCTCACCAAGATGACGTCGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTGGGC  >  NC_002947/4220719‑4221015

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