Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,170,795 G→A A1057V (GCG→GTG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,7950GA98.4% 226.6 / NA 61A1057V (GCG→GTG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (37/23);  minor base C (0/1);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGA  >  NC_000913/1170656‑1170937
                                                                                                                                           |                                                                                                                                              
caAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCt                                                                                                                                        <  1:608310/148‑1 (MQ=255)
  aaCCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGgcag                                                                                                                                   >  1:773397/1‑151 (MQ=255)
    ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCg                                                                                                                                             >  1:1030609/1‑139 (MQ=255)
    ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCg                                                                                                                                             <  2:1030981/139‑1 (MQ=255)
    ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGgcagca                                                                                                                                 >  1:313503/1‑151 (MQ=255)
    ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGgcagca                                                                                                                                 >  1:360471/1‑151 (MQ=255)
                 ccGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                                                    <  2:285822/151‑1 (MQ=255)
                 ccGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                                                    <  1:519545/151‑1 (MQ=255)
                    gATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGcgatcgat                                                                                                                 >  2:396921/1‑151 (MQ=255)
                     aTTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATa                                                                                                                <  1:322398/151‑1 (MQ=255)
                         aCGTGATTCTTCTCGGCAAATTCGATCTCCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCt                                                                                                            <  1:558368/151‑1 (MQ=255)
                         aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAg                                                                                                              >  1:193455/1‑149 (MQ=255)
                         aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCt                                                                                                            <  1:990410/151‑1 (MQ=255)
                          cGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCt                                                                                                            >  2:918716/1‑150 (MQ=255)
                              aTTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGa                                                                                                         >  2:35257/1‑149 (MQ=255)
                               ttcttcTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACtt                                                                                                      <  1:62741/151‑1 (MQ=255)
                               tccttctCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACtt                                                                                                      <  1:978512/149‑1 (MQ=255)
                                   tctcGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAt                                                                                                  >  2:461533/1‑151 (MQ=255)
                                      cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAATCTCGACTTTGAtct                                                                                                >  1:902932/1‑150 (MQ=255)
                                        gCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAAACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTGGATctct                                                                                              <  2:885361/150‑1 (MQ=255)
                                             ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGt                                                                                        >  2:65173/1‑151 (MQ=255)
                                             ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGt                                                                                        >  1:807439/1‑151 (MQ=255)
                                             ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGt                                                                                        >  1:412146/1‑151 (MQ=255)
                                               cGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGtt                                                                                       >  2:577558/1‑150 (MQ=255)
                                                         gCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGtt                                                                            >  2:34676/1‑151 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  1:794867/1‑150 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  2:424660/1‑150 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  1:318530/1‑150 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  2:522256/1‑150 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  2:1058284/1‑150 (MQ=255)
                                                              tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTgg                                                                        >  2:393918/1‑150 (MQ=255)
                                                                          ccAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                                                           <  1:3164/151‑1 (MQ=255)
                                                                          ccAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                                                           <  1:834766/151‑1 (MQ=255)
                                                                          ccAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                                                           <  1:925662/151‑1 (MQ=255)
                                                                          ccAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                                                           <  1:495789/151‑1 (MQ=255)
                                                                          ccAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                                                           <  2:519014/151‑1 (MQ=255)
                                                                                 ccTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAg                                                    <  1:837720/151‑1 (MQ=255)
                                                                                        cccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACa                                             >  1:160306/1‑151 (MQ=255)
                                                                                            gTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTCTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcg                                         >  2:989557/1‑151 (MQ=255)
                                                                                             tttCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcg                                         >  2:792288/1‑150 (MQ=255)
                                                                                              ttCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgcg                                       <  1:605067/151‑1 (MQ=255)
                                                                                                  gcgcTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGtt                                   >  1:781589/1‑151 (MQ=255)
                                                                                                       tGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACTACAGACGCGTGTTCACGt                              >  1:606885/1‑151 (MQ=255)
                                                                                                         cTGGCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTc                             <  2:683331/150‑1 (MQ=255)
                                                                                                           ggCGCAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAgg                          >  1:460785/1‑151 (MQ=255)
                                                                                                            gcgcAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAggg                         >  2:496289/1‑151 (MQ=255)
                                                                                                             cgcAGACGGGCAATATCCAGCAGGGTACGCCCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGa                        <  2:349889/151‑1 (MQ=255)
                                                                                                             cgcAGACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTCCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGa                        <  2:518968/151‑1 (MQ=255)
                                                                                                                 gACGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGaa                    >  2:406000/1‑151 (MQ=255)
                                                                                                                  aCGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGGCTAAGAGATCGATAAGCTCGACTTTGATCTCCTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGaaa                   <  2:219745/151‑1 (MQ=255)
                                                                                                                   cGGGCAATATCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAat                  <  1:461349/151‑1 (MQ=255)
                                                                                                                         atatCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGCTAAGCTCGACTTTGTTCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGAAGAACGACAGACGCGTGTTCACGTCAGGGATGAAatcat               >  2:211123/1‑148 (MQ=255)
                                                                                                                         atatCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAatcat               >  1:492129/1‑148 (MQ=255)
                                                                                                                         atatCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAatcat               >  2:241779/1‑148 (MQ=255)
                                                                                                                           atCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTgg           >  2:805787/1‑150 (MQ=255)
                                                                                                                           atCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTgg           >  2:313567/1‑150 (MQ=255)
                                                                                                                           atCCAGCAGGGTACGCACCGGATCTGGCAGCAGGCCCAAGCGATCGATAAGCTCGACTTTGCTCTCTTCCAGTTCGTTTTCCGTTTTGTCGCTGGCAATACGTTACTAGAACGACCGACGCGTGTTCACGTCAGGGATGAAATCACCTgg           >  1:984250/1‑150 (MQ=255)
                                                                                                                                gcagGGTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGAAATAg     >  2:237637/1‑151 (MQ=255)
                                                                                                                                    ggTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGc    <  1:460983/148‑1 (MQ=255)
                                                                                                                                    ggTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGc    <  1:1020808/148‑1 (MQ=255)
                                                                                                                                    ggTACGCACCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGa  >  2:228584/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGA  >  NC_000913/1170656‑1170937

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