Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,364,708 C→T 5.5% A159A (GCC→GCT)  puuB → gamma‑glutamylputrescine oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,364,7080CT 5.5% 229.7 / 2.7 73A159A (GCC→GCT) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACG  >  NC_000913/1364563‑1364855
                                                                                                                                                 |                                                                                                                                                   
tCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAtt                                                                                                                                                  >  3:27398/1‑149 (MQ=255)
  aaaCGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTTGACGCCAACGCCATTc                                                                                                                                                 <  5:347516/148‑1 (MQ=255)
   aaCGTTATCAGATTGACTGCGACTA‑CGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGcc                                                                                                                                              <  5:284064/149‑1 (MQ=255)
         aTCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTg                                                                                                                                           >  7:84262/1‑147 (MQ=255)
            aGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTa                                                                                                                                      <  2:124607/149‑1 (MQ=255)
             gATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAg                                                                                                                                     >  4:309559/1‑149 (MQ=255)
             gATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAAAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAg                                                                                                                                     >  2:310722/1‑149 (MQ=255)
              aTTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGc                                                                                                                                    >  4:290002/1‑149 (MQ=255)
               ttGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                                                                                                                                   >  8:348269/1‑149 (MQ=255)
                    tGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGa                                                                                                                              >  6:187845/1‑149 (MQ=255)
                      cGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACCGCTCGCCACCCTTGAAGAGCAGCAAGAGAACTGGGAACGCGACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATGCGCCGTGAAGTAGCCCGCGATc                                                                                                                            >  1:262488/1‑149 (MQ=255)
                      cGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATc                                                                                                                            >  6:291226/1‑149 (MQ=255)
                      cGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATc                                                                                                                            >  4:29256/1‑149 (MQ=255)
                      cGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAAAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATc                                                                                                                            <  1:310722/149‑1 (MQ=255)
                        aCTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc                                                                                                                          <  8:393556/149‑1 (MQ=255)
                             cGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAc                                                                                                                     >  2:136473/1‑149 (MQ=255)
                                 ccGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGt                                                                                                                 <  8:84262/149‑1 (MQ=255)
                                      ggcCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCgct                                                                                                            >  2:304010/1‑149 (MQ=255)
                                         cTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc                                                                                                                          >  1:391832/1‑132 (MQ=255)
                                         cTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc                                                                                                                          <  2:391832/132‑1 (MQ=255)
                                             ttGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGAt                                                                                                     <  3:290002/149‑1 (MQ=255)
                                              tgtgGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGAt                                                                                                     <  5:269528/148‑1 (MQ=255)
                                               gtgGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATca                                                                                                   >  4:157734/1‑149 (MQ=255)
                                                  gCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAg                                                                                                >  8:272837/1‑149 (MQ=255)
                                                        aatgaTAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTgg                                                                                          <  1:304010/149‑1 (MQ=255)
                                                          tgaTAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGt                                                                                                                 >  3:167940/1‑124 (MQ=255)
                                                          tgaTAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGt                                                                                                                 <  4:167940/124‑1 (MQ=255)
                                                            aTAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCa                                                                                       >  8:236799/1‑148 (MQ=255)
                                                             tAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCata                                                                                     >  1:390680/1‑149 (MQ=255)
                                                               aaCAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCatat                                                                                    <  3:309559/148‑1 (MQ=255)
                                                                aCAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTc                                                                                  <  3:29256/149‑1 (MQ=255)
                                                                aCAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCCGTGCGCTGCTGGATCACAGCGGTGGGCATAtt                                                                                   >  1:312212/1‑148 (MQ=255)
                                                                 cAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCa                                                                                 <  1:136473/149‑1 (MQ=255)
                                                                  aGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCTCAGCGGTGGGCATATTCAt                                                                                <  6:285219/149‑1 (MQ=255)
                                                                  aGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCAt                                                                                >  7:33718/1‑149 (MQ=255)
                                                                   gCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATc                                                                               <  2:390680/149‑1 (MQ=255)
                                                                   gCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTCCTGGATCACAGCGGTGGGCATATTCATc                                                                               >  1:40841/1‑149 (MQ=255)
                                                                    cTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATcc                                                                              <  5:222117/149‑1 (MQ=255)
                                                                       gCCACAGTTGAAGAGCAGACAGAGACCTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCCTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGc                                                                            <  3:273553/148‑1 (MQ=255)
                                                                       gCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCTATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGc                                                                            >  6:73662/1‑148 (MQ=255)
                                                                       gCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCTATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGc                                                                            >  6:68376/1‑148 (MQ=255)
                                                                         cacaCTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTa                                                                          >  8:196255/1‑148 (MQ=255)
                                                                              ttGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCtt                                                                    <  6:344605/149‑1 (MQ=255)
                                                                                gAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCTATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCtt                                                                    <  5:68376/147‑1 (MQ=255)
                                                                                gAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCTATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCtt                                                                    <  5:73662/147‑1 (MQ=255)
                                                                                gAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTg                                                                   <  7:272837/148‑1 (MQ=255)
                                                                                      cAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            >  7:42906/1‑149 (MQ=255)
                                                                                      cAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            <  3:364411/149‑1 (MQ=255)
                                                                                        gAAAGAGAACTGGGAACGCTACTGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCGGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCg                                                          <  2:262488/149‑1 (MQ=255)
                                                                                        gAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCAACCGCTAAACCTTGCCCTGGGc                                                           >  4:15000/1‑148 (MQ=255)
                                                                                        gAAAGAGAACTGGGAACGCTACGGCAATAAACACCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATCCCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAACCCTTGCCATTGGc                                                           >  7:153716/1‑148 (MQ=255)
                                                                                            gagaACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGc                                                      <  5:187845/149‑1 (MQ=255)
                                                                                            gagaACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGAAGCCAACTAAAAAAGCCGTGAAGTAGAAAGCGATCTCTATACCGGTGCGATGCTGGATCACAGCCGTGGGCATATTCATACGCTAAAACTTGCCATTGGCGAAGc                                                      >  3:133247/1‑149 (MQ=255)
                                                                                               aaCTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCgctgc                                                                                                          >  6:377031/1‑94 (MQ=255)
                                                                                               aaCTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCgctgc                                                                                                          <  5:377031/94‑1 (MQ=255)
                                                                                                    ggAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCa                                              >  1:221945/1‑149 (MQ=255)
                                                                                                    ggAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGAATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCa                                              >  3:229328/1‑149 (MQ=255)
                                                                                                         gCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGc                                         <  7:348269/149‑1 (MQ=255)
                                                                                                             cGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTc                                      >  5:35379/1‑148 (MQ=255)
                                                                                                              ggCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            <  1:33499/125‑1 (MQ=255)
                                                                                                              ggCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            <  1:39313/125‑1 (MQ=255)
                                                                                                              ggCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            >  2:33499/1‑125 (MQ=255)
                                                                                                              ggCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTgg                                                            >  2:39313/1‑125 (MQ=255)
                                                                                                                   tAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAAcggcg                               >  2:70774/1‑149 (MQ=255)
                                                                                                                   tAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAAcggcg                               <  8:42906/149‑1 (MQ=255)
                                                                                                                      aCAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACggcgggc                            >  5:292006/1‑149 (MQ=255)
                                                                                                                          cTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCgt                        >  5:113243/1‑149 (MQ=255)
                                                                                                                          cTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCgt                        <  1:70774/149‑1 (MQ=255)
                                                                                                                          cTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCgt                        <  6:35379/149‑1 (MQ=255)
                                                                                                                              aaTTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCTGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTAt                    <  4:133247/149‑1 (MQ=255)
                                                                                                                              aaTTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTAt                    <  3:15000/149‑1 (MQ=255)
                                                                                                                                    tGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACtt              >  3:240572/1‑149 (MQ=255)
                                                                                                                                      gACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAAGCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGTGTGTATGAACTTTc            <  3:157734/149‑1 (MQ=255)
                                                                                                                                       acgccaacgccaTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCt           <  7:236799/149‑1 (MQ=255)
                                                                                                                                         gccaacgccaTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCatat                                                                                    <  7:308060/74‑1 (MQ=255)
                                                                                                                                         gccaacgccaTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCatat                                                                                    >  8:308060/1‑74 (MQ=255)
                                                                                                                                            aacgccaTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGt      >  1:44595/1‑149 (MQ=255)
                                                                                                                                                 caTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACg  >  4:248274/1‑148 (MQ=255)
                                                                                                                                                 caTTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACg  >  8:332234/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCGCTGCTGGATCACAGCGGTGGGCATATTCATCCGCTAAACCTTGCCATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACG  >  NC_000913/1364563‑1364855

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