Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,973 C→T 8.3% R302Q (CGG→CAG)  PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,9730CT 8.3% 170.9 / 3.4 60R302Q (CGG→CAG) PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

GCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAA  >  NC_002947/3865832‑3866121
                                                                                                                                              |                                                                                                                                                    
gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                               <  2:337413/149‑1 (MQ=255)
gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                               <  4:52553/149‑1 (MQ=255)
gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                               >  4:123210/1‑149 (MQ=255)
 cgcgAGAAGGCCTTGGCGGGCACTCTCAGGTTCTCTCGCGCACCATCAGCTCAAAGACCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCACCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtggg                                                                                                                                              <  4:142677/149‑3 (MQ=255)
        aGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATgg                                                                                                                                       <  8:29843/149‑1 (MQ=255)
          gacatGGCGGGCACTCTCAGGTGCTCTCGCCCACCAACACCACAAAGCCCAGATCCTGCTGCGCATTCGTCACCAGCTTCCCATCCAGCAACGCCACCACCGCCTGTGCCGCCCCCCGCCCCACCGCATCCCGCCGTGTGCGAATGGAg                                                                                                                                     >  4:311929/5‑149 (MQ=255)
                  gggCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCg                                                                                                                             >  7:264284/1‑149 (MQ=255)
                       ctctCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCa                                                                                                                        <  5:95157/149‑1 (MQ=255)
                        tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                                       >  5:341176/1‑149 (MQ=255)
                        tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                                       >  5:73213/1‑149 (MQ=255)
                         ctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                                       >  3:130795/1‑148 (MQ=255)
                            aGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAcc                                                                                                                      <  3:311929/146‑1 (MQ=255)
                                  tctcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc                                                                                                             >  2:1715/1‑149 (MQ=255)
                                  tctcGCGCACCATCAGCTCAAAGCCAACAACCTGCTGCGCAATCGTCACCCGCTTCCCATCCAGCAACGCAAGCAGCCCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATCCGcc                                                                                                             >  5:271765/1‑149 (MQ=255)
                                      gcgcACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAgg                                                                                                         >  8:207777/1‑149 (MQ=255)
                                         cACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCg                                                                                                      >  1:166469/1‑149 (MQ=255)
                                          aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                     <  3:123210/149‑1 (MQ=255)
                                            cATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                     >  1:143206/1‑147 (MQ=255)
                                              tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGgaga                                                                                                 <  7:207777/149‑1 (MQ=255)
                                                   tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGCTCGtt                                                                                            >  5:273411/1‑149 (MQ=255)
                                                   tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                            >  8:116834/1‑149 (MQ=255)
                                                   tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                             <  3:350511/99‑4 (MQ=255)
                                                   tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                             >  4:350511/1‑96 (MQ=255)
                                                    cAAAGCCCAGATCCTGCTGCGCATTCGGCACCCGCTTCCCTTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGTTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                            <  6:271765/148‑1 (MQ=255)
                                                     aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTTCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTg                                                                                           <  1:276732/148‑1 (MQ=255)
                                                     aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCTCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                            <  7:314236/147‑1 (MQ=255)
                                                     aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                             <  7:77074/98‑4 (MQ=255)
                                                     aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                             >  8:77074/1‑95 (MQ=255)
                                                     aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTg                                                                                           >  5:147809/1‑148 (MQ=255)
                                                       aGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaaa                                                                                        >  8:176709/1‑149 (MQ=255)
                                                        gCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaaa                                                                                        >  1:96571/1‑149 (MQ=255)
                                                                         cATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTgg                                                                      >  8:284187/1‑149 (MQ=255)
                                                                          aTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGc                                                                     >  6:67254/1‑149 (MQ=255)
                                                                          aTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGc                                                                     <  6:273411/149‑1 (MQ=255)
                                                                               tCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCACAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                     <  1:165205/112‑1 (MQ=255)
                                                                               tCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCACAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                     >  2:165205/1‑112 (MQ=255)
                                                                                  cccGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCAcc                                                             >  5:248398/1‑149 (MQ=255)
                                                                                    cGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTc                                                           <  8:264284/149‑1 (MQ=255)
                                                                                          ccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                                     <  1:205186/149‑1 (MQ=255)
                                                                                          ccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                                     >  2:205186/1‑149 (MQ=255)
                                                                                           cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                                                                                                    >  2:265108/1‑101 (MQ=255)
                                                                                           cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                                                                                                    <  1:265108/101‑1 (MQ=255)
                                                                                           cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCt                                                    >  5:46363/1‑149 (MQ=255)
                                                                                            aTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCtt                                                   >  8:5665/1‑149 (MQ=255)
                                                                                             tCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTg                                                  <  4:130795/149‑1 (MQ=255)
                                                                                              ccAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGa                                                 >  6:9799/1‑149 (MQ=255)
                                                                                               cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGAc                                                <  7:284187/149‑1 (MQ=255)
                                                                                               cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGAc                                                <  1:1715/149‑1 (MQ=255)
                                                                                                 gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                                     >  4:35331/1‑142 (MQ=255)
                                                                                                 gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                                     <  3:35331/142‑1 (MQ=255)
                                                                                                    aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                           >  3:286918/1‑149 (MQ=255)
                                                                                                    aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                           <  4:286918/149‑1 (MQ=255)
                                                                                                          gcagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgc                                     <  6:73213/149‑1 (MQ=255)
                                                                                                           cagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTg                                      >  1:5026/1‑147 (MQ=255)
                                                                                                                gcCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCAc                                <  7:116834/149‑1 (MQ=255)
                                                                                                                       ccGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCAGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGccc                         >  8:334777/1‑149 (MQ=255)
                                                                                                                          cgccccgcccCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGa                      <  5:67254/149‑1 (MQ=255)
                                                                                                                               cgccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAg                 <  6:46363/149‑1 (MQ=255)
                                                                                                                                gccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGACCTggctcagg                <  2:166469/149‑8 (MQ=255)
                                                                                                                                     aCCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtcgtc           <  6:248398/149‑1 (MQ=255)
                                                                                                                                              cgcgGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACagaa  <  2:96571/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAA  >  NC_002947/3865832‑3866121

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: