Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,954 G→A 9.2% I167I (ATC→ATT PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,9540GA 9.2% 176.3 / 7.0 65I167I (ATC→ATTPP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (3/3);  total (34/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 = 2.80e-01

GCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCC  >  NC_002947/4051813‑4052101
                                                                                                                                              |                                                                                                                                                   
gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAgg                                                                                                                                               >  6:217870/1‑148 (MQ=255)
gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAg                                                                                                                                                <  3:134260/147‑1 (MQ=255)
 ccAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGc                                                                                                                                              <  6:166101/148‑1 (MQ=255)
  cAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGGTACGCCGAATGGCACGCTGCGCACTCCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGCGTTCACCGATCAGGca                                                                                                                                             >  1:249541/1‑149 (MQ=255)
   aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGca                                                                                                                                             >  1:44263/1‑147 (MQ=255)
                    aCTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGttt                                                                                                                                 >  5:271804/1‑142 (MQ=255)
                    aCTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGttt                                                                                                                                 <  6:271804/142‑1 (MQ=255)
                        gATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCAGAGCCGTGGTTTACACGGgc                                                                                                                         <  6:154063/146‑1 (MQ=255)
                        gATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGgc                                                                                                                         >  5:154063/1‑146 (MQ=255)
                              gTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                                                                                <  1:265434/149‑1 (MQ=255)
                                   ccccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGGATTGCCAGCTTGCCCTTGTTCGCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGTCGCTGTTGGAgggga                                                                                                           <  2:249541/149‑3 (MQ=255)
                                     ccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagc                                                                                                         <  4:141539/149‑1 (MQ=255)
                                       aTACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagc                                                                                                         >  3:68059/1‑147 (MQ=255)
                                        tACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagccag                                                                                                      >  2:245794/1‑149 (MQ=255)
                                              gcgTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATc                                                                                                  >  2:108213/1‑147 (MQ=255)
                                                     gCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCAc                                                                                          <  5:37755/149‑1 (MQ=255)
                                                     gCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCAc                                                                                          >  6:37755/1‑149 (MQ=255)
                                                       cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGc                                                                                        >  6:218592/1‑149 (MQ=255)
                                                               aCGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAaacaa                                                                                >  7:190735/1‑149 (MQ=255)
                                                                 gCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAAcc                                                                              <  1:245794/149‑1 (MQ=255)
                                                                   tGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                                                                            >  2:41893/1‑149 (MQ=255)
                                                                        ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGActtct                                                                       >  4:212564/1‑149 (MQ=255)
                                                                          tGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACttcttc                                                                     >  4:31793/1‑149 (MQ=255)
                                                                            cACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACttcttc                                                                     <  6:142496/147‑1 (MQ=255)
                                                                            cACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACttcttc                                                                     <  5:217870/147‑1 (MQ=255)
                                                                             aCCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGt                                                                   >  8:206534/1‑148 (MQ=255)
                                                                               cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa                                                                 >  6:162187/1‑148 (MQ=255)
                                                                               cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa                                                                 >  7:316434/1‑148 (MQ=255)
                                                                                         aTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGc                                                      <  5:218592/149‑1 (MQ=255)
                                                                                              cAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 >  6:188937/1‑149 (MQ=255)
                                                                                              cAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 <  8:316434/149‑1 (MQ=255)
                                                                                                     cccTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCg                                           <  8:190735/148‑1 (MQ=255)
                                                                                                       cTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCa                                        >  2:149735/1‑149 (MQ=255)
                                                                                                             ccccAGAACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCAGgcgc                                   >  6:302531/1‑148 (MQ=255)
                                                                                                              cccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCaa                                 <  5:188937/149‑1 (MQ=255)
                                                                                                              cccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCaa                                 >  6:329063/1‑149 (MQ=255)
                                                                                                              cccAGTACCGGTTTCAGTTCGGTGGGTTCACCGATCAGGCACAGCCCTGGTTTACACGGGCGCTCTTGGAGGGCAGACAGCATCGAGCGCACGCCCAAACAACCGACTTCTACGACATAGGAAAACGCCAGGTGTACCGGCATGCGCaa                                 >  6:84575/1‑149 (MQ=255)
                                                                                                               ccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCatc                                <  5:84575/149‑3 (MQ=255)
                                                                                                               ccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAAc                                >  5:83788/1‑149 (MQ=255)
                                                                                                               ccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAAc                                >  8:280521/1‑149 (MQ=255)
                                                                                                               ccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAAc                                <  5:162187/149‑1 (MQ=255)
                                                                                                               ccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAAc                                <  1:149735/149‑1 (MQ=255)
                                                                                                                     ccGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGc                                         <  8:275407/134‑1 (MQ=255)
                                                                                                                     ccGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGc                                         >  7:275407/1‑134 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGTGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTgcgc                        <  1:108213/149‑1 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 >  7:188399/1‑124 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 >  7:193385/1‑124 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 <  8:193385/124‑1 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 <  8:188399/124‑1 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTGACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 >  7:193387/1‑124 (MQ=255)
                                                                                                                       ggTTTGAGTTCGGTGGGTTCACCAATCAGGCACAGCCGTGGTTAACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgt                                                 <  8:193387/124‑1 (MQ=255)
                                                                                                                         tttGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCga                      <  7:280521/149‑1 (MQ=255)
                                                                                                                              gTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCAt                                                                >  2:184081/1‑102 (MQ=255)
                                                                                                                              gTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCAt                                                                <  1:184081/102‑1 (MQ=255)
                                                                                                                              gTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg                                                              >  5:54658/1‑104 (MQ=255)
                                                                                                                              gTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg                                                              <  6:54658/104‑1 (MQ=255)
                                                                                                                              gTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAAc                 >  8:271198/1‑149 (MQ=255)
                                                                                                                                  ggTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAAcgcg              <  3:212564/148‑1 (MQ=255)
                                                                                                                                      ggTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGca          <  6:83788/148‑1 (MQ=255)
                                                                                                                                          cACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAg                                                   >  6:256401/1‑103 (MQ=255)
                                                                                                                                          cACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAg                                                   <  5:256401/103‑1 (MQ=255)
                                                                                                                                          cACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCAcagc      <  4:68059/148‑1 (MQ=255)
                                                                                                                                           aCCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGcc                                                     <  7:203991/100‑1 (MQ=255)
                                                                                                                                           aCCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGcc                                                     >  8:203991/1‑100 (MQ=255)
                                                                                                                                            ccGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACagcagc   >  2:176033/1‑149 (MQ=255)
                                                                                                                                              gATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGcc  >  3:290721/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCC  >  NC_002947/4051813‑4052101

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: