Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,308,600 C→G 5.5% A137P (GCC→CCC)  fucA ← L‑fuculose phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,308,6000CG 5.5% 238.5 / 6.2 73A137P (GCC→CCC) fucAL‑fuculose phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (2/2);  total (38/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.42e-01

CAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGGCCATGCCATGGCCTTGC  >  ACB122/2308469‑2308733
                                                                                                                                   |                                                                                                                                     
cAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAg                                                                                                                                    <  5:245440/135‑1 (MQ=255)
 aGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGc                                                                                                                                   <  4:58370/135‑1 (MQ=255)
       gcgcgcTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACt                                                                                                                             >  3:330175/1‑135 (MQ=255)
          cgcTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTcc                                                                                                                          <  4:330175/135‑1 (MQ=255)
          cgcTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTcc                                                                                                                          <  4:191433/135‑1 (MQ=255)
          cgcTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTcc                                                                                                                          <  4:179276/135‑1 (MQ=255)
               gTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCAt                                                                                                                     <  7:373717/135‑1 (MQ=255)
                tCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATc                                                                                                                    >  5:77123/1‑135 (MQ=255)
                tCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTAAGGAATGCACAGTCGTCGATCAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATc                                                                                                                    >  3:289882/1‑135 (MQ=255)
                     aaTCGCGTTCGCGATGATTTCGCCTTCTTATTTGCCCACAGGTACACCTGGAAATTAGTTGAGGAATGCACATTCGTAGAAGAGGGGGCACAGGTAAATGTGCGAGATTGCCAGCGTAACTTACAGCATCTACAAg                                                                                                              <  8:315007/136‑1 (MQ=255)
                         gcgTTCGCGATGATTTCGCCTTCTTCGGTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAggcg                                                                                                           <  8:107897/135‑1 (MQ=255)
                          cgTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAggcg                                                                                                           <  4:210978/134‑1 (MQ=255)
                          cgTTCGCGATGATTTCGCCTTCTTCGTGGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAggcgg                                                                                                          <  6:390171/135‑1 (MQ=255)
                            ttCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGgcggcg                                                                                                        <  6:273175/135‑1 (MQ=255)
                               gcgATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCtg                                                                                                     >  3:178397/1‑135 (MQ=255)
                                 gatgatTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCtgtg                                                                                                   <  2:334553/135‑1 (MQ=255)
                                      tttCGCCTTCTTCGTTGCCGACCGGTACAGCGGGCCATTGCTTGAGGAATGCAAAGTCGTCTAACTGGTGTCACAGTTCCATGTGCGATATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGAAGc                                                                                               <  6:244535/134‑1 (MQ=255)
                                          gCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCgggtg                                                                                          <  7:35884/135‑1 (MQ=255)
                                            cttcttCGTTGCCCACCGGTCCACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCgggtggg                                                                                        <  6:391708/135‑1 (MQ=255)
                                                    ttGCCCACCGGGAGTCCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACTGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGAGAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgatg                                                                                <  6:359302/135‑1 (MQ=255)
                                                        ccACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGc                                                                            <  6:67212/135‑1 (MQ=255)
                                                        ccACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGc                                                                            <  5:22991/135‑1 (MQ=255)
                                                         cACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCa                                                                           >  2:243813/1‑135 (MQ=255)
                                                                       gCCATTCCTTGAGGAATGCACAGTCGTCGAACATGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggc                                                             >  1:272545/1‑135 (MQ=255)
                                                                       gCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggc                                                             >  5:11128/1‑135 (MQ=255)
                                                                       gCCATTCCTTGAGGAATGCACAGTCGGCGCACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggc                                                             <  1:191254/135‑1 (MQ=255)
                                                                         cATTCCTTGAGGAATGGACAGTCGTCGCACAGGGGGCACACGTCCAAGTGCGAGATGGCCAGCGGAACTTCAAGCATAGACAAGGCGGCGCTCTGAAGCGCGTGGGTGTGGATGATGCAGATGACGTCGgggcgg                                                           >  6:211325/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCCTGTGCGAGATGGCCGGCGGAACTTCCATCATGCACCAGGCGGCGCGGTGCAGCGGGTGGGTGTGGATGATGCAGATGCCGTCGgggctg                                                           >  6:6945/1‑133 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCCTCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcggg                                                          >  8:69873/1‑136 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  5:133600/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  6:166525/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  5:222020/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  2:183463/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  3:141018/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgg                                                           >  4:133373/1‑135 (MQ=255)
                                                                         cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGAAGGCCAGCGGAACTTCCCGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGCTGTGGATGATGCAGTTGCCGTCGgggcgg                                                           >  1:305687/1‑135 (MQ=255)
                                                                           ttCCTTGAGGAATGCACAGGCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggcgggc                                                         <  3:78919/135‑1 (MQ=255)
                                                                            tCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACtt                                                                                                                            >  4:335446/1‑67 (MQ=255)
                                                                            tCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACtt                                                                                                                            <  3:335446/67‑1 (MQ=255)
                                                                              cTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGgcgc                                                       <  4:178397/134‑1 (MQ=255)
                                                                                 gAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGt                                                    >  4:89390/1‑134 (MQ=255)
                                                                                   ggAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgatg                                                                                >  4:252515/1‑104 (MQ=255)
                                                                                   ggAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgatg                                                                                <  3:252515/104‑1 (MQ=255)
                                                                                   ggAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGAGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTAcag                                                >  8:259058/1‑136 (MQ=255)
                                                                                         cacaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTAcagccagc                                           >  2:13281/1‑135 (MQ=255)
                                                                                          acaTTCGTCGAAAACGGGGAACAGGTCCATGTGCCACATGTCCAGCGGAACTTCCAGCATAGACAATGCCGCGCAGTGCACCGGCTGGGTGTCGATGATGCATTTGACGTAGCGGAGGGCGCGTTACAGACACCt                                          >  6:30836/1‑135 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAggcgg                                                                                                          >  1:347114/1‑71 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAggcgg                                                                                                          <  2:347114/71‑1 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCgggtggg                                                                                        <  6:258770/89‑1 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCgggtggg                                                                                        >  5:258770/1‑89 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtg                                         >  8:212535/1‑136 (MQ=255)
                                                                                          acaGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtg                                         >  7:64754/1‑136 (MQ=255)
                                                                                           caGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtgt                                        >  8:30563/1‑136 (MQ=255)
                                                                                            aGTCGTCGACCAGGGGGCACAGGTCCATGTGCGAGCTGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGGGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtgt                                        <  1:148127/135‑1 (MQ=255)
                                                                                            aGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCCGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGCCGTCGGGGCGGGCGCGGTACCGCCAGCtgt                                        >  1:393468/1‑135 (MQ=255)
                                                                                            aGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtgt                                        >  6:117724/1‑135 (MQ=255)
                                                                                                    aaCAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCgg                                >  6:337493/1‑135 (MQ=255)
                                                                                                    aaCAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGt                               >  7:21676/1‑136 (MQ=255)
                                                                                                         ggggCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGc                           <  6:173492/135‑1 (MQ=255)
                                                                                                          gggCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGcc                          <  5:337493/135‑1 (MQ=255)
                                                                                                              acaGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGt                      >  5:162592/1‑135 (MQ=255)
                                                                                                               caGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTg                    <  7:47701/136‑1 (MQ=255)
                                                                                                                aGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTg                    <  1:13281/135‑1 (MQ=255)
                                                                                                                 ggTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTgg                   >  1:376190/1‑135 (MQ=255)
                                                                                                                      ttgtgCGCGATGGCCAGCGCAACGTCCAGCAGCGCCAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGgccatg              <  2:393468/134‑1 (MQ=255)
                                                                                                                      aTGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCGGTGCAGCGGGTGGGTGTGGATGTTGCAGTTGACGTCGGGGCGGGCGGGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGgccat               <  5:211325/134‑1 (MQ=255)
                                                                                                                        gtgCGAGATGGGCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAg                                                                                                >  6:310149/1‑51 (MQ=255)
                                                                                                                        gtgCGAGATGGGCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAg                                                                                                >  6:305101/1‑51 (MQ=255)
                                                                                                                        gtgCGAGATGGGCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAg                                                                                                <  5:305101/51‑1 (MQ=255)
                                                                                                                        gtgCGAGATGGGCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAg                                                                                                <  5:310149/51‑1 (MQ=255)
                                                                                                                          gCGAGATGGCCAGCGGAACTTCCAGAATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCTGGGCGGGCTCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGgccatgccat          >  4:305308/1‑135 (MQ=255)
                                                                                                                              gaTGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGCCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGGACATGCCATGGCCt     >  8:15139/1‑136 (MQ=255)
                                                                                                                                 ggCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGGCCATGCCATGGCCTTgc  <  7:45026/136‑1 (MQ=255)
                                                                                                                                   ccAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGGCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGGCCATGCCATGGCCTTgc  <  3:124830/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
CAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGTGGAAGCGGTTGGCCGGGTTGGCCATGCCATGGCCTTGC  >  ACB122/2308469‑2308733

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