Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 640,263 A→T 8.2% L161Q (CTG→CAG)  PP_0551 ← MORN domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947640,2630AT 8.2% 169.4 / 5.6 61L161Q (CTG→CAG) PP_0551MORN domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (2/3);  total (30/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 = 8.62e-01

CGTCGGCCAGTTCCAGTTGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCC  >  NC_002947/640117‑640409
                                                                                                                                                  |                                                                                                                                                    
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)
                                                                                                                                                  |                                                                                                                                                    
CGTCGGCCAGTTCCAGTTGGCCTGCACCGTCGTACAGGTCGTCCTTGAACTGGCCGCGGTATGTCTGGTCGGCTTGCTTGAGGGTGCCCTCGCCGTCGCGGCGGCCGTGGCGGAAGGTGCCGCTATAGTGGCTGCCTGGGGTGATCAGGTCA‑CCCAGGCCCTGAAACAGCCCCTCGGCGAACTGGCCGCGGTAGACCTCGCCGTTCTGGCCATGCCATTCGCCCTGGCC‑GTGCCACTGGCCGTCCTTGAAGCTGCCGGCGTACCAGCTGCCGTTGGGGTAGTCGATCCGGCCC  >  NC_002947/640117‑640409

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