Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,960,052 A→T 11.3% noncoding (73/103 nt) PP_mr63 ← rnk_pseudo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,960,0520AT11.3% 198.1 / 16.2 71noncoding (73/103 nt)PP_mr63rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base T (4/4);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

CGTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959905‑5960197
                                                                                                                                                   |                                                                                                                                                 
cGTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAgg                                                                                                                                                 <  1:330256/150‑1 (MQ=255)
     tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            >  1:566471/1‑150 (MQ=255)
      ctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCgggg                                                                                                                                           <  2:816308/150‑1 (MQ=255)
       tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGt                                                                                                                                          >  1:334636/1‑150 (MQ=255)
        cgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtt                                                                                                                                         <  2:480814/150‑1 (MQ=255)
         gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCgg                                                                                                                                             >  1:801014/1‑145 (MQ=255)
         gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCgg                                                                                                                                             <  2:801014/145‑1 (MQ=255)
           cgAGGCTGTCGATCATACGCTCCAGACGCTGTACGTCCAATCGGGTCACGATTACGAATTGCTCGCTCCTCGCGATCCAGACAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                      <  1:678843/150‑1 (MQ=255)
            gAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTgag                                                                                                                                     <  1:301167/150‑1 (MQ=255)
              ggCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGaa                                                                                                                                   <  1:869309/150‑1 (MQ=255)
                     tatcggccGCTCCAGACGCTGGACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                            <  2:1105934/143‑1 (MQ=255)
                       tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                          <  2:1425180/150‑1 (MQ=255)
                       tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGACTAGTTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTgga                                                                                                                          >  1:1105934/1‑149 (MQ=255)
                                  aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGaa                                                                                                                                   >  1:1238517/1‑130 (MQ=255)
                                  aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGaa                                                                                                                                   <  2:1238517/130‑1 (MQ=255)
                                    aCGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAg                                                                                                             >  1:527353/1‑150 (MQ=255)
                                       cTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAAc                                                                                                          <  1:273937/150‑1 (MQ=255)
                                        tGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAGCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAAcc                                                                                                         >  2:1088302/1‑150 (MQ=255)
                                         gTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCTGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAcccc                                                                                                        <  1:36777/150‑1 (MQ=255)
                                          tACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                      >  1:754507/1‑119 (MQ=255)
                                          tACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCg                                                                                                       <  2:334636/150‑1 (MQ=255)
                                          tACGTCCAATCGGGTGAGGATGAGGGAAGGCTTCGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                      <  2:754507/119‑1 (MQ=255)
                                             gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGtt                                                                                                                                         <  2:1065843/113‑1 (MQ=255)
                                             gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGtt                                                                                                                                         >  1:1065843/1‑113 (MQ=255)
                                             gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGaca                                                                                                    >  2:228256/1‑150 (MQ=255)
                                                cAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAc                                                                                                <  2:566471/150‑1 (MQ=255)
                                                      ggTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAg                                                                                           <  2:1362192/150‑1 (MQ=255)
                                                            gATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                     >  2:1377281/1‑150 (MQ=255)
                                                            gATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                     >  1:1306379/1‑150 (MQ=255)
                                                                  ggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATAc                                                                               >  2:782458/1‑150 (MQ=255)
                                                                    aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGt                                                                                                                                          <  1:494586/89‑1 (MQ=255)
                                                                    aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGt                                                                                                                                          >  2:494586/1‑89 (MQ=255)
                                                                    aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTc                                                                             <  2:527353/150‑1 (MQ=255)
                                                                       gCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                      <  1:44241/90‑1 (MQ=255)
                                                                       gCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                      >  2:44241/1‑90 (MQ=255)
                                                                                tCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            >  1:1338736/1‑75 (MQ=255)
                                                                                tCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            <  2:1338736/75‑1 (MQ=255)
                                                                                    gATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAg                                                             >  2:1208555/1‑150 (MQ=255)
                                                                                         aGTCAGCCTCCTTTGAACAGGCAGGTTCCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGc                                                        <  1:191548/150‑1 (MQ=255)
                                                                                         aGTCAGACTCCTTTGAACAGGCAGTTTTCCATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGc                                                        >  1:545263/1‑150 (MQ=255)
                                                                                               aCTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCACAAGTGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGc                                                  >  1:1034549/1‑150 (MQ=255)
                                                                                                 tCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCt                                                >  1:643820/1‑150 (MQ=255)
                                                                                                 tCCTTTGAACAGGCAGTTGTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGATTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAACACATCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCt                                                >  2:293305/1‑150 (MQ=255)
                                                                                                   cTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTgg                                                                                                                           <  1:7140/73‑1 (MQ=255)
                                                                                                   cTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTgg                                                                                                                           >  2:7140/1‑73 (MQ=255)
                                                                                                    tttGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGTGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCt                                             >  1:871367/1‑150 (MQ=255)
                                                                                                    tttGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                      >  1:72303/1‑93 (MQ=255)
                                                                                                    tttGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                      <  2:72303/93‑1 (MQ=255)
                                                                                                     ttGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTc                                            >  2:239640/1‑150 (MQ=255)
                                                                                                       gaccAGGCAGCTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTCCGCTCACTGTGCCTTAGCCTGCTcac                                          >  2:533099/4‑148 (MQ=255)
                                                                                                        aaCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTcgcg                                         >  1:303540/1‑150 (MQ=255)
                                                                                                          cAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCgg                                                                                       >  1:286864/1‑102 (MQ=255)
                                                                                                          cAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCgg                                                                                       <  2:286864/102‑1 (MQ=255)
                                                                                                           aGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCaa                                      <  1:782458/150‑1 (MQ=255)
                                                                                                                  tttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTgg                               <  1:1377281/150‑1 (MQ=255)
                                                                                                                     cAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTg                                >  2:851723/1‑146 (MQ=255)
                                                                                                                     cAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTg                                <  1:851723/146‑1 (MQ=255)
                                                                                                                      aaTGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGc                                                                         >  2:625042/1‑104 (MQ=255)
                                                                                                                      aaTGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGc                                                                         <  1:625042/104‑1 (MQ=255)
                                                                                                                         gtgCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgc                        <  1:239640/150‑1 (MQ=255)
                                                                                                                             aaaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                    >  1:703180/1‑150 (MQ=255)
                                                                                                                                   aGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTg              >  1:608387/1‑150 (MQ=255)
                                                                                                                                      aaaaCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCGCAGTGTGCCTTACCCTGCTCTCGCACGACCTGGGCCTCGCGGCAGATCTGccc           >  2:648001/1‑149 (MQ=255)
                                                                                                                                        agccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgc         <  1:1088302/148‑1 (MQ=255)
                                                                                                                                         aCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcg        >  1:641816/1‑150 (MQ=255)
                                                                                                                                         aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcg        >  2:1369621/1‑150 (MQ=255)
                                                                                                                                           cccGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCGCCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc      >  1:1151729/1‑150 (MQ=255)
                                                                                                                                           cccGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                           <  1:1186718/81‑1 (MQ=255)
                                                                                                                                           cccGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                           >  2:1186718/1‑81 (MQ=255)
                                                                                                                                           cccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTAAACGCCGAGGCTCGACATACTACCACAGCCGGGCAAAAACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc      >  1:1096918/1‑150 (MQ=255)
                                                                                                                                               cccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTc  <  1:1208555/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959905‑5960197

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