Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,159,394 A→T 6.3% W585R (TGG→AGG)  ygiQ ← radical SAM superfamily protein YgiQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,159,3940AT 6.3% 195.3 / 3.2 63W585R (TGG→AGG) ygiQradical SAM superfamily protein YgiQ
Reads supporting (aligned to +/- strand):  ref base A (23/36);  new base T (2/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCA  >  NC_000913/3159245‑3159543
                                                                                                                                                     |                                                                                                                                                     
gTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCca                                                                                                                                                       >  2:68829/1‑150 (MQ=255)
  cGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCcaca                                                                                                                                                     >  1:781629/1‑150 (MQ=255)
   gCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAg                                                                                                                                                    >  2:1391656/1‑150 (MQ=255)
   gCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAg                                                                                                                                                    <  1:1391656/150‑1 (MQ=255)
    cccTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAgc                                                                                                                                                   >  1:1097220/1‑150 (MQ=255)
     ccTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGTGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAgcg                                                                                                                                                  >  1:1174458/1‑150 (MQ=255)
          ggTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGa                                                                                                                                             >  1:369017/1‑150 (MQ=255)
          ggTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGa                                                                                                                                             >  2:1404851/1‑150 (MQ=255)
           gTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGAt                                                                                                                                            >  2:492500/1‑150 (MQ=255)
               gAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCAc                                                                                                                                        >  1:921485/1‑150 (MQ=255)
                aaGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCAcc                                                                                                                                       <  2:781629/150‑1 (MQ=255)
                   aCGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCata                                                                                                                                    <  2:498087/150‑1 (MQ=255)
                     gTCTTCGCTCTTATAACCAATCTTCGCCATCGTGTTTTTGCCGTTGTAATACATTGTGTTTGAGTTCTCCAGCGGCGACGTATAGAAGTTCTTCACCTGGTCGAGGCGGAAGCGGTTCTTTTTCAGCCACAGCGCCAGATTCACCATATc                                                                                                                                  <  2:705464/150‑1 (MQ=255)
                      tCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCt                                                                                                                                 >  1:1428498/1‑150 (MQ=255)
                       cTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCtt                                                                                                                                >  1:684995/1‑150 (MQ=255)
                           gCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTtcatc                                                                                                                            >  2:587474/1‑150 (MQ=255)
                              ctTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCAcg                                                                                                                         <  1:328323/150‑1 (MQ=255)
                                 ataACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGt                                                                                                                      >  2:1179667/1‑150 (MQ=255)
                                  taACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTa                                                                                                                     <  1:662144/150‑1 (MQ=255)
                                  taACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTa                                                                                                                     >  1:749773/1‑150 (MQ=255)
                                     ccAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCg                                                                                                                  >  1:150432/1‑150 (MQ=255)
                                     ccAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGCAccc                                                                                                                  >  2:439962/1‑149 (MQ=255)
                                      cAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCgg                                                                                                                 >  1:32641/1‑150 (MQ=255)
                                      cAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCgg                                                                                                                 <  1:1096459/150‑1 (MQ=255)
                                        aTCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCgggg                                                                                                               <  1:835165/150‑1 (MQ=255)
                                            tCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTgcg                                                                                                           >  2:494583/1‑150 (MQ=255)
                                             cGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTgcgc                                                                                                          >  2:1147094/1‑150 (MQ=255)
                                                cAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGga                                                                                                       >  2:416048/1‑150 (MQ=255)
                                                  gCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGgaga                                                                                                     <  2:228562/150‑1 (MQ=255)
                                                  gCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGCACCCCGGTAcgcgtgat                                                                                                     <  1:439962/150‑5 (MQ=255)
                                                    gggTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGGGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATg                                                                                                   <  1:756766/150‑1 (MQ=255)
                                                        ttttCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAAt                                                                                               <  2:1428498/150‑1 (MQ=255)
                                                                 tgtAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCa                                                                                      <  2:518993/150‑1 (MQ=255)
                                                                   tAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCa                                                                                                                                         >  1:801319/1‑97 (MQ=255)
                                                                   tAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCa                                                                                                                                         <  2:801319/97‑1 (MQ=255)
                                                                   tAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGa                                                                                    <  1:943658/150‑1 (MQ=255)
                                                                    aaTACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGAt                                                                                   <  1:1404851/150‑1 (MQ=255)
                                                                     aTACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCTCAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATa                                                                                  >  1:892937/1‑150 (MQ=255)
                                                                     aTACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCTCAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATa                                                                                  >  1:887746/1‑150 (MQ=255)
                                                                      tACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATAc                                                                                 <  1:666087/150‑1 (MQ=255)
                                                                         aTGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGc                                                                              >  2:244662/1‑150 (MQ=255)
                                                                           ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGctc                                                                            >  1:1206575/1‑150 (MQ=255)
                                                                           ggtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGctc                                                                            <  2:749773/150‑1 (MQ=255)
                                                                            gtggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGctct                                                                           <  2:1423291/150‑1 (MQ=255)
                                                                             tggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCtt                                                                          <  2:1420872/150‑1 (MQ=255)
                                                                              ggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCttt                                                                         <  2:1097220/150‑1 (MQ=255)
                                                                              ggtTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCttt                                                                         <  1:549744/150‑1 (MQ=255)
                                                                                     ttCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTAcctgcc                                                                  <  2:400793/150‑1 (MQ=255)
                                                                                     ttCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTAcctgcc                                                                  <  1:1377450/150‑1 (MQ=255)
                                                                                         ccAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGtt                                                              >  2:252839/1‑150 (MQ=255)
                                                                                           aGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGtttt                                                            <  1:1000410/150‑1 (MQ=255)
                                                                                               gcgACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGagt                                                        <  2:1031066/150‑1 (MQ=255)
                                                                                                   cGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGagtaagt                                                    <  2:1058318/150‑1 (MQ=255)
                                                                                                       tAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCg                                                <  1:244662/150‑1 (MQ=255)
                                                                                                          aaGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAAc                                             <  1:300890/150‑1 (MQ=255)
                                                                                                            gTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCTCAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAg                                           <  2:892937/150‑1 (MQ=255)
                                                                                                            gTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCTCAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAg                                           <  2:887746/150‑1 (MQ=255)
                                                                                                                      tGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAg                                 <  2:684995/150‑1 (MQ=255)
                                                                                                                                ggAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGc                       <  1:492500/150‑1 (MQ=255)
                                                                                                                                    gCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGTTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGcc                   <  1:588186/150‑1 (MQ=255)
                                                                                                                                    gCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGcc                   <  2:147510/150‑1 (MQ=255)
                                                                                                                                             tttCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCgg          >  1:3043/1‑150 (MQ=255)
                                                                                                                                                   ccACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtcat    <  1:68829/150‑1 (MQ=255)
                                                                                                                                                   ccACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtcat    <  1:1377352/150‑1 (MQ=255)
                                                                                                                                                   ccACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtcat    <  2:7939/150‑1 (MQ=255)
                                                                                                                                                    cacaGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGtttt                                                            >  1:238714/1‑93 (MQ=255)
                                                                                                                                                    cacaGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGtttt                                                            <  2:238714/93‑1 (MQ=255)
                                                                                                                                                    cacaGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtcatc   <  2:150432/150‑1 (MQ=255)
                                                                                                                                                     acaGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTtcatca  <  2:1179224/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCA  >  NC_000913/3159245‑3159543

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