Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,391,329 A→G 8.4% E82E (GAA→GAG)  PP_5707 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,391,3290AG 8.4% 198.0 / 7.3 71E82E (GAA→GAG) PP_5707hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base G (3/3);  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 = 9.59e-01

CTTGCCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCAA  >  NC_002947/5391181‑5391476
                                                                                                                                                    |                                                                                                                                                    
cTTGCCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGaa                                                                                                                                                      >  2:53399/1‑149 (MQ=255)
 ttGCCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGaa                                                                                                                                                      >  2:99148/1‑148 (MQ=255)
    ccAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTc                                                                                                                                                   <  2:238975/148‑1 (MQ=255)
     cAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCgg                                                                                                                                                 <  3:42574/149‑1 (MQ=255)
      aaaaCGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGc                                                                                                                                                <  1:28332/149‑1 (MQ=255)
       aaaCGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCa                                                                                                                                               <  1:108122/149‑1 (MQ=255)
             ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatcc                                                                                                                                         <  4:28820/149‑3 (MQ=255)
             ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatcc                                                                                                                                         <  4:175928/149‑3 (MQ=255)
             ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGc                                                                                                                                         >  8:32606/1‑149 (MQ=255)
             ccGCGAACAGGCCTCGGAACTGATGCGCTATGCGCGGGACCTGACCATCGACACACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGCGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGc                                                                                                                                         >  1:310846/1‑149 (MQ=255)
               gcgAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGt                                                                                                                                       <  6:381050/149‑1 (MQ=255)
                          tCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCa                                                                                                                            >  1:187386/1‑149 (MQ=255)
                             gAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAACCTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCCCGCACcg                                                                                                                         <  2:310846/149‑1 (MQ=255)
                                  gATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATg                                                                                                                    >  5:133778/1‑149 (MQ=255)
                                   aTGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATg                                                                                                                    <  2:391542/148‑1 (MQ=255)
                                     gggCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCTCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCt                                                                                                                  <  8:58291/148‑1 (MQ=255)
                                        cTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              <  7:52747/149‑1 (MQ=255)
                                           tGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATACGTTGAATCACGCACCGCGATGCTTAAAggg                                                                                                           <  6:203126/149‑1 (MQ=255)
                                            gcgcGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAgggg                                                                                                          >  8:286547/1‑149 (MQ=255)
                                               cgGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCa                                                                                                       >  4:240883/1‑149 (MQ=255)
                                                    ccTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGt                                                                                                   >  6:76271/1‑148 (MQ=255)
                                                        aCCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAg                                                                                                <  6:96310/147‑1 (MQ=255)
                                                           aTCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              >  2:91014/1‑130 (MQ=255)
                                                           aTCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              <  1:91014/130‑1 (MQ=255)
                                                                 ggACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGccc                                                                                     <  8:194052/149‑1 (MQ=255)
                                                                    cTGATGGAAGACAAGCCCAAACTGATCTGGGCGACGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCatgat                                                                                                                                           >  2:299928/1‑92 (MQ=255)
                                                                       aTGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCtttt                                                                                <  7:286547/148‑1 (MQ=255)
                                                                         ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGAt                                                                                                                     >  1:20183/1‑109 (MQ=255)
                                                                         ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGAt                                                                                                                     <  2:20183/109‑1 (MQ=255)
                                                                         ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             >  8:217562/1‑149 (MQ=255)
                                                                          gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGTATGATGCCGAACACGGAATGATGCGTTGAATCACCCACCGCGATGCTAAAAGGGGTAATCGACAGCCGATTGGCCCCTTTTAccc                                                                            >  5:344188/1‑147 (MQ=255)
                                                                          gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             <  3:240883/148‑1 (MQ=255)
                                                                          gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             <  6:133778/148‑1 (MQ=255)
                                                                            aGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCca                                                                          >  1:157117/1‑149 (MQ=255)
                                                                                  gCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGAt                                                                   >  2:115872/1‑149 (MQ=255)
                                                                                   cccAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATc                                                                  <  5:76271/149‑1 (MQ=255)
                                                                                       aaCTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACg                                                              >  3:181843/1‑149 (MQ=255)
                                                                                            aTCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGa                                                          <  7:217562/148‑1 (MQ=255)
                                                                                                gggCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     >  1:353537/1‑149 (MQ=255)
                                                                                                gggCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     >  1:347937/1‑149 (MQ=255)
                                                                                                 ggcggcGCACTACCTCTGTTCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGTTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAgg                                                        <  1:99148/145‑1 (MQ=255)
                                                                                                 ggcggcGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     <  7:32606/148‑1 (MQ=255)
                                                                                                 ggcggcGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATg                                                    <  7:61881/149‑1 (MQ=255)
                                                                                                   cggcgCACTACCTCTGTGCTCTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  >  5:142388/1‑149 (MQ=255)
                                                                                                      cgcACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCaa                                                                >  4:382109/1‑132 (MQ=255)
                                                                                                      cgcACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCaa                                                                <  3:382109/132‑1 (MQ=255)
                                                                                                         aCTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  <  4:370542/143‑1 (MQ=255)
                                                                                                         aCTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  >  3:370542/1‑143 (MQ=255)
                                                                                                            aCCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTtct                                          >  3:203766/1‑148 (MQ=255)
                                                                                                             ccTCTGTGCTCTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCa                                        <  6:142388/149‑1 (MQ=255)
                                                                                                              ctctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTtct                                          >  3:250092/1‑146 (MQ=255)
                                                                                                              ctctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAg                                       >  3:44276/1‑149 (MQ=255)
                                                                                                                      tttGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTc                               >  5:128593/1‑149 (MQ=255)
                                                                                                                      tttGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTc                               >  6:157557/1‑149 (MQ=255)
                                                                                                                           gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTc                                                                                                        <  4:98522/72‑1 (MQ=255)
                                                                                                                           gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTc                                                                                                        >  3:98522/1‑72 (MQ=255)
                                                                                                                           gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAtt                          >  5:348026/1‑149 (MQ=255)
                                                                                                                            cAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAttt                         >  6:276844/1‑149 (MQ=255)
                                                                                                                            cAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACCCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAtt                          >  7:211672/1‑148 (MQ=255)
                                                                                                                             aaaGCGCTACTTGATGATGCCTAACTCGTCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTa                        <  6:344188/149‑1 (MQ=255)
                                                                                                                             aaaGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTa                        <  6:128593/149‑1 (MQ=255)
                                                                                                                               aGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGc                                                                                       >  7:153859/1‑85 (MQ=255)
                                                                                                                               aGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGc                                                                                       <  8:153859/85‑1 (MQ=255)
                                                                                                                                gcgcTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATgc                     <  1:53399/149‑1 (MQ=255)
                                                                                                                                 cgcTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATgcg                    >  8:36017/1‑149 (MQ=255)
                                                                                                                                      acgcatcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGcgccg                 >  5:357175/8‑147 (MQ=255)
                                                                                                                                       tGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCgcag              >  5:161280/1‑149 (MQ=255)
                                                                                                                                         gatgatGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACCACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCgcagc             <  1:115872/149‑1 (MQ=255)
                                                                                                                                          atcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCcaca               <  4:288477/142‑1 (MQ=255)
                                                                                                                                          atcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCcaca               >  3:288477/4‑145 (MQ=255)
                                                                                                                                           tgatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          <  2:353537/149‑1 (MQ=255)
                                                                                                                                           tgatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          <  2:347937/149‑1 (MQ=255)
                                                                                                                                           tcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          >  3:380018/3‑149 (MQ=255)
                                                                                                                                                   aaCTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCaa  <  4:203766/149‑1 (MQ=255)
                                                                                                                                                   aaCTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCaa  >  8:48410/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGCCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCAA  >  NC_002947/5391181‑5391476

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