Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 396,862:1 +G 6.1% coding (84/924 nt) PP_0329 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947396,8621.G 6.1% 265.9 / 11.1 66coding (84/924 nt)PP_0329hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (33/29);  new base G (2/2);  total (35/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 = 9.69e-01

AAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGATTATCT  >  NC_002947/396715‑397008
                                                                                                                                                    |                                                                                                                                                  
aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                                     >  2:215475/1‑148 (MQ=255)
  aGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGaa                                                                                                                                                 >  6:283884/1‑149 (MQ=255)
             gtgtGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATc                                                                                                                                      <  4:115861/149‑1 (MQ=255)
                 gATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGtt                                                                                                                                   >  4:199858/1‑148 (MQ=255)
                  aTTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGttt                                                                                                                                  <  4:185646/148‑1 (MQ=255)
                               ccGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGCAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTgg                                                                                                                     >  3:215516/1‑148 (MQ=255)
                                cGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGt                                                                                                                    <  1:163652/148‑1 (MQ=255)
                                 gTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGt                                                                                                                    >  2:221384/1‑147 (MQ=255)
                                       cGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                            <  5:283884/149‑1 (MQ=255)
                                         tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                                          <  1:215475/149‑1 (MQ=255)
                                           cTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCGGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggt                                                                                                        >  1:245638/1‑149 (MQ=255)
                                              aCGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                                  >  2:2504/1‑136 (MQ=255)
                                                gACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGtt                                                                                                                                   >  2:290517/1‑117 (MQ=255)
                                                   aaGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                                 >  8:2230/1‑148 (MQ=255)
                                                     gCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaaa                                                                                              >  4:287943/1‑149 (MQ=255)
                                                            gTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAATAGTGAGCGTTCCAATCCGGAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                       <  4:311496/149‑1 (MQ=255)
                                                               catctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                                     <  5:149343/148‑1 (MQ=255)
                                                                  ctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGt                                                                                                                <  3:222419/118‑1 (MQ=255)
                                                                  ctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGt                                                                                                                >  4:222419/1‑118 (MQ=255)
                                                                       ctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                                  >  2:201283/1‑111 (MQ=255)
                                                                       ctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                                  <  1:201283/111‑1 (MQ=255)
                                                                          ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                                 <  6:215126/125‑1 (MQ=255)
                                                                          ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                                 >  5:215126/1‑125 (MQ=255)
                                                                          ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGa                                                                            <  3:199858/146‑1 (MQ=255)
                                                                          ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTg                                                                         <  7:2230/149‑1 (MQ=255)
                                                                          ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTg                                                                         >  8:251703/1‑149 (MQ=255)
                                                                                aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                                              >  2:174440/1‑138 (MQ=255)
                                                                                aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                                              <  1:174440/138‑1 (MQ=255)
                                                                                 aaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                                  <  3:287943/149‑1 (MQ=255)
                                                                                     tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                                          <  7:97986/105‑1 (MQ=255)
                                                                                     tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                                          >  8:97986/1‑105 (MQ=255)
                                                                                       ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                             <  1:221384/148‑1 (MQ=255)
                                                                                       ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                             <  7:251703/148‑1 (MQ=255)
                                                                                         aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTg                                                          >  3:268545/1‑149 (MQ=255)
                                                                                            ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                            >  8:204146/1‑96 (MQ=255)
                                                                                            ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                            <  7:204146/96‑1 (MQ=255)
                                                                                              ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaa                                                                                               >  1:293238/1‑107 (MQ=255)
                                                                                              ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaa                                                                                               <  2:293238/107‑1 (MQ=255)
                                                                                              ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCAGCAGGATTGACCGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCCCCGACGTGCTccc                                                     >  6:261886/1‑149 (MQ=255)
                                                                                                tGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGt                                                   >  5:187968/1‑149 (MQ=255)
                                                                                                tGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGGCAGTTTGGTGGTCTGGTAGAAATATTTGGGGGTATCAGGAATGTTTTTCACCGCCGTGGTcaggt                                                   >  4:20511/1‑145 (MQ=255)
                                                                                                  cttcttTGCCATCGGACTCAGCAAGAGGGATCTTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTTGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                                  <  6:187968/148‑1 (MQ=255)
                                                                                                  cttcttTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCATTCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                                  <  6:81518/148‑1 (MQ=255)
                                                                                                        tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                                     <  6:108503/91‑1 (MQ=255)
                                                                                                        tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                                     >  5:108503/1‑91 (MQ=255)
                                                                                                          ccAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTggg                                                                                      >  2:309490/1‑104 (MQ=255)
                                                                                                          ccAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTggg                                                                                      <  1:309490/104‑1 (MQ=255)
                                                                                                           cAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                                        <  4:215516/149‑1 (MQ=255)
                                                                                                             gCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                                      >  3:77929/1‑149 (MQ=255)
                                                                                                                      gCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGGTCCCGTTGTCTGCAGTTGGCCCCGGaaa                             >  1:280136/1‑149 (MQ=255)
                                                                                                                      gCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                             <  5:68077/149‑1 (MQ=255)
                                                                                                                       cAAGATGGAGCGTTCAAATCCGCATGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTg                                             <  7:125434/132‑1 (MQ=255)
                                                                                                                       cAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTATGGTATAAATATTTGTGGGTATCAGTACTGTTTTTCACCGACGTTCTCCCGTTGTCTg                                             >  8:125434/1‑132 (MQ=255)
                                                                                                                          gagGGATCGTTCCAATCCGCCGGATTGAACGCTTCTGATCGTTTCCAATGCCAGTGGTCAGTTTGCTGGTCTGGTTGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTTTATGCAGTTGGCCCCGGAAATAGa                         <  7:48810/149‑1 (MQ=255)
                                                                                                                           cagaaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                                  <  3:297252/103‑1 (MQ=37)
                                                                                                                           cagaaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                                  >  4:297252/5‑107 (MQ=37)
                                                                                                                             ggAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg                      >  5:16322/1‑149 (MQ=255)
                                                                                                                             ggAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg                      >  6:196248/1‑149 (MQ=255)
                                                                                                                              aaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                                      <  4:280657/131‑1 (MQ=255)
                                                                                                                              aaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                                      >  3:280657/2‑132 (MQ=255)
                                                                                                                                  gTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCAGCAGTTGGCCCCGCAAATAGATGGCCTgg                 >  7:289732/1‑149 (MQ=255)
                                                                                                                                     ccAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACAGACGTGCTCCAGTTGTCTGCAGTTGGCCCCTGAAATAGATGGCATGGGc               >  2:216847/1‑148 (MQ=255)
                                                                                                                                        aTCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGTGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTa           <  2:245638/149‑1 (MQ=255)
                                                                                                                                             cAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACt            <  6:271188/143‑1 (MQ=255)
                                                                                                                                             cAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACt            >  5:271188/1‑143 (MQ=255)
                                                                                                                                               ggATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGATTa     >  6:28606/1‑148 (MQ=255)
                                                                                                                                                 aTTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGATTATCt  >  7:95410/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGATTATCT  >  NC_002947/396715‑397008

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