Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,532 A→C 5.7% intergenic (+9/+1666) PP_3326 → / ← PP_3328 hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,5320AC 5.7% 200.1 / 4.3 70intergenic (+9/+1666)PP_3326/PP_3328hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base C (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCG  >  NC_002947/3763386‑3763668
                                                                                                                                                  |                                                                                                                                        
gCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACAt                                                                                                                                       >  1:1302234/1‑150 (MQ=255)
    gTATTAGCGTACGCAATGAAGGACTCAACATTATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCGCGTGGCCCAGCAGAGCGTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACAACGAGTAAACTAAGCGACATACCt                                                                                                                                   >  1:475581/1‑150 (MQ=255)
     tATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTc                                                                                                                                  >  2:1363628/1‑150 (MQ=255)
            gTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATg                                                                                                                           >  1:870957/1‑150 (MQ=255)
             tACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                                                                                                          >  2:339932/1‑150 (MQ=255)
                gCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACgg                                                                                                                       >  1:209903/1‑150 (MQ=255)
                    tGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                   <  2:704534/150‑1 (MQ=255)
                     gAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAg                                                                                                                  <  1:360007/150‑1 (MQ=255)
                           cTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATgg                                                                                                            <  1:256942/150‑1 (MQ=255)
                            tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGc                                                                                                           >  1:503807/1‑150 (MQ=255)
                              aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAg                                                                                                         >  2:735886/1‑150 (MQ=255)
                               aCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGt                                                                                                        <  2:1168651/150‑1 (MQ=255)
                               aCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGt                                                                                                        <  1:849569/150‑1 (MQ=255)
                                  tcatcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTgg                                                                                                     >  2:199562/1‑150 (MQ=255)
                                   catcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGa                                                                                                    <  1:1363628/150‑1 (MQ=255)
                                   catcGCGGCCAAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGa                                                                                                    <  2:456797/150‑1 (MQ=255)
                                    atcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                   >  2:1000513/1‑150 (MQ=255)
                                         ggCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGtt                                                                                              <  2:1393178/150‑1 (MQ=255)
                                          gCCCAGGCTACCCTGCCAATGGCAGCCGGCCGACATGTGGCTCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTc                                                                                             <  2:1302234/150‑1 (MQ=255)
                                           cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                   <  2:1116728/127‑1 (MQ=255)
                                           cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                   >  1:1116728/1‑127 (MQ=255)
                                           cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                   >  1:484414/1‑143 (MQ=255)
                                           cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                   <  2:484414/143‑1 (MQ=255)
                                              aGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGaa                                                                                         <  1:199562/150‑1 (MQ=255)
                                                 cTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGAACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGa                                                                                      <  1:98952/150‑1 (MQ=255)
                                                  tACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAg                                                                                     >  2:834369/1‑150 (MQ=255)
                                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGAc                                                                                                                                         <  1:396502/95‑1 (MQ=255)
                                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGAc                                                                                                                                         >  2:396502/1‑95 (MQ=255)
                                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCa                                                                                                          >  2:515631/1‑126 (MQ=255)
                                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCa                                                                                                          <  1:515631/126‑1 (MQ=255)
                                                      cTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGcc                                                                                 >  1:588207/1‑150 (MQ=255)
                                                         ccAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGc                                                                              >  2:431422/1‑150 (MQ=255)
                                                           aaTGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                   >  1:184062/1‑111 (MQ=255)
                                                           aaTGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                   <  2:184062/111‑1 (MQ=255)
                                                                                   ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                <  1:557109/122‑1 (MQ=255)
                                                                                   ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                >  2:557109/1‑122 (MQ=255)
                                                                                   ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTg                                                    >  2:1395872/1‑150 (MQ=255)
                                                                                    cagcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTgg                                                   <  1:339932/150‑1 (MQ=255)
                                                                                     agcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGCCATACCTCGGCCGTCATGGCCGAg                                                                                                                  <  1:45316/86‑1 (MQ=37)
                                                                                     agcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGCCATACCTCGGCCGTCATGGCCGAg                                                                                                                  >  2:45184/1‑86 (MQ=37)
                                                                                     agcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGCCATACCTCGGCCGTCATGGCCGAg                                                                                                                  <  1:45184/86‑1 (MQ=38)
                                                                                     agcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGCCATACCTCGGCCGTCATGGCCGAg                                                                                                                  >  2:45316/1‑86 (MQ=37)
                                                                                      gcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAg                                                                                                                  <  2:1305172/85‑1 (MQ=38)
                                                                                      gcagAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAg                                                                                                                  >  1:1305172/1‑85 (MQ=38)
                                                                                            ccTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGcc                                           <  2:772320/150‑1 (MQ=255)
                                                                                             cTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGccc                                          <  1:735886/150‑1 (MQ=255)
                                                                                                   agttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGc                                    >  1:1067489/1‑150 (MQ=255)
                                                                                                    gttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCg                                   <  1:1000513/150‑1 (MQ=255)
                                                                                                     ttcaCTGTCGTTCCTGGGTAAGCCGCACAACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGCCACAAGCCCGGCTGCgg                                  <  2:475581/150‑1 (MQ=255)
                                                                                                          tGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt                                                                                    >  2:179546/1‑95 (MQ=255)
                                                                                                          tGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt                                                                                    <  1:179546/95‑1 (MQ=255)
                                                                                                              gTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGACTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGc                         >  2:1112944/1‑150 (MQ=255)
                                                                                                                  cTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTGGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACgg                          >  1:1357704/1‑145 (MQ=255)
                                                                                                                  cTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGGGTTCGGAATGAGTGCCTGCAACCa                                                                         >  2:809100/1‑98 (MQ=255)
                                                                                                                  cTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCGCGGCTGCGGCCACACGGCAATc                     >  1:1009707/1‑150 (MQ=255)
                                                                                                                  cTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACgg                          <  2:1357704/145‑1 (MQ=255)
                                                                                                                  cTGGGTAAGCCGCACCACGAGTAAACCAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGGGTTCGGAATGAGTGCCTGCAACCa                                                                         <  1:809100/98‑1 (MQ=255)
                                                                                                                     ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGaa                  <  2:1067489/150‑1 (MQ=255)
                                                                                                                       tAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGCGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAaca                >  2:165504/1‑150 (MQ=255)
                                                                                                                       tAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAaca                <  1:1276358/150‑1 (MQ=255)
                                                                                                                           ccGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGAt            >  1:1336338/1‑150 (MQ=255)
                                                                                                                            cGCACGACGAGTAAAGTGAGCGACATACCTCGGCAATGACGGCCGAGGTATGGCAGTTTGAGGCGTTCGGAATGAGTGCCCGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCCCACACGATg           <  2:492365/150‑1 (MQ=255)
                                                                                                                             gCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATCGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGAATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATgg          <  2:209903/150‑1 (MQ=255)
                                                                                                                                ccacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCa       <  1:431422/150‑1 (MQ=255)
                                                                                                                                 cacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCAt      >  2:737532/1‑150 (MQ=255)
                                                                                                                                   cGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGt    <  1:834369/150‑1 (MQ=255)
                                                                                                                                    gAGTAAAATAAAAGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTc   >  2:770199/1‑150 (MQ=255)
                                                                                                                                     aGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCAt      >  2:250023/1‑146 (MQ=255)
                                                                                                                                     aGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCAt      <  1:250023/146‑1 (MQ=255)
                                                                                                                                     aGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCg  >  2:1418478/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GCACGTATTACCGTACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCG  >  NC_002947/3763386‑3763668

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