Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,999,564 C→A 9.6% A33A (GCC→GCA PP_5256 → cyclic nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,999,5640CA 9.6% 203.5 / 10.5 73A33A (GCC→GCAPP_5256cyclic nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/28);  new base A (4/3);  total (42/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 = 4.94e-01

AAACCTGCTCTATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCGGCTGGTGGCTGT  >  NC_002947/5999416‑5999712
                                                                                                                                                     |                                                                                                                                                    
aaaCCTGCTCTATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGcc                                                                                                                                                      >  5:83469/1‑149 (MQ=255)
      gCTCTATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGgcg                                                                                                                                                 >  5:148323/1‑148 (MQ=255)
        tctATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCAt                                                                                                                                              >  7:110630/1‑149 (MQ=255)
        tctATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGATCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGACATGCCTGGCGCa                                                                                                                                               >  2:84892/1‑148 (MQ=255)
         ctATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCAt                                                                                                                                              <  8:249406/148‑1 (MQ=255)
         ctATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCAt                                                                                                                                              >  7:249406/1‑148 (MQ=255)
          tATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCAt                                                                                                                                              <  3:131401/147‑1 (MQ=255)
           aTTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCACCGCCATGCCTGGCGCATc                                                                                                                                             <  7:182703/147‑1 (MQ=255)
           aTTACTACCCACAGCCCTCCCCTCACCCCCGAGATCAACATGCTGGCGTTCATCCAGTCGTACCAGTTGTTGGTCGGTTGCGCCTTGATCCTGCTAGACCTCGTGCTTTGGCAACTGAGACCCATTAAGCGCCATGCCTGGCGCAACaa                                                                                                                                           >  8:182703/1‑148 (MQ=255)
            ttACTACCCTCAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGAGCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCAGGCATGgcgctgaatg                                                                                                                                          <  8:195736/149‑7 (MQ=255)
            ttACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCa                                                                                                                                            <  8:110630/147‑1 (MQ=255)
                  cccACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAgcgcgcgc                                                                                                                                    >  6:274214/1‑149 (MQ=255)
                    cacaGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCt                                                                                                                                  >  8:231016/1‑149 (MQ=255)
                         ccctcccctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTg                                                                                                                              <  2:90302/148‑1 (MQ=255)
                         ccctcccctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTg                                                                                                                              >  1:90302/1‑148 (MQ=255)
                         ccctcccctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGa                                                                                                                             >  8:281012/1‑149 (MQ=255)
                           ctcccctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATc                                                                                                                           >  1:96149/1‑149 (MQ=255)
                              ccctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCt                                                                                                                          <  4:207877/147‑1 (MQ=255)
                               cctcACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTcc                                                                                                                       <  3:264662/149‑1 (MQ=255)
                                      cccGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGtt                                                                                                                >  1:44460/1‑149 (MQ=255)
                                            aTCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCt                                                                                                           >  1:53884/1‑148 (MQ=255)
                                               aCCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAg                                                                                                        <  3:65873/148‑1 (MQ=255)
                                               aCCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGt                                                                                                       <  7:231016/149‑1 (MQ=255)
                                                ccATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTg                                                                                                      >  4:40075/1‑149 (MQ=255)
                                                     cTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGa                                                                                                  <  2:53884/148‑1 (MQ=255)
                                                        gCGTTCATCCATTCGTCCCCGTTTTTGCTCGGTTGGGCCTTGATCCTGCTCGACCTCTTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCAGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATgg                                                                                               <  1:84892/148‑1 (MQ=255)
                                                          gTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCg                                                                                            >  1:302897/1‑149 (MQ=255)
                                                           ttCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGc                                                                                           >  1:295774/1‑149 (MQ=255)
                                                           ttCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGc                                                                                           >  1:291918/1‑149 (MQ=255)
                                                              aTCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTgg                                                                                        >  6:319442/1‑149 (MQ=255)
                                                                ccAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTGTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTgg                                                                                        >  1:89522/1‑148 (MQ=255)
                                                                ccAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCa                                                                                      <  3:143039/149‑1 (MQ=255)
                                                                 cAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCa                                                                                      <  3:40075/148‑1 (MQ=255)
                                                                   gtcgtcCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCAt                                                                                     >  1:47644/1‑147 (MQ=255)
                                                                   gtcgtcCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGa                                                                                   <  7:281012/149‑1 (MQ=255)
                                                                         cccGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGc                                                                             >  3:242365/1‑149 (MQ=255)
                                                                          ccGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGTTGATGGCCGCTGGCATGAGCCCGCt                                                                            <  7:27518/149‑1 (MQ=255)
                                                                          ccGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCt                                                                            >  2:162344/1‑149 (MQ=255)
                                                                             ttgttgCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCa                                                                         <  5:319442/149‑1 (MQ=255)
                                                                               gttgCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCaa                                                                        <  5:274214/148‑1 (MQ=255)
                                                                                    tCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTg                                                                                                      <  3:7790/113‑1 (MQ=255)
                                                                                    tCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTg                                                                                                      >  4:7790/1‑113 (MQ=255)
                                                                                          gcgcCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTg                                                             <  2:275806/148‑1 (MQ=255)
                                                                                           cgcCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGc                                                           >  4:115949/1‑149 (MQ=255)
                                                                                            gcCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGcc                                                          >  8:111298/1‑149 (MQ=255)
                                                                                              cTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGccc                                                         >  3:295997/1‑148 (MQ=255)
                                                                                              cTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGAGCGCCTGGTGATATTCCTGCAGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGccc                                                         >  1:99984/1‑148 (MQ=255)
                                                                                                 gATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCcgac                                                      <  1:162344/148‑1 (MQ=255)
                                                                                                      tGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGtc                                                >  8:26746/1‑149 (MQ=255)
                                                                                                      tGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGtc                                                >  4:151814/1‑149 (MQ=255)
                                                                                                         tCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTcg                                             >  5:190320/1‑149 (MQ=255)
                                                                                                           gACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTcgcg                                           <  4:242365/149‑1 (MQ=255)
                                                                                                           gACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTcg                                             <  4:295997/147‑1 (MQ=255)
                                                                                                           gACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTcg                                             <  2:47644/147‑1 (MQ=255)
                                                                                                                 gTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATCAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAAc                                       <  1:228079/147‑1 (MQ=255)
                                                                                                                     tttGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATg                                 >  4:280993/1‑149 (MQ=255)
                                                                                                                        tGGCAACTGATCCCCATTCAGCGCCATGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGc                               >  7:47979/1‑149 (MQ=255)
                                                                                                                           cAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCAc                            >  1:147003/1‑149 (MQ=255)
                                                                                                                           cAAATGATCCCCATTCAGCGGCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTAAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCTTCCTGATGGCCAc                            <  2:99984/149‑1 (MQ=255)
                                                                                                                            aaCTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACg                           <  6:148323/149‑1 (MQ=255)
                                                                                                                                 aTCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTa                        <  2:89522/147‑1 (MQ=255)
                                                                                                                                 aTCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTAtt                      >  2:97744/1‑149 (MQ=255)
                                                                                                                                 aTCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTAtt                      >  2:47892/1‑149 (MQ=255)
                                                                                                                                  tCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTAtt                      <  2:147003/148‑1 (MQ=255)
                                                                                                                                  tCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGCTCATCTTCCTCCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTg                     >  6:71822/1‑149 (MQ=255)
                                                                                                                                       aTTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCAt                >  6:164861/1‑149 (MQ=255)
                                                                                                                                       aTTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCAt                >  5:65968/1‑149 (MQ=255)
                                                                                                                                        tgatgcgcCAGGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATc               >  4:86208/5‑149 (MQ=255)
                                                                                                                                        tgatgcgcCAGGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATc               >  2:243757/5‑149 (MQ=255)
                                                                                                                                          atgcgcCAGGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCgg             >  6:135255/3‑149 (MQ=255)
                                                                                                                                           aGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCGGc            >  5:292645/1‑149 (MQ=255)
                                                                                                                                             cgcCAGGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCGGCtg          <  3:86208/149‑1 (MQ=255)
                                                                                                                                              gcCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGCTGGCCACGGTATTGGCCATCGGCtgg         >  2:96990/1‑149 (MQ=255)
                                                                                                                                                 aTGCATGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCTATCGGCtggtgg      <  8:47979/149‑1 (MQ=255)
                                                                                                                                                     cTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCGGCTGGTGGCTGt  >  1:237228/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACCTGCTCTATTACTACCCACAGCCCTCCCCTCACCCCCGAGATCACCATGCTGGCGTTCATCCAGTCGTCCCCGTTGTTGCTCGGTTGCGCCTTGATCCTGCTCGACCTCGTGCTTTGGCAACTGATCCCCATTCAGCGCCATGCCTGGCGCATCAGCGCGCGCCTGGTGATCTTCCTGCTGTTCAGCTCAGTGCTGATGGCCGCTGGCATGAGCCCGCTGCAACCACCGCCGTGGCCCGACGATGTCTCGCGCAACCTGATGGCCACGGTATTGGCCATCGGCTGGTGGCTGT  >  NC_002947/5999416‑5999712

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: