Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,508,682 T→A 15.6% pseudogene (261/267 nt) insM ← pseudogene, IS600 transposase;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,508,6820TA15.6% 274.2 / 38.6 109pseudogene (261/267 nt)insMpseudogene, IS600 transposase;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
Reads supporting (aligned to +/- strand):  ref base T (44/48);  new base A (6/11);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 113 total aligned reads displayed.

GGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGG  >  NC_000913/4508585‑4508780
                                                                                                 |                                                                                                  
ggCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGc                                                                                                  >  2:2269228/1‑100 (MQ=32)
 gCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCt                                                                                                 >  1:2509777/1‑100 (MQ=33)
  cTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCtt                                                                                                <  1:1377379/100‑1 (MQ=35)
    ttttCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCtttt                                                                                              >  2:949035/1‑100 (MQ=38)
      ttCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCtttttt                                                                                            <  1:102135/100‑1 (MQ=38)
       tCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCttttttt                                                                                           <  1:2361726/100‑1 (MQ=38)
         ctctGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTaa                                                                                         >  2:881879/1‑100 (MQ=38)
          tctGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAat                                                                                        <  2:1244448/100‑1 (MQ=38)
          tctGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAat                                                                                        <  1:224174/100‑1 (MQ=38)
          tatgctgTAAATGTTCGTATCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTa                                                                                          <  1:2015326/96‑1 (MQ=255)
              ctgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAatatat                                                                                    >  2:2264163/1‑100 (MQ=14)
              ctgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAAACCCATGTTTTCACATAACTGTTGCGTTGACCAATTGACTCTACAGTTGCTTTTTTTAatatat                                                                                    >  1:2555156/1‑100 (MQ=255)
               tgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtc                                                                                   >  1:916190/1‑100 (MQ=40)
               tgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtc                                                                                   >  1:1217310/1‑100 (MQ=255)
               tgTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtc                                                                                   <  2:2233424/100‑1 (MQ=255)
                gTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                                                  >  1:2269394/1‑100 (MQ=255)
                gTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtc                                                                                   >  2:1574985/1‑99 (MQ=40)
                gTAAATGTTCGTCTCATTCAAAAACCTCCGCAATCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtctc                                                                                 <  2:353119/100‑1 (MQ=255)
                 tAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATAtctc                                                                                 <  2:626862/100‑1 (MQ=25)
                  aaaTGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTTCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCg                                                                                <  1:948199/100‑1 (MQ=255)
                   aaTGTTCTTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGc                                                                               <  1:1458398/100‑1 (MQ=255)
                   aaTGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGc                                                                               <  2:822369/100‑1 (MQ=255)
                   aaTGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGc                                                                               <  1:1647067/100‑1 (MQ=255)
                   aaTGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACAAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGc                                                                               <  2:1100842/100‑1 (MQ=255)
                      gTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCa                                                                            <  2:2572074/100‑1 (MQ=255)
                       ttCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCaa                                                                           >  1:633378/1‑100 (MQ=255)
                        tCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAg                                                                          <  2:2041177/100‑1 (MQ=255)
                         cGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAg                                                                          <  2:585324/99‑1 (MQ=255)
                           tctcATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCg                                                                       >  2:81899/1‑100 (MQ=255)
                           tctcATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCCACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCg                                                                       <  2:1437773/100‑1 (MQ=255)
                            ctcATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGa                                                                      <  2:1977509/100‑1 (MQ=255)
                              cATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGc                                                                    >  1:71031/1‑100 (MQ=255)
                              cATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTCTTTTAATATATCTCGCTCAAGGCGa                                                                      >  1:825347/1‑98 (MQ=255)
                               aTTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCt                                                                   >  2:572313/1‑100 (MQ=255)
                               aTTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTAGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATCTATCTCTCTCAAGGCgtg                                                                     >  1:922615/1‑96 (MQ=17)
                               aTTCAAAAAACCTCAGCAACCACATGTTTTCACATAACTGCTGTGTTGACCAATAGAATATACAGTTGCTTTTTTTAATATGTCTCGCTCAACGCCAgat                                                                   >  1:1170166/1‑98 (MQ=14)
                               aTTCAAAAAACCGCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCt                                                                   >  2:1699304/1‑100 (MQ=255)
                                ttCAAAAATCCTCCGCAACCCCATGTTTTCACATAACTGGTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCt                                                                   <  2:2098699/99‑1 (MQ=255)
                                ttCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCtt                                                                  <  2:916193/100‑1 (MQ=255)
                                  cAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCa                                                                <  2:1687012/100‑1 (MQ=255)
                                  cAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCa                                                                >  1:2587698/1‑100 (MQ=255)
                                   aaaaaaCCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAt                                                               <  1:566472/100‑1 (MQ=255)
                                   aaaaaaCCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAt                                                               >  2:2572594/1‑100 (MQ=255)
                                   aaaaaaCCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAt                                                               >  1:1314622/1‑100 (MQ=255)
                                     aaaaCCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAttt                                                             <  2:2124009/100‑1 (MQ=255)
                                       aaCCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTaa                                                           >  2:1474573/1‑100 (MQ=255)
                                         ccTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAAc                                                          <  2:1539489/99‑1 (MQ=255)
                                         ccTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACg                                                         >  2:2084744/1‑100 (MQ=255)
                                         ccTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACg                                                         <  2:1956451/100‑1 (MQ=255)
                                           tCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTATTATATCTCGCTCAAGGCGAGCTTCATTTAACGcc                                                       >  1:388182/1‑100 (MQ=255)
                                              gCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTa                                                    <  1:395691/100‑1 (MQ=255)
                                               cAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTAc                                                   <  1:698817/100‑1 (MQ=255)
                                               cAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTAc                                                   <  2:737819/100‑1 (MQ=255)
                                               cAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTAc                                                   >  2:720336/1‑100 (MQ=255)
                                                 aCCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGc                                                 >  2:291709/1‑100 (MQ=255)
                                                 aCCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACg                                                  >  1:2486144/1‑99 (MQ=255)
                                                 aCCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGc                                                 >  1:1574423/1‑100 (MQ=255)
                                                  ccccATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCa                                                <  1:2084730/100‑1 (MQ=255)
                                                  ccccATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCa                                                <  2:1626480/100‑1 (MQ=255)
                                                   cccATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTAGCTTTTTTTTTTTTTTTTTTCTCAAAACGACCCCCATTTAAAGCCTTACgta                                                >  1:1943916/1‑97 (MQ=2)
                                                   cccATGTTTTCACATAACTGTTGCGTTGACAAATTGAATATACAGTTGCTTTTTTTAATATATCTCGTTCAAGGCGAGCTTCATTTAATGCCTTACGCAg                                               >  1:1043830/1‑100 (MQ=25)
                                                    ccATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGt                                              >  2:337090/1‑100 (MQ=255)
                                                     cATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGtt                                             >  2:174778/1‑100 (MQ=255)
                                                      aTGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTg                                            <  1:1138375/100‑1 (MQ=255)
                                                      aTGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTg                                            <  2:2139288/100‑1 (MQ=255)
                                                              acaTAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTTTTTTTTTTCGATCAAGGCAAGCTTCATTTAATGTCTTACGCAGTTg                                            >  1:128085/1‑92 (MQ=18)
                                                               caTAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTc                                   <  2:706480/100‑1 (MQ=255)
                                                                aTAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCa                                  <  1:945738/100‑1 (MQ=255)
                                                                 tAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAg                                 >  2:347094/1‑100 (MQ=255)
                                                                 tAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAg                                 <  1:1402870/100‑1 (MQ=255)
                                                                 tAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTAAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAg                                 >  2:236103/1‑100 (MQ=255)
                                                                    cTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTAAGAt                               >  1:1494905/1‑99 (MQ=255)
                                                                      gTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTcc                            >  1:1024004/1‑100 (MQ=255)
                                                                      gTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTcc                            <  1:1763715/100‑1 (MQ=255)
                                                                       ttGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCa                           >  1:93420/1‑100 (MQ=255)
                                                                        tGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAg                          <  2:2468439/100‑1 (MQ=255)
                                                                         gCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGt                         <  1:1650720/100‑1 (MQ=255)
                                                                         gCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGt                         >  2:570888/1‑100 (MQ=255)
                                                                         gCGTTGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGt                         <  1:1060020/100‑1 (MQ=255)
                                                                           gTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTc                       >  1:766308/1‑100 (MQ=255)
                                                                            ttGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGATTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCa                      >  1:482648/1‑100 (MQ=255)
                                                                             tGACCAATTGAATCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTCAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAg                     <  1:63117/100‑1 (MQ=255)
                                                                              gACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGc                    >  1:231975/1‑100 (MQ=255)
                                                                                  aaTTGAATCTACAGTTGCTTTTTTTAATATATCTCTCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCAcc                >  1:1526438/1‑100 (MQ=255)
                                                                                        aTCTACAGTTGCTTTTTTTTTTATTTTTTGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCCCCGTGc            >  1:1054054/1‑98 (MQ=255)
                                                                                        aTCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCgg          <  2:1554648/100‑1 (MQ=255)
                                                                                        aTCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCgg          >  2:51872/1‑100 (MQ=255)
                                                                                        aTCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCgg          >  1:356697/1‑100 (MQ=255)
                                                                                        aTCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCgg          <  2:87551/100‑1 (MQ=255)
                                                                                        aTCTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCgg          >  2:179137/1‑100 (MQ=255)
                                                                                         tCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCggg         <  2:221768/100‑1 (MQ=255)
                                                                                          cTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCggg         <  1:2062969/99‑1 (MQ=255)
                                                                                          cTACAGTAGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGa        <  1:280651/100‑1 (MQ=255)
                                                                                            aCAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGaga      >  2:2062982/1‑98 (MQ=255)
                                                                                             cAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAAcc     <  2:971130/100‑1 (MQ=255)
                                                                                             cAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAAcc     <  1:321867/100‑1 (MQ=255)
                                                                                              aGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGTTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCa    >  2:2065884/1‑100 (MQ=255)
                                                                                                ttGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAgg  <  1:358561/100‑1 (MQ=255)
                                                                                                 tGCTTTTTTTTTTTTTTTTCGCTCAAGGCGAGCTTCATTTAACGTCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGttttt          >  1:1953041/1‑87 (MQ=255)
                                                                                                 tGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCACAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAgg  >  1:625892/1‑99 (MQ=255)
                                                                                                 |                                                                                                  
GGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTTCAGATTCCAGTTCAGCCACCGTGCGGGAACCAGG  >  NC_000913/4508585‑4508780

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