Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,130,281 C→G 7.4% L19V (CTG→GTG) ‡ PP_2747 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,130,2810CG 7.4% 204.1 / 6.6 68L19V (CTG→GTG) ‡PP_2747hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/3);  total (32/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 = 8.57e-01

CACCCTGGCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCA  >  NC_002947/3130132‑3130428
                                                                                                                                                     |                                                                                                                                                    
caccCTGGCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGccc                                                                                                                                                      >  1:519445/1‑150 (MQ=255)
 accCTGGCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCCTACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCt                                                                                                                                                     <  1:534929/150‑1 (MQ=255)
  cccTGGCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCAGACATGAAGGCTATCCGCCGTTGCCCGACATAGTGGCCGTCAGGGCGCTg                                                                                                                                                    <  2:422122/150‑1 (MQ=255)
       gcgcCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGgc                                                                                                                                               <  1:882261/150‑1 (MQ=255)
         gcCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGgcgc                                                                                                                                             <  1:1441602/150‑1 (MQ=255)
              cTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtcg                                                                                                                                        <  1:623715/150‑4 (MQ=255)
              cTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCaa                                                                                                                                        <  2:285207/150‑1 (MQ=255)
                 tCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGCCATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGt                                                                                                                                     >  2:629057/1‑150 (MQ=255)
                   gcgACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   <  1:1673803/150‑1 (MQ=255)
                    cgACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCCTGCCCATACATGAATGCTATCACAAGGCGGCTCTTGCAGTGGGCGACAGCGCCCTGACGGCGCTGCAACGTTgg                                                                                                                                  <  1:1790967/150‑1 (MQ=255)
                          cagccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCaga                                                                                                                            <  2:674728/150‑1 (MQ=255)
                            gccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAgag‑a                                                                                                                         >  1:981378/1‑150 (MQ=255)
                             ccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCACAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑aa                                                                                                                        >  2:752684/1‑150 (MQ=255)
                                             gTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCACTCCAGTACGGAGTGTCCGACAGCGCCCTGGCGGAGGCGCAAAGATCGG‑CAGAGCAACGCCTGACACAAATCg                                                                                                        >  1:35060/1‑150 (MQ=255)
                                               cAGTGCCCGCCCGGCAGTGCCGGGCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      <  1:429543/150‑1 (MQ=255)
                                                       gCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCt                                                                                                                                                     >  2:551227/1‑96 (MQ=255)
                                                       gCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCACGACATAGTGGCCGACAGCGCCCt                                                                                                                                                     <  1:551227/96‑1 (MQ=255)
                                                            gCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGcc                                                                                         <  2:438806/150‑1 (MQ=255)
                                                                    ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         >  1:735007/1‑93 (MQ=255)
                                                                    ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         >  1:736507/1‑93 (MQ=255)
                                                                    ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         <  2:735007/95‑3 (MQ=255)
                                                                    ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         <  2:736507/95‑3 (MQ=255)
                                                                    ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGt                                                                                 <  2:1394060/150‑1 (MQ=255)
                                                                      ccTGAATCGGCAGAGCAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCGCACTCAGGGTGTGACAGCGCCATGACGGCCCGGCAACGTTGGGCAAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGc                                                                               <  1:629057/150‑1 (MQ=255)
                                                                       cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACTGCGCGGGAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCg                                                                              >  1:1204290/1‑150 (MQ=255)
                                                                       cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATc                                                                                            >  2:1360537/1‑136 (MQ=255)
                                                                       cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATc                                                                                            <  1:1360537/136‑1 (MQ=255)
                                                                               gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCg                                                                                                        >  1:803946/1‑116 (MQ=255)
                                                                               gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCg                                                                                                        <  2:803946/116‑1 (MQ=255)
                                                                               gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGc                                                                      >  1:1724785/1‑150 (MQ=255)
                                                                               gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGc                                                                      >  2:35892/1‑150 (MQ=255)
                                                                                cAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGcc                                                                     >  1:1272282/1‑150 (MQ=255)
                                                                                  gaggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTg                                                                   >  2:364520/1‑150 (MQ=255)
                                                                                    ggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   <  1:168701/85‑1 (MQ=255)
                                                                                    ggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   >  2:168701/1‑85 (MQ=255)
                                                                                     gagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGAACCAGAAGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGgcg                                                                >  1:1504955/1‑150 (MQ=255)
                                                                                        tCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTg                                                             <  2:1063988/150‑1 (MQ=255)
                                                                                          aTCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGa                                                           <  1:844601/150‑1 (MQ=255)
                                                                                          aTCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGa                                                           <  1:687099/150‑1 (MQ=255)
                                                                                              cATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGa                                                       >  1:762485/1‑150 (MQ=255)
                                                                                               aTGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGAt                                                      >  2:1041803/1‑150 (MQ=255)
                                                                                               aTGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGAt                                                      >  1:1035247/1‑150 (MQ=255)
                                                                                                tGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGCGCGTGGAAGGCCTGGCGCTGGATCGAtg                                                     <  1:1810140/150‑1 (MQ=255)
                                                                                                     aTACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGtgctg                                                >  1:1104865/1‑150 (MQ=255)
                                                                                                        cATGAATGCCATCCGCCGTTGGCCGACATAGTGGACGAAAGCGCTCTGACGGCGCTGCAACGTTGGGCCAGAG‑ACCGCCCGACACAACTCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGt                                             <  2:705889/150‑1 (MQ=255)
                                                                                                         aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGaca                                                                                                               <  2:794032/83‑1 (MQ=255)
                                                                                                         aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGaca                                                                                                               >  1:794032/1‑83 (MQ=255)
                                                                                                         aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTg                                            >  2:1250698/1‑150 (MQ=255)
                                                                                                         aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTg                                            >  1:1266932/1‑150 (MQ=255)
                                                                                                            aaTGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAAAGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTGATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAgg                                         >  2:489405/1‑150 (MQ=255)
                                                                                                              tGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGTGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTc                                       >  2:1376992/1‑150 (MQ=255)
                                                                                                                   tCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCt                                                                                                                                                     >  2:1157077/1‑36 (MQ=255)
                                                                                                                   tCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCt                                                                                                                                                     <  1:1157077/36‑1 (MQ=255)
                                                                                                                   tCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAAGCGCCCTGCGTCATCGCCGTCGCGGTGCGTGGAATAGCGGTCGCTGGTGCGAGGTGCTGCGTGTGGTCGAGCg                                  >  2:1859945/1‑150 (MQ=255)
                                                                                                                    ccgccgTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 <  2:1134368/150‑1 (MQ=255)
                                                                                                                     cgccgTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCggg                                >  1:392621/1‑150 (MQ=255)
                                                                                                                         gTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGc                            <  1:1250698/150‑1 (MQ=255)
                                                                                                                         gTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGc                            <  2:762485/150‑1 (MQ=255)
                                                                                                                             cccGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGAt                        >  2:1155084/1‑150 (MQ=255)
                                                                                                                              ccGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATg                       <  2:1104865/150‑1 (MQ=255)
                                                                                                                               cGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      <  1:919487/70‑1 (MQ=255)
                                                                                                                               cGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      >  2:919487/1‑70 (MQ=255)
                                                                                                                                 aCATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACCATGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCt                    <  1:1859945/150‑1 (MQ=255)
                                                                                                                                   acaGTGGCCGACCGCGCCCTGACGGCGCTGCAACGTTGGCCCAGAG‑AACGCCCGACACAAATCCCCCTGCGTGATCGCCGTCTCGGTGCGTGGCAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGCGTTGCCGATGGCTTg                  <  2:193840/148‑1 (MQ=255)
                                                                                                                                   aTAGCGGCCGACAGCGCCCTGGAGGCGCTGCAACGTTGGGCCAGTG‑AACGCCCGACAAAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTg                  <  2:192650/150‑1 (MQ=255)
                                                                                                                                          ccGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCcagc           >  2:1206338/1‑150 (MQ=255)
                                                                                                                                               aGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 <  1:923513/123‑1 (MQ=255)
                                                                                                                                               aGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 >  2:923513/1‑123 (MQ=255)
                                                                                                                                                  gcCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTc   <  2:1724785/150‑1 (MQ=255)
                                                                                                                                                   cccTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCa  <  2:35060/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACCCTGGCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCA  >  NC_002947/3130132‑3130428

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