Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 640,265 G→A 8.3% D160D (GAC→GAT)  PP_0551 ← MORN domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947640,2650GA 8.3% 169.1 / 3.1 60D160D (GAC→GAT) PP_0551MORN domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/3);  total (30/30)
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.93e-01

CGTCGGCCAGTTCCAGTTGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCCTCA  >  NC_002947/640117‑640412
                                                                                                                                                    |                                                                                                                                                     
cGTCGGCCAGTTCCAGTTGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAgg                                                                                                                                                       <  1:94402/149‑1 (MQ=255)
       cagttccagttAGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCa                                                                                                                                               <  8:88606/149‑1 (MQ=255)
        agttccagttGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCT‑GGGTGACCTGATCA‑ccccag                                                                                                                                             <  7:14066/149‑4 (MQ=255)
          ttccagttGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCT‑GGGTGACCTGATCA‑ccccaggc                                                                                                                                           <  1:110163/149‑6 (MQ=255)
          ttccagttGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCT‑GGGTGACCTGATCA‑ccccaggc                                                                                                                                           <  4:47567/149‑6 (MQ=255)
            ccagttGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGcc                                                                                                                                           <  7:56767/148‑1 (MQ=255)
             cagttGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCt                                                                                                                                         <  3:317016/149‑1 (MQ=255)
                   gCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAAc                                                                                                                                    <  4:139644/148‑1 (MQ=255)
                           cgtcgtACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCg                                                                                                                           >  5:39910/1‑149 (MQ=255)
                            gtcgtACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCgg                                                                                                                          <  8:90437/149‑1 (MQ=255)
                                 aCAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGAGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAAc                                                                                                                     >  4:32211/1‑149 (MQ=255)
                                 aCAGGTCGTCCTTGAACTGGCCGCGGTATGACTCGTCGGAATGCTTGAGCGTGACCACCCCGTCGCGGCGCACGTCGCGGAAGGTGCAGCAATAGTAGCAGCCTGGGGTGATCAGGTCC‑CCCAGCCCAGGAAACAGCCACTCGGCGaa                                                                                                                      >  2:240362/1‑148 (MQ=255)
                                   aGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTg                                                                                                                   >  6:269289/1‑149 (MQ=255)
                                       cgtcCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCAGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCcg                                                                                                               >  1:142032/1‑149 (MQ=255)
                                        gtcCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGACCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAGCAGCCCCTCGGCGAACTGGCcgc                                                                                                              >  2:215342/1‑149 (MQ=255)
                                           cTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCgg                                                                                                            <  1:288518/148‑1 (MQ=255)
                                           cTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCgg                                                                                                            <  2:49552/148‑1 (MQ=255)
                                            ttGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCgg                                                                                                            >  1:318046/1‑147 (MQ=255)
                                            ttGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTAGGCGAACTGGCCGCgg                                                                                                            >  5:238287/1‑147 (MQ=255)
                                            ttGAACTGGCCGCGGTATGTCTGGACGGATTGCTTGAGGGTGACCTCGCCGTCGCGGCGGCAGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATAAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCgg                                                                                                            >  2:209741/1‑147 (MQ=255)
                                              gAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTa                                                                                                          <  4:186687/147‑1 (MQ=255)
                                                      cgcgGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCCGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGCCCTCGCCg                                                                                                >  2:212883/1‑149 (MQ=255)
                                                      cgcgGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCg                                                                                                >  7:214445/1‑149 (MQ=255)
                                                             tGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGgc                                                                                         <  4:31839/149‑1 (MQ=255)
                                                                   gTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCTGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGgccat                                                                                      <  2:142032/146‑1 (MQ=255)
                                                                        cttgcttgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAgg                                                                                                                                             <  6:139455/86‑1 (MQ=255)
                                                                        cttgcttgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAgg                                                                                                                                             >  5:139455/1‑86 (MQ=255)
                                                                        cttgcttgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCg                                                                              <  6:120663/149‑1 (MQ=255)
                                                                          tgcttgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGcc                                                                            >  7:268064/1‑149 (MQ=255)
                                                                             ttgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTTGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTg                                                                         >  8:156419/1‑149 (MQ=255)
                                                                              tgAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTgg                                                                        >  4:311545/1‑149 (MQ=255)
                                                                               gAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGc                                                                       >  5:43241/1‑149 (MQ=255)
                                                                                  ggTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑Gt                                                                   <  1:121628/149‑1 (MQ=255)
                                                                                  ggTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑Gt                                                                   >  6:281668/1‑149 (MQ=255)
                                                                                    tGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGc                                                                 <  5:281668/149‑1 (MQ=255)
                                                                                       ccTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCAc                                                              <  1:209741/149‑1 (MQ=255)
                                                                                          cGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTgg                                                           <  1:215342/149‑1 (MQ=255)
                                                                                            ccGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCATGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGcc                                                         <  6:238287/149‑1 (MQ=255)
                                                                                               tCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTc                                                      >  2:55272/1‑149 (MQ=255)
                                                                                                   ggcggcCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTg                                                  <  8:214445/149‑1 (MQ=255)
                                                                                                         cGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTgg                                                           <  2:290157/134‑1 (MQ=255)
                                                                                                         cGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTgg                                                           >  1:290157/1‑134 (MQ=255)
                                                                                                               ggAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCg                                      >  2:88595/1‑149 (MQ=255)
                                                                                                                             taGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGtt                        <  1:88595/149‑1 (MQ=255)
                                                                                                                               gTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTg                       <  3:32211/148‑1 (MQ=255)
                                                                                                                                tGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGTGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTggg                     >  5:103597/1‑149 (MQ=255)
                                                                                                                                tGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTggg                     <  1:55272/149‑1 (MQ=255)
                                                                                                                                tGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTggg                     >  7:86063/1‑149 (MQ=255)
                                                                                                                                  tctggctggGGTGCTCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAAATGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGGCCTGGCCGGTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTggg                     <  1:240362/147‑1 (MQ=255)
                                                                                                                                   cTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTa                  <  6:39910/149‑1 (MQ=255)
                                                                                                                                    tGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGGTGGGGTAg                 >  3:25313/1‑149 (MQ=255)
                                                                                                                                    gggcctgggTGACCTGATCACCCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTgggg                    >  2:217077/7‑147 (MQ=255)
                                                                                                                                     ggcctgggTGACCTGATCACCCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCAGGCGTACCAGCTGCCGTTGGGGTAg                 >  4:120938/6‑149 (MQ=255)
                                                                                                                                      ccTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTTCCGGCGTACCAGCTTCCCTTGGGGTAGTc               >  2:232629/1‑149 (MQ=255)
                                                                                                                                      ccTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGt                <  6:143527/148‑1 (MQ=255)
                                                                                                                                       cTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCg              >  8:64609/1‑149 (MQ=255)
                                                                                                                                          gggTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATc           >  6:214475/1‑149 (MQ=255)
                                                                                                                                               aTCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGTGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGcc      <  6:103597/149‑1 (MQ=255)
                                                                                                                                               aTCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGcc      <  8:268064/149‑1 (MQ=255)
                                                                                                                                               aTCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCCTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGcc      >  1:135329/1‑149 (MQ=255)
                                                                                                                                                tCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGccc     >  7:90040/1‑149 (MQ=255)
                                                                                                                                                   ggTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCCTCa  >  7:293442/1‑149 (MQ=255)
                                                                                                                                                   ggTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCCTCa  >  8:63340/1‑149 (MQ=255)
                                                                                                                                                   ggTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCATGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCCTCa  >  4:249752/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTCGGCCAGTTCCAGTTGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCCTCA  >  NC_002947/640117‑640412

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