Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 100,113 C→A 9.7% G99V (GGG→GTG) ‡ PP_0094 ← 5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947100,1130CA 9.7% 200.4 / 9.1 78G99V (GGG→GTG) ‡PP_00945'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base A (6/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

TCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAG  >  NC_002947/99970‑100261
                                                                                                                                               |                                                                                                                                                    
tCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGtt                                                                                                                                                  >  1:255894/1‑148 (MQ=255)
   gCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAgg                                                                                                                                              >  6:326758/1‑149 (MQ=255)
   gCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAgg                                                                                                                                              >  1:3830/1‑149 (MQ=255)
    cTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAggg                                                                                                                                             <  2:255894/149‑1 (MQ=255)
    cTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAggg                                                                                                                                             >  4:263188/1‑149 (MQ=255)
          cTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTATGGAGGTGGCAATGTACAGCGTGCGGCCCTGCACGTTCAGGGCCTc                                                                                                                                         >  8:274840/1‑147 (MQ=255)
             gttgGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTTCCCGTTCAGGGCCTCAAGc                                                                                                                                     <  5:157401/148‑1 (MQ=255)
             gttgGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGc                                                                                                                                     <  5:326758/148‑1 (MQ=255)
              gtgGTGCGCGTGCCTTCCCGGTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGGGGCAATGTACAGCGTGCGGCCCTGAACGGGCAGGGCCGCACGct                                                                                                                                    <  6:98031/147‑2 (MQ=255)
                 gTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGcc                                                                                                                                           >  6:352771/1‑138 (MQ=255)
                 gTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGcc                                                                                                                                           <  5:352771/138‑1 (MQ=255)
                 gTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTc                                                                                                                                >  3:149008/1‑149 (MQ=255)
                 gTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTc                                                                                                                                <  8:249484/149‑1 (MQ=255)
                 gTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTc                                                                                                                                >  3:154660/1‑149 (MQ=255)
                     gcgGGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGc                                                                                                                            <  3:42588/149‑1 (MQ=255)
                        tGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGcaca                                                                                                                         >  8:97654/1‑149 (MQ=255)
                           cGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGAACGGGCAGGGccgc                                                                                                                                         >  4:115309/1‑128 (MQ=255)
                           cGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGAACGGGCAGGGccgc                                                                                                                                         <  3:115309/130‑3 (MQ=255)
                            gTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCAGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCTTGCCTCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGCCACACCta                                                                                                                     >  4:49363/1‑148 (MQ=255)
                                    cACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAacac                                                                                                              >  5:92824/1‑148 (MQ=255)
                                     aCTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAacac                                                                                                              >  5:120642/1‑147 (MQ=255)
                                     aCTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGa                                                                                                            <  3:247849/149‑1 (MQ=255)
                                     aCTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGa                                                                                                            <  6:62262/149‑1 (MQ=255)
                                     aCTGCCGTAGATCACCTAGTAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGa                                                                                                            <  8:333260/149‑1 (MQ=255)
                                       tGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAgg                                                                                                          >  6:126397/1‑149 (MQ=255)
                                         ccGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGccc                                                                                                                                                     <  2:269177/104‑1 (MQ=255)
                                         ccGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAgg                                                                                                          <  7:97654/147‑1 (MQ=255)
                                         ccGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCATGccc                                                                                                                                                     >  1:269177/1‑104 (MQ=255)
                                              gATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGAGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCATGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTc                                                                                                     <  1:314171/147‑1 (MQ=255)
                                              gATCACCTTGAAGTGATGATCGAACGCAAAGTGCAGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTAAGGGCCTCAAGCAACTAAGGCACACCTTCGAACACGAGGtttta                                                                                                     >  2:314171/1‑146 (MQ=255)
                                                tCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGGCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTa                                                                                                  <  1:42183/148‑1 (MQ=255)
                                                 cACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTa                                                                                                  >  1:13531/1‑147 (MQ=255)
                                                    cTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGAc                                                                                             >  5:36865/1‑149 (MQ=255)
                                                    cTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGAc                                                                                             >  4:199525/1‑149 (MQ=255)
                                                     ttGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTAAAGCGTGCGGCCCTGCCCCTTCAGGGCCTAACGCAACTCCGGCACACCTTCGAGCACGAGGTTTTCGTAAAGAcc                                                                                            >  5:278931/1‑149 (MQ=255)
                                                       gAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCg                                                                                           <  2:3830/148‑1 (MQ=255)
                                                           tgatgaTCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCAc                                                                                      >  5:56488/1‑149 (MQ=255)
                                                                 tCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGaa                                                                                <  2:208446/149‑1 (MQ=255)
                                                                 tCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAAGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCgg                                                                                  >  2:195615/1‑147 (MQ=255)
                                                                  cGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAg                                                                               >  7:93216/1‑149 (MQ=255)
                                                                      cGCTAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCt                                                                           <  1:288181/149‑1 (MQ=255)
                                                                           aGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCgg                                                                      >  5:180789/1‑149 (MQ=255)
                                                                                 gAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAg                                                                >  8:339219/1‑149 (MQ=255)
                                                                                       tttCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCAcc                                                          <  3:263188/149‑1 (MQ=255)
                                                                                           gcgGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCAt                                                      <  6:56488/149‑1 (MQ=255)
                                                                                               gCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATcc                                                    >  2:278527/1‑147 (MQ=255)
                                                                                                cGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAg                                                 <  7:339219/149‑1 (MQ=255)
                                                                                                    cTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCt                                             >  4:140696/1‑149 (MQ=255)
                                                                                                      cccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTg                                           >  6:69683/1‑149 (MQ=255)
                                                                                                       ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCtt                                            <  5:126397/147‑1 (MQ=255)
                                                                                                       ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                          >  3:115274/1‑149 (MQ=255)
                                                                                                       ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                          >  7:89821/1‑149 (MQ=255)
                                                                                                       ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                          >  8:134014/1‑149 (MQ=255)
                                                                                                         aTGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATACCAGGCCTTGGCt                                        >  6:305296/1‑149 (MQ=255)
                                                                                                          tGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCtt                                       <  6:36865/149‑1 (MQ=255)
                                                                                                             tttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgt                                    <  6:92824/149‑1 (MQ=255)
                                                                                                             tttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgt                                    <  2:13531/149‑1 (MQ=255)
                                                                                                             tttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgt                                    <  4:154660/149‑1 (MQ=255)
                                                                                                             tttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgt                                    <  4:149008/149‑1 (MQ=255)
                                                                                                                 gAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgtcg                                  >  2:337210/1‑147 (MQ=255)
                                                                                                                 gAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGaa                                <  1:278527/149‑1 (MQ=255)
                                                                                                                  aGGTGGCAATGTACAGCGTGCGGCCCTGCGCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAg                               <  7:274840/149‑1 (MQ=255)
                                                                                                                  aGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCTTAAAGACCGGTCACTCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAg                               <  4:115274/149‑1 (MQ=255)
                                                                                                                  aGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAg                               <  1:147105/149‑1 (MQ=255)
                                                                                                                         aaTGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAg                         >  5:63019/1‑148 (MQ=255)
                                                                                                                             tACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgtcg                                  <  4:113727/135‑1 (MQ=255)
                                                                                                                             tACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgtcg                                  >  3:113727/1‑135 (MQ=255)
                                                                                                                               cAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCa                  <  6:180789/149‑1 (MQ=255)
                                                                                                                               cAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCa                  <  8:93216/149‑1 (MQ=255)
                                                                                                                                 gcttgaggCCCTGAACGGGCAGGGCCGCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCAt                 >  6:79001/7‑148 (MQ=255)
                                                                                                                                  cttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATg                >  3:5212/6‑148 (MQ=255)
                                                                                                                                   ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGa               >  7:193297/5‑148 (MQ=255)
                                                                                                                                   ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGa               >  5:14867/5‑148 (MQ=255)
                                                                                                                                    tGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGaaa             <  3:199525/149‑1 (MQ=255)
                                                                                                                                      cGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGc           <  1:337210/149‑1 (MQ=255)
                                                                                                                                      cGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGc           <  6:120642/149‑1 (MQ=255)
                                                                                                                                       ttgggtGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGc           <  3:344938/143‑1 (MQ=255)
                                                                                                                                       ggCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTCCATGAAAGcc          >  5:2379/1‑149 (MQ=255)
                                                                                                                                          ccTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTgc       <  5:69683/149‑1 (MQ=255)
                                                                                                                                           cTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTg        >  3:254378/1‑147 (MQ=255)
                                                                                                                                               ccGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGcagcag  >  7:224195/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAG  >  NC_002947/99970‑100261

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