Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 19.9% intergenic (+24/+79) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3110CG19.9% 333.9 / 80.5 143intergenic (+24/+79)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (50/63);  new base G (17/11);  total (68/75)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 147 total aligned reads displayed.

GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCA  >  NC_000913/2286212‑2286406
                                                                                                   |                                                                                                
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTGAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                  >  2:882489/1‑100 (MQ=255)
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                  >  2:368924/1‑100 (MQ=255)
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                  <  1:2945807/100‑1 (MQ=255)
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                  <  1:1477609/100‑1 (MQ=255)
  cACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCtt                                                                                                >  1:1618321/1‑100 (MQ=255)
     gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                              <  2:539/99‑1 (MQ=255)
     gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                             >  1:2834504/1‑100 (MQ=255)
      tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgg                                                                                            <  1:2223218/100‑1 (MQ=255)
       aGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggt                                                                                           >  1:1136388/1‑100 (MQ=255)
         cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                         <  2:3227065/100‑1 (MQ=255)
         cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTGCAAATCCTCTCGTGCCGACCAACAATCCCAAGAAAAAACCAACCGTAAGggttg                                                                                         >  1:1166319/1‑100 (MQ=255)
          gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                        >  1:1379064/1‑100 (MQ=255)
          gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                        >  1:2582350/1‑100 (MQ=255)
          gcAGCCTGGTAGCGCACCGTCAGGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgggt                                                                                          >  2:2113696/1‑96 (MQ=255)
          gcAGCCTGGTAGCGCACCGTCAAGGGGCGTCGGGGGTCGGAGGCTCAAACCCCCTCGTGCCGACAAAAACTCCCAAGAAAAAACCCACCCTTACggttgg                                                                                        >  2:2116738/1‑100 (MQ=255)
            tccgatctAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggt‑t                                                                                     <  2:466884/93‑1 (MQ=255)
            aGCCAGGTAGCGAACCGTCATAGGGAGTCGGGGGTCGGAGGTTCAAATCCTCACGTGCCGAGCAAAAACCCCTAGAAAAAACCAACCCTTACggn                                                                                           <  2:2153246/95‑2 (MQ=255)
             gCCTGGTAGCGCACCGTCCTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGa‑tt                                                                                    >  2:2646914/1‑100 (MQ=255)
                tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGt‑ttttt                                                                                 <  1:732592/100‑1 (MQ=255)
                tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGt‑ttttt                                                                                 <  1:742780/100‑1 (MQ=255)
                 ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGt‑tttttt                                                                                >  2:1047579/1‑100 (MQ=255)
                 ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGt‑tttttt                                                                                <  2:625645/100‑1 (MQ=255)
                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGt‑tttttct                                                                               <  1:2245738/100‑3 (MQ=38)
                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTta                                                                               <  2:1136379/100‑1 (MQ=255)
                   tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTtat                                                                              >  1:1947822/1‑100 (MQ=255)
                   tAGAGCGCCTTCATGGGGTATCGGAGGTCGGATGGTCAAATACTCGCGTGGCGATCAAAAATCCCAAGAAGAAACCAACCCTTACGGTTGGt‑tttttt                                                                                <  2:624035/98‑1 (MQ=255)
                    aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggt‑tagatcgga                                                                             >  1:466886/1‑92 (MQ=255)
                    aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggt‑tagatcgg                                                                              >  1:2295517/1‑92 (MQ=38)
                     gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGt‑tttttcttgg                                                                            <  1:1588248/100‑6 (MQ=37)
                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCt                                                                          <  2:1347839/100‑1 (MQ=255)
                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCt                                                                          >  1:25912/1‑100 (MQ=255)
                       gcACAGTCATGGGGTCCCGGGAGTCGCAGACTCAAATCCCCTCGTTCCGACCAAAAATCCCAAGAAAACACCAACCGTAAGGGTTAGt‑ttcttct                                                                               <  2:2358583/95‑3 (MQ=255)
                        cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTg                                                                         >  1:2071462/1‑100 (MQ=255)
                            gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaat                                                                     >  1:1671570/1‑100 (MQ=255)
                            gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaat                                                                     >  1:1122109/1‑100 (MQ=255)
                             ttaTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaatt                                                                    <  2:2350775/98‑1 (MQ=255)
                             tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaatt                                                                    <  2:2532128/100‑1 (MQ=255)
                              cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaatta                                                                   <  1:3499551/100‑1 (MQ=255)
                                tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaattaat                                                                 >  1:3074380/1‑100 (MQ=255)
                                tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGA‑TTTTTTATATCTGCaattaat                                                                 >  2:3414466/1‑100 (MQ=255)
                                 ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCaattaatt                                                                >  1:1314623/1‑100 (MQ=255)
                                  gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTc                                                               >  1:493597/1‑100 (MQ=255)
                                  gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTc                                                               <  1:2681223/100‑1 (MQ=255)
                                         gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAAc                                                        <  1:1993652/100‑1 (MQ=255)
                                          ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACa                                                       <  1:1964111/100‑1 (MQ=255)
                                                ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                 >  2:3100338/1‑100 (MQ=255)
                                                ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                 >  1:1309042/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                              >  1:85385/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACACACCGTGac                                              >  2:2772916/1‑100 (MQ=255)
                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                            >  1:807851/1‑100 (MQ=255)
                                                      taaaaTCCTCTCGTGCCTACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                              <  2:419631/95‑1 (MQ=255)
                                                       cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAGCCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                          >  2:1949293/1‑100 (MQ=255)
                                                       cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAACCCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                          >  2:1949276/1‑100 (MQ=255)
                                                        aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                         >  2:3151236/1‑100 (MQ=255)
                                                        aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                         <  2:2834484/100‑1 (MQ=255)
                                                         aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATca                                        <  2:2618630/100‑1 (MQ=255)
                                                          aTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcac                                       <  2:2059015/100‑1 (MQ=255)
                                                             ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                    <  1:2960263/100‑1 (MQ=255)
                                                                tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                 <  2:85385/100‑1 (MQ=255)
                                                                  gTGCCGACCAACAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                               <  1:1699664/100‑1 (MQ=255)
                                                                  gTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                               >  1:2005328/1‑100 (MQ=255)
                                                                   tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                     <  2:942202/93‑1 (MQ=255)
                                                                     ccGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                            >  1:1872787/1‑100 (MQ=255)
                                                                        aCCAACAATCCCAAGAAAAAACCAACCATTACGGTTGGT‑TTTTTTATATCTGCAATTAATTAGATAAACAAACAGTGACACATCACAATCTGTTTAt                            >  2:2416143/1‑97 (MQ=255)
                                                                         ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                        >  1:355921/1‑100 (MQ=255)
                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                       <  2:1370817/100‑1 (MQ=255)
                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                          <  2:2724468/97‑1 (MQ=255)
                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                      <  2:1504987/100‑1 (MQ=255)
                                                                            aaaaTCCCAGGAAAAGACCTACCCTTCCGGTTGGT‑TGTTTTATAACTGGAATTATTTCGATAAACAGGCCGTGACACATCACAGCCTGTTTACTTTCTg                      <  2:3268450/99‑1 (MQ=255)
                                                                            aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                     <  1:65891/100‑1 (MQ=255)
                                                                             aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                    <  1:916709/100‑1 (MQ=255)
                                                                               aTCCCAAGAAAAAACCAACCCTTACTGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAt                  <  2:93675/100‑1 (MQ=255)
                                                                                tCCCAAGAAAAAACCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                 <  1:851738/100‑1 (MQ=38)
                                                                                agatataaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                 >  1:14958/8‑100 (MQ=37)
                                                                                 cccAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                <  1:3098130/100‑1 (MQ=255)
                                                                                 cccAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                <  1:3238/100‑1 (MQ=255)
                                                                                  ccAAGAAACAACCAACCCTTACGGCTGGT‑TTTTTTATATCTGCAATTAATTCGATATACAGACCGTGACTCATCATAGCATGTTTATTTTCTGTTATCAg               >  2:1657767/1‑100 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg               >  2:1409463/1‑100 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg               >  1:118651/1‑100 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              >  1:1394344/1‑100 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              <  1:1926471/100‑1 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              >  2:1889431/1‑100 (MQ=255)
                                                                                    ataaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa             >  2:2815370/4‑100 (MQ=38)
                                                                                     taaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAGACAGACCGTTACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  2:301490/3‑100 (MQ=38)
                                                                                     taaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  1:593702/3‑100 (MQ=38)
                                                                                     taaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  1:1890527/3‑100 (MQ=38)
                                                                                      naaaaaaCCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTTTTTTCTGTTATCAGAACg           >  1:3265918/2‑100 (MQ=255)
                                                                                      gAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg           >  1:3108356/1‑100 (MQ=255)
                                                                                      gAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg           <  2:1169231/100‑1 (MQ=255)
                                                                                      aaaaaaaCCAACCGTCAGGGGTGGT‑TTTTTTATATCTGCAATCAATACGATAAACAGACCGTGACCCCTCCCAGCCTGTTTATTTTCTGTAAt                  >  2:740945/2‑93 (MQ=25)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  2:796510/2‑99 (MQ=38)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg           >  1:320910/2‑100 (MQ=38)
                                                                                       aaaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  2:1457232/1‑98 (MQ=255)
                                                                                        aaaaaCCAACCGTAAGGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc         >  1:2059456/1‑100 (MQ=39)
                                                                                         aaaaCCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc        >  1:2700724/1‑100 (MQ=255)
                                                                                            tctaaccCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGa     <  1:1199012/97‑1 (MQ=255)
                                                                                             ccaaccCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc    >  2:1383336/1‑100 (MQ=255)
                                                                                               aaccCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca  <  1:2428936/100‑1 (MQ=255)
                                                                                               aaccCTTACGGTTGGT‑ATTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc    <  2:3467623/98‑1 (MQ=255)
                                                                                                accTTTACGGTTGGT‑TTTTTTTTATCTGCAATTAATTAGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca  <  2:2280387/99‑1 (MQ=255)
                                                                                                   |                                                                                                
GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGT‑TTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCA  >  NC_000913/2286212‑2286406

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