Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,417 C→T 6.2% G42S (GGC→AGC)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,4170CT 6.2% 198.3 / 3.0 64G42S (GGC→AGC) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/32)
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.41e-01

ACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTC  >  NC_000913/2147269‑2147565
                                                                                                                                                    |                                                                                                                                                    
aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                                                                                      <  7:249524/149‑1 (MQ=255)
   aTTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAc                                                                                                                                                   <  4:187136/149‑1 (MQ=255)
    tttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAcc                                                                                                                                                  <  5:308327/149‑1 (MQ=255)
    tttGGGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAcc                                                                                                                                                  >  5:85165/1‑149 (MQ=255)
               cATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTAc                                                                                                                                        <  4:321550/148‑1 (MQ=255)
                aTTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACgg                                                                                                                                      <  1:192449/149‑1 (MQ=255)
                 ttGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTGCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACggg                                                                                                                                     <  7:139422/149‑1 (MQ=255)
                     aGATCACACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAACCCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTTTATCCGGAATAGCAGTCACCCCGCCGCGCTATTCGCCCCCCGCCAGACTACGGGCat                                                                                                                                  >  4:257991/1‑148 (MQ=255)
                          aaactcgcggcTCATTTTCGCCATACAGTCTGCGGCAACAGGTTCTATTCCTGCACTTAACTTTATGTGCAGAGTATGGCCGGCTATATCCTGAATTGCAGTCAACTCGCCGCGTTATTCGCCCACCGCCTGTGTACGGGCAtattaa                                                                                                                             <  6:260209/142‑1 (MQ=255)
                          aaaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                                                                                                                             <  5:296846/148‑1 (MQ=255)
                            aCAGGCGACTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTg                                                                                                                          >  6:200284/1‑149 (MQ=255)
                              aGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCTCACCGCCAGACTACGGGCATAATAACTGTc                                                                                                                        >  1:201616/1‑149 (MQ=255)
                                     gcATTTTCGCCAGACACTCTGCGGCAACAGGTTCAAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAcc                                                                                                                 <  1:212128/148‑1 (MQ=255)
                                           tCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCTCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTc                                                                                                            <  2:201616/148‑1 (MQ=255)
                                             gCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCAGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCca                                                                                                          <  8:353019/148‑1 (MQ=255)
                                                  acacTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGc                                                                                                    >  7:146007/1‑149 (MQ=255)
                                                      tctGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGCCTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCAc                                                                                                >  6:152949/1‑149 (MQ=255)
                                                      tctGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCAc                                                                                                >  8:284140/1‑149 (MQ=255)
                                                      tctGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCAc                                                                                                >  4:210430/1‑149 (MQ=255)
                                                           ggCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCAc                                                                                           <  5:200284/149‑1 (MQ=255)
                                                             cAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACgg                                                                                         >  3:142396/1‑149 (MQ=255)
                                                                cAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcg                                                                                      >  7:133117/1‑149 (MQ=255)
                                                                        aaaCCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGCCTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaaa                                                                              >  2:306668/1‑149 (MQ=255)
                                                                         aaCCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaaa                                                                              <  7:284140/148‑1 (MQ=255)
                                                                         aaCCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAAt                                                                             <  8:133117/149‑1 (MQ=255)
                                                                             tGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc                                                                          >  5:137885/1‑148 (MQ=255)
                                                                               cACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                         <  6:137885/147‑1 (MQ=255)
                                                                                    aaaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCa                                                                                            <  7:237215/123‑1 (MQ=255)
                                                                                    aaaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCa                                                                                            >  8:237215/1‑123 (MQ=255)
                                                                                     aaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaa                                                                               >  4:180317/1‑135 (MQ=255)
                                                                                     aaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaa                                                                               <  3:180317/135‑1 (MQ=255)
                                                                                     aaTTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCAc                                                                 >  3:60077/1‑149 (MQ=255)
                                                                                           tgtgCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGTCGAGAAATCCCTACATCCACGACCa                                                            <  5:152949/148‑1 (MQ=255)
                                                                                             tgCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGt                                                          <  8:177750/148‑1 (MQ=255)
                                                                                               cAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCCCGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCa                                                        >  1:307543/1‑148 (MQ=255)
                                                                                               cAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCa                                                        <  4:142396/148‑1 (MQ=255)
                                                                                                agagTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATa                                                      >  1:201561/1‑149 (MQ=255)
                                                                                                 gagTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAc                                                     >  8:335960/1‑149 (MQ=255)
                                                                                                         cGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                             >  1:111928/1‑149 (MQ=255)
                                                                                                              tataTCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGt                                        >  6:164440/1‑149 (MQ=255)
                                                                                                                  tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTTATACGGCGGCTGCCAGTTCAg                                    >  2:357619/1‑149 (MQ=255)
                                                                                                                   ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAgg                                   >  2:385691/1‑149 (MQ=255)
                                                                                                                    cGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAgg                                   <  6:157922/148‑1 (MQ=255)
                                                                                                                           gCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACa                            <  2:111928/148‑1 (MQ=255)
                                                                                                                           gCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACAt                           >  8:232600/1‑149 (MQ=255)
                                                                                                                              gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCaa                                                                        >  5:135538/1‑101 (MQ=255)
                                                                                                                              gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCaa                                                                        <  6:135538/101‑1 (MQ=255)
                                                                                                                               tCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGATTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                       >  4:101339/1‑149 (MQ=255)
                                                                                                                                   cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                              >  2:121340/1‑106 (MQ=255)
                                                                                                                                   cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                              <  1:121340/106‑1 (MQ=255)
                                                                                                                                   cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAGCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAgctc                   >  6:56731/1‑149 (MQ=255)
                                                                                                                                   cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtctc                   >  4:151627/1‑149 (MQ=255)
                                                                                                                                   cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtctc                   >  6:145388/1‑149 (MQ=255)
                                                                                                                                       ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTa               <  2:307543/149‑1 (MQ=255)
                                                                                                                                       ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTa               <  2:263542/149‑1 (MQ=255)
                                                                                                                                        cgcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              <  1:385691/149‑1 (MQ=255)
                                                                                                                                        cgcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              <  8:146007/149‑1 (MQ=255)
                                                                                                                                          cgATATTCGCTCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc                                                                          >  5:385505/1‑87 (MQ=255)
                                                                                                                                          cgATATTCGCTCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc                                                                          <  6:385505/87‑1 (MQ=255)
                                                                                                                                          cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt             >  7:279942/1‑148 (MQ=255)
                                                                                                                                            atatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt          <  3:210430/149‑1 (MQ=255)
                                                                                                                                            ataaTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGATATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt          <  6:85165/149‑1 (MQ=255)
                                                                                                                                             tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAACTCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTACTTATTAATc         >  5:104281/1‑149 (MQ=255)
                                                                                                                                               ttCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTTATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc        <  1:357619/148‑1 (MQ=255)
                                                                                                                                                  gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCt                 <  7:258909/136‑1 (MQ=255)
                                                                                                                                                    ccACCGCCAGACTACGGGGATAATAACTGTACGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGTTATACATCGCATCTCCTTATTCATCCGCTTTc  <  5:284299/149‑1 (MQ=255)
                                                                                                                                                    ccACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCttt   >  3:350204/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTC  >  NC_000913/2147269‑2147565

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: