Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,297 A→T 10.5% intergenic (‑236/+19) PP_3989 ← / ← PP_t60 DNA‑cytosine methyltransferase/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,501,2970AT10.5% 155.1 / 8.6 60intergenic (‑236/+19)PP_3989/PP_t60DNA‑cytosine methyltransferase/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base T (4/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

GAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACAC  >  NC_002947/4501150‑4501442
                                                                                                                                                   |                                                                                                                                                 
gAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGaa                                                                                                                                                   >  5:290643/1‑148 (MQ=255)
     gCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCa                                                                                                                                             <  2:226812/149‑1 (MQ=255)
       tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAgg                                                                                                                                           >  3:233866/1‑149 (MQ=255)
       tctAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAgg                                                                                                                                           >  1:20004/1‑149 (MQ=255)
         tAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGc                                                                                                                                         >  7:106799/1‑149 (MQ=255)
            gggAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                      >  8:248246/1‑149 (MQ=255)
            gggAAGCGCGTCACGCCAAGCTGAAGCAGCTCACACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                      >  3:309715/1‑149 (MQ=255)
              gAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGTAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCATGGCttttt                                                                                                                                    <  4:233866/149‑3 (MQ=255)
              gAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGGGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTg                                                                                                                                     >  5:294974/1‑148 (MQ=255)
              gAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGATTATTCAGGGCttttt                                                                                                                                    <  3:268818/149‑3 (MQ=255)
                aGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAat                                                                                                                                  <  8:205313/149‑1 (MQ=255)
                aGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAat                                                                                                                                  <  3:333706/149‑1 (MQ=255)
                  cgcgTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAatat                                                                                                                                >  2:67723/1‑149 (MQ=255)
                     gTCACGCCAAGCTGAAGCAGCTCACACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGc                                                                                                                             <  4:309715/149‑1 (MQ=255)
                        aCGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATAtggtgg                                                                                                                           >  7:153599/1‑148 (MQ=255)
                        aCGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCgg                                                                                                                           >  7:111303/1‑148 (MQ=255)
                         cGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                         <  2:185531/149‑1 (MQ=255)
                                  gAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                    >  3:118549/1‑145 (MQ=255)
                                     gcagcTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGGATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAGGCGTAGAGATTc                                                                                                             >  2:133865/1‑149 (MQ=255)
                                     gcagcTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTc                                                                                                             <  6:290643/149‑1 (MQ=255)
                                        gcTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGaa                                                                                                          <  3:182705/149‑1 (MQ=255)
                                          tCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAAc                                                                                                         >  2:21227/1‑148 (MQ=255)
                                            gCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAActc                                                                                                       <  3:189254/148‑1 (MQ=255)
                                                  aGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATa                                                                                                <  8:111303/149‑1 (MQ=255)
                                                     aTGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGtt                                                                                             >  8:172232/1‑149 (MQ=255)
                                                        tgtCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCATGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGc                                                                                           <  8:239869/148‑1 (MQ=255)
                                                                 cGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACg                                                                                  <  7:248246/148‑1 (MQ=255)
                                                                   gAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                         >  6:241229/1‑107 (MQ=255)
                                                                   gAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                         <  5:241229/107‑1 (MQ=255)
                                                                    aaaCGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGCCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGttt                                                                              >  7:247496/1‑149 (MQ=255)
                                                                      aCGATCAGAGAATTGCGGAAACGACCCGTATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTc                                                                            <  7:271745/149‑1 (MQ=255)
                                                                      aCGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGTGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTc                                                                            <  8:153599/149‑1 (MQ=255)
                                                                        gATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCaa                                                                          <  1:21227/149‑1 (MQ=255)
                                                                             gagaATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATc                                                                                                                                              >  5:279680/1‑76 (MQ=255)
                                                                             gagaATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATc                                                                                                                                              <  6:279680/76‑1 (MQ=255)
                                                                             gagaATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCg                                                                     >  7:10731/1‑149 (MQ=255)
                                                                                 aTTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGc                                                                 <  7:172232/149‑1 (MQ=255)
                                                                                    gCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCt                                                               <  2:20004/148‑1 (MQ=255)
                                                                                     cGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTa                                                             <  1:67723/149‑1 (MQ=255)
                                                                                       gAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTaaa                                                           <  6:294974/149‑1 (MQ=255)
                                                                                            gACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAGGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACt                                                      <  1:133865/149‑1 (MQ=255)
                                                                                              cccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCt                                                               <  3:249065/138‑1 (MQ=255)
                                                                                              cccGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCt                                                               >  4:249065/1‑138 (MQ=255)
                                                                                                        aGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTc                                          >  7:138020/1‑149 (MQ=255)
                                                                                                              gAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTc                                    >  4:79925/1‑149 (MQ=255)
                                                                                                               aaaaaaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGATTATGGCGGAAGCGTAGAGATTCGAACTCTAGGATTGTTGCCCATCGACGGTTTTCAAGACCGTTGCGTTAATCCACTCGGCCACGCTTCCATCTCg                                   >  2:112060/1‑149 (MQ=37)
                                                                                                               aaaaaaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCg                                   >  8:165167/1‑149 (MQ=255)
                                                                                                                   aaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGtttt                               <  4:118549/149‑1 (MQ=255)
                                                                                                                   aaCTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAAAGATTAGAACTCTAGGATATTTTGCCATAGACGGTTTTGAAGACCGTTTGCTTTAAACACTCGGCCACGCTTCCAGCTCGtttt                               >  1:275329/1‑149 (MQ=25)
                                                                                                                            cTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCGATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgc                       >  6:292251/1‑148 (MQ=255)
                                                                                                                               aGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGGCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGgccgccta                    >  7:61338/1‑146 (MQ=255)
                                                                                                                                    tattcaaaGCCCTGAGTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGGCCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATAcc              >  1:128987/6‑149 (MQ=255)
                                                                                                                                       aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGACCTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGGCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGTCCGCCATAATACCtt            >  5:240855/1‑148 (MQ=255)
                                                                                                                                       aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCtt            >  7:109941/1‑148 (MQ=255)
                                                                                                                                          aaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATg        <  8:10731/149‑1 (MQ=255)
                                                                                                                                           aGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATg        >  3:111261/1‑148 (MQ=255)
                                                                                                                                           aGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCAATCGGCCACCCTTCCAGCTCGTTTTGCGCCCGCAATAATACCTTAATg        >  6:25494/1‑148 (MQ=38)
                                                                                                                                            gCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAg                                       >  6:110584/1‑116 (MQ=255)
                                                                                                                                               cTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAaca   >  7:177822/1‑149 (MQ=255)
                                                                                                                                               cTGAATAATCAGGGCTGTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAaca   >  6:226460/1‑149 (MQ=255)
                                                                                                                                                tGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAacac  >  6:313003/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAAAGCTCTAGGGGAAGCGCGTCACGCCAAGCTGAAGCAGCTCGCACACAGCATGTGTCGGATGCGGAAACGATCAGAGAATTGCGGAAACGACCCGGATGACAGAAGGGAAAAAACTAGAGCCTCAGAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCGAACTCTAGGATAGTTGCCCATCGACGGTTTTCAAGACCGTTGCCTTAAACCACTCGGCCACGCTTCCAGCTCGTTTTGCGGCCGCCATAATACCTTAATGAAACAC  >  NC_002947/4501150‑4501442

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: