Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 20.2% intergenic (+28/+75) 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,3150CG20.2% 317.1 / 79.1 139intergenic (+28/+75)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/61);  new base G (17/11);  total (67/72)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 141 total aligned reads displayed.

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCA  >  NC_000913/2286217‑2286406
                                                                                                  |                                                                                            
gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                              <  2:539/99‑1 (MQ=255)
gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                              <  2:683532/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:2582350/1‑100 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCAGGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgggt                                                                                          >  2:2113696/1‑96 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCAAGGGGCGTCGGGGGTCGGAGGCTCAAACCCCCTCGTGCCGACAAAAACTCCCAAGAAAAAACCCACCCTTACggttgg                                                                                        >  2:2116738/1‑100 (MQ=255)
       aGCCAGGTAGCGAACCGTCATAGGGAGTCGGGGGTCGGAGGTTCAAATCCTCACGTGCCGAGCAAAAACCCCTAGAAAAAACCAACCCTTACggn                                                                                           <  2:2153246/95‑2 (MQ=255)
        gCCTGGTAGCGCACCGTCCTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGatt                                                                                    >  2:2646914/1‑100 (MQ=255)
        gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttt                                                                                    >  1:3635012/1‑100 (MQ=255)
         ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                   <  2:2814279/100‑1 (MQ=39)
         ccTGGTAGAGCACCGGCAAGGGGTGTCGGGGGTCGGAGGTGCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                    <  2:2838206/99‑1 (MQ=37)
           tGGTAGCGCACCGTCATGGGGTGTCGGGTGGCGGAGGTTAAAATCCTCTCGTGCCGACCAAACATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttttt                                                                                 <  2:2413583/100‑1 (MQ=37)
           tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                 <  1:742780/100‑1 (MQ=255)
           tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                 <  1:732592/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCGAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                <  1:3456491/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                <  1:3001199/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                >  2:1047579/1‑100 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                <  2:625645/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                <  1:2286944/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                               <  2:1136379/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                               <  2:595501/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                               <  1:1933964/100‑1 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggttagatcgga                                                                             >  1:466886/1‑92 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggttagatcgg                                                                              >  1:2295517/1‑92 (MQ=38)
               aGCGAACCGTCATGGGGTGTTAGGGGTCGGAGGTGCAAATCCTCTCATGACGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttt                                                                                  <  2:2594403/95‑1 (MQ=25)
                gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                            <  1:2012172/100‑6 (MQ=37)
                gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                            <  1:1588248/100‑6 (MQ=37)
                 cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTatatn                                                                           <  1:1731452/100‑2 (MQ=255)
                  gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                          >  1:25912/1‑100 (MQ=255)
                  gcACAGTCATGGGGTCCCGGGAGTCGCAGACTCAAATCCCCTCGTTCCGACCAAAAATCCCAAGAAAACACCAACCGTAAGGGTTAGtttcttct                                                                               <  2:2358583/95‑3 (MQ=255)
                   cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTg                                                                         >  1:2071462/1‑100 (MQ=255)
                       gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                     >  1:1122109/1‑100 (MQ=255)
                       gTCATGGGGTGTCGGGAGTCGGATGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                       <  2:822007/98‑1 (MQ=255)
                        ttaTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                    <  2:2350775/98‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                    <  2:2532128/100‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                    >  2:2681205/1‑100 (MQ=255)
                         cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                   <  1:3499551/100‑1 (MQ=255)
                           tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGATTTTTTATATCTGCaattaat                                                                 >  2:3414466/1‑100 (MQ=255)
                            ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                >  1:1314623/1‑100 (MQ=255)
                            ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                >  1:1504995/1‑100 (MQ=255)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                               <  1:2681223/100‑1 (MQ=255)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                               <  2:1379056/100‑1 (MQ=255)
                                     ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                       <  1:1964111/100‑1 (MQ=255)
                                      gggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAg                                                      >  1:520227/1‑100 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                 >  1:1309042/1‑100 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                 >  2:3340960/1‑100 (MQ=255)
                                            gAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                <  1:1966418/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                              >  1:85385/1‑100 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                            >  1:807851/1‑100 (MQ=255)
                                                 taaaaTCCTCTCGTGCCTACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                              <  2:419631/95‑1 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAGCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                          >  2:1949293/1‑100 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAACCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                          >  2:1949276/1‑100 (MQ=255)
                                                   aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                         >  2:3151236/1‑100 (MQ=255)
                                                     aTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcac                                       <  2:2059015/100‑1 (MQ=255)
                                                        ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                    <  1:2960263/100‑1 (MQ=255)
                                                        ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                    >  1:1640456/1‑100 (MQ=255)
                                                         tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                   <  1:2527778/100‑1 (MQ=255)
                                                           tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                 <  2:85385/100‑1 (MQ=255)
                                                           tcGTGCCGAACAAAAATCACAAGAAAAAACCAACCCTGACGGTTGGTTTTTTTATATCGGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                 <  2:3124059/100‑1 (MQ=255)
                                                             gTGCCGACCAACAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                               <  1:1699664/100‑1 (MQ=255)
                                                             gTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                               >  1:2005328/1‑100 (MQ=255)
                                                              tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                     <  2:942202/93‑1 (MQ=255)
                                                                    ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                        <  1:3100361/100‑1 (MQ=255)
                                                                     cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                       <  2:1370817/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                          <  2:2724468/97‑1 (MQ=255)
                                                                       ataaTCCCAAGAAAAAAACAACCCTTACTGTTAGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATTACAGCCTGGTTAGTGTCTGt                     <  2:43757/98‑1 (MQ=255)
                                                                       aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                     <  1:65891/100‑1 (MQ=255)
                                                                        aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                    <  1:916709/100‑1 (MQ=255)
                                                                          aTCCCAAGAAAAAACCAACCCTTACTGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAt                  <  2:93675/100‑1 (MQ=255)
                                                                           tCCCAAGAAAAAACCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                 <  1:851738/100‑1 (MQ=38)
                                                                           tCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                 <  1:1374421/100‑1 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                <  1:3098130/100‑1 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                <  1:3238/100‑1 (MQ=255)
                                                                             ccAAGAAACAACCAACCCTTACGGCTGGTTTTTTTATATCTGCAATTAATTCGATATACAGACCGTGACTCATCATAGCATGTTTATTTTCTGTTATCAg               >  2:1657767/1‑100 (MQ=255)
                                                                             ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg               >  1:118651/1‑100 (MQ=255)
                                                                             ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg               >  2:1409463/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              >  2:1889431/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              >  1:1394344/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa              <  1:1926471/100‑1 (MQ=255)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa             >  2:2815370/4‑100 (MQ=38)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa             >  2:346779/4‑100 (MQ=38)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAGACAGACCGTTACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  2:301490/3‑100 (MQ=38)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  1:1890527/3‑100 (MQ=38)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  1:593702/3‑100 (MQ=38)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg           <  2:1169231/100‑1 (MQ=255)
                                                                                 aaaaaaaCCAACCGTCAGGGGTGGTTTTTTTATATCTGCAATCAATACGATAAACAGACCGTGACCCCTCCCAGCCTGTTTATTTTCTGTAAt                  >  2:740945/2‑93 (MQ=25)
                                                                                 aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg           >  1:320910/2‑100 (MQ=38)
                                                                                  aaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc            >  2:1457232/1‑98 (MQ=255)
                                                                                   aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAAATCGATAACCAGACCGTAACACAACAAAGCAAGTTTATTTACTGTTAAAAGAACGt          >  1:2121930/1‑99 (MQ=255)
                                                                                    aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc        >  1:2700724/1‑100 (MQ=255)
                                                                                      aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg      >  1:1443760/1‑100 (MQ=255)
                                                                                       tctaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGa     <  1:1199012/97‑1 (MQ=255)
                                                                                         caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc   <  1:3209199/100‑1 (MQ=255)
                                                                                          aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGaaga  >  2:1199004/1‑97 (MQ=255)
                                                                                           accTTTACGGTTGGTTTTTTTTTATCTGCAATTAATTAGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca  <  2:2280387/99‑1 (MQ=255)
                                                                                           accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc   <  2:1535688/98‑1 (MQ=255)
                                                                                                  |                                                                                            
GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCA  >  NC_000913/2286217‑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: