Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,815,531 C→T 7.8% intergenic (+15/+76) PP_5090 → / ← PP_5091 cell division protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,815,5310CT 7.8% 186.0 / 4.4 64intergenic (+15/+76)PP_5090/PP_5091cell division protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base T (3/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TCAAGGATGAAACCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCAC  >  NC_002947/5815383‑5815675
                                                                                                                                                    |                                                                                                                                                
tCAAGGATGAAACCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTAGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGc                                                                                                                                                  <  5:366771/149‑1 (MQ=255)
         aaaCCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGa                                                                                                                                          <  1:220193/148‑1 (MQ=255)
         aaaCCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGAt                                                                                                                                         <  7:386794/149‑1 (MQ=255)
          aaCCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATg                                                                                                                                        >  3:59382/1‑149 (MQ=255)
           aCCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGc                                                                                                                                       >  1:339414/1‑149 (MQ=255)
            ccTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGc                                                                                                                                       >  5:317971/1‑148 (MQ=255)
             cTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGa                                                                                                                                     <  6:134704/149‑1 (MQ=255)
                 tACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACACAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                 <  1:335001/149‑1 (MQ=255)
                   ccGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCAAg                                                                                                                                 <  4:37999/147‑1 (MQ=255)
                   ccGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCAAg                                                                                                                                 >  3:37999/1‑147 (MQ=255)
                         gCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTCGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCAGCAGCGATGCGACAAGGCcttctt                                                                                                                         >  1:372449/1‑149 (MQ=255)
                             gTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtata                                                                                                                     <  2:339414/149‑1 (MQ=255)
                                   ccGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCTGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACa                                                                                                               <  3:313615/149‑1 (MQ=255)
                                     gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                 <  5:205516/129‑1 (MQ=255)
                                     gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                 >  6:205516/1‑129 (MQ=255)
                                     gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                             <  5:325906/149‑1 (MQ=255)
                                         aGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTg                                                                                                         >  5:229021/1‑149 (MQ=255)
                                          gCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGc                                                                                                        >  6:159399/1‑149 (MQ=255)
                                              ccGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGt                                                                                                    <  6:75210/149‑1 (MQ=255)
                                                 cgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                             <  8:325585/137‑1 (MQ=255)
                                                 cgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                             >  7:325585/1‑137 (MQ=255)
                                                  gAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                <  2:183232/149‑1 (MQ=255)
                                                   aGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                <  3:152472/148‑1 (MQ=255)
                                                   aGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                >  4:152472/1‑148 (MQ=255)
                                                       gTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCt                                                                                           <  4:59382/149‑1 (MQ=255)
                                                             aGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAag                                                                                     <  6:317971/149‑1 (MQ=255)
                                                                   acaGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGccc                                                                               >  5:213309/1‑149 (MQ=255)
                                                                          cagcTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCAcc                                                                        >  5:80103/1‑149 (MQ=255)
                                                                                gctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgc                                                                  <  6:80103/149‑1 (MQ=255)
                                                                                      gCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccgcc                                                            <  6:229021/149‑1 (MQ=255)
                                                                                      gCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccgcc                                                            >  6:305342/1‑149 (MQ=255)
                                                                                         gggTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGccccgccccg                                                         >  2:207483/1‑149 (MQ=255)
                                                                                          ggTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGccccgccccg                                                         >  2:354664/1‑148 (MQ=255)
                                                                                                 cAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCtt                                                 >  1:358150/1‑149 (MQ=255)
                                                                                                  aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                >  8:83533/1‑149 (MQ=255)
                                                                                                  aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                <  1:367301/149‑1 (MQ=255)
                                                                                                   cgcTGTTGCAGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTc                                                    <  5:204705/142‑1 (MQ=255)
                                                                                                   cgcTGTTGCAGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTc                                                    >  6:204705/3‑144 (MQ=255)
                                                                                                     cTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCg                                              >  5:426231/1‑148 (MQ=255)
                                                                                                        gtgctgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGa                                            <  4:110546/146‑1 (MQ=255)
                                                                                                           ctgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACAcc                                        <  5:305342/148‑1 (MQ=255)
                                                                                                              cAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                     <  1:207483/148‑1 (MQ=255)
                                                                                                               aaCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGa                                                                                                  <  2:38396/86‑1 (MQ=255)
                                                                                                               aaCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGa                                                                                                  >  1:38396/1‑86 (MQ=255)
                                                                                                                    cGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTc                               >  8:305267/1‑148 (MQ=255)
                                                                                                                       cAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACg                           >  5:60521/1‑149 (MQ=255)
                                                                                                                          aCCCGTCAGTAACGAAGAAGGCCTTGCAGCAGCGATGCCACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCACCGCCCCGCATTCCTTTCGGACCCCTGGCTGTTCAACGTc                         >  4:342868/1‑148 (MQ=255)
                                                                                                                           cccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTcc                        <  6:213309/148‑1 (MQ=255)
                                                                                                                           cccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                         >  7:300149/1‑147 (MQ=255)
                                                                                                                            ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                         >  5:16685/1‑146 (MQ=255)
                                                                                                                            ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCAGTGAAGGGGCTGTAAAGCAGCCCCGCCCCCGTCAGCCCCGCCCCGCATTACTTTCGGACACCTGGCTGTTCAACGta                         >  5:77523/1‑145 (MQ=255)
                                                                                                                                     aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCCCCGTCAGCCCCGCCCCGCATTCATTTCGGCCAACAGGCGGTTCAAAGGCCAGAAGGCATaa             >  8:385707/1‑148 (MQ=255)
                                                                                                                                     aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCCGCCCCCCCCCGCCTTCCTTTCGGACACCTGGCTGTTCCACGTCCAGAAGTCCTAc             >  3:26589/1‑149 (MQ=255)
                                                                                                                                     aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATAc             <  2:358150/149‑1 (MQ=255)
                                                                                                                                      cGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACa            <  6:426231/149‑1 (MQ=255)
                                                                                                                                         agaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCa         >  1:6083/1‑149 (MQ=255)
                                                                                                                                          gaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCACGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                <  5:257346/109‑1 (MQ=255)
                                                                                                                                          gaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCACGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                >  6:257346/1‑109 (MQ=255)
                                                                                                                                            agGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCAcc       <  7:305267/148‑1 (MQ=255)
                                                                                                                                               ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                         >  3:399739/1‑127 (MQ=255)
                                                                                                                                               ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                         <  4:399739/127‑1 (MQ=255)
                                                                                                                                               ccTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCcca   <  6:16685/149‑1 (MQ=255)
                                                                                                                                                 ttGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGACCAGAAGTCATACAGCACCCccac  >  7:111884/1‑148 (MQ=255)
                                                                                                                                                 ttGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                     >  6:61680/1‑113 (MQ=255)
                                                                                                                                                 ttGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                     <  5:61680/113‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAAGGATGAAACCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCAC  >  NC_002947/5815383‑5815675

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