Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,505 G→A 6.2% intergenic (+484/‑40) yfjW → / → ypjI CP4‑57 prophage; uncharacterized protein YfjW/putative UPF0401 protein YpjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,775,5050GA 6.2% 302.9 / 6.5 96intergenic (+484/‑40)yfjW/ypjICP4‑57 prophage; uncharacterized protein YfjW/putative UPF0401 protein YpjI
Reads supporting (aligned to +/- strand):  ref base G (44/46);  new base A (3/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGA  >  NC_000913/2775356‑2775654
                                                                                                                                                     |                                                                                                                                                     
gTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAg                                                                                                                                                       <  1:134499/150‑1 (MQ=255)
 tAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAgg                                                                                                                                                      >  2:1497859/1‑150 (MQ=255)
   aaTGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGAt                                                                                                                                                    >  2:241566/1‑150 (MQ=255)
   aaTGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGAt                                                                                                                                                    <  1:1394180/150‑1 (MQ=255)
        cccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGc                                                                                                                                               <  2:40962/150‑1 (MQ=255)
          cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCg                                                                                                                                             >  2:497083/1‑150 (MQ=255)
          cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCg                                                                                                                                             <  2:37377/150‑1 (MQ=255)
          cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCg                                                                                                                                             <  2:1065412/150‑1 (MQ=255)
           tACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCgt                                                                                                                                            <  1:484868/150‑2 (MQ=255)
                aaTGAACCATAATATTTATAGGTCGTCTAGCACAACAACTTCATCTATTGTAACACAACGGCCTATCAGAAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCttttt                                                                                                                                       >  2:496744/1‑150 (MQ=255)
                 aTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttttt                                                                                                                                      <  1:121128/150‑1 (MQ=255)
                  tGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTtta                                                                                                                                     <  2:1308304/150‑1 (MQ=255)
                    aaCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatt                                                                                                                                   <  2:1683389/150‑1 (MQ=255)
                     aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                                  <  1:644430/150‑1 (MQ=255)
                       cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                                >  2:276995/1‑150 (MQ=255)
                        attattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttt                                                                                                                               >  1:678954/1‑150 (MQ=255)
                         ttattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAtttt                                                                                                                              >  1:514898/1‑150 (MQ=255)
                         ttattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAtttt                                                                                                                              >  1:1342879/1‑150 (MQ=255)
                         ttattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAtttt                                                                                                                              <  2:829007/150‑1 (MQ=255)
                            tttttAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCa                                                                                                                           >  1:1665283/1‑150 (MQ=255)
                             ttttAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAg                                                                                                                          >  2:1370532/1‑150 (MQ=255)
                                 agagGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGa                                                                                                                      >  2:1659185/1‑150 (MQ=255)
                                  gagGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGag                                                                                                                     >  2:869699/1‑150 (MQ=255)
                                     ggCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGagtag                                                                                                                  <  1:667805/150‑1 (MQ=255)
                                     ggCGTCTGGCACTCCAACTACCGCATTTCGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATAACATTTCTACAAACATACCCGATTAACGGCACAATCTCAGCATTCTGCCGCTATTTTAATATATTTCAGGGGagtag                                                                                                                  >  2:1524063/1‑150 (MQ=255)
                                          cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                             <  1:1370532/150‑1 (MQ=255)
                                            ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAACAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                           <  1:1301756/150‑1 (MQ=255)
                                            ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCGGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                           >  2:544172/1‑150 (MQ=255)
                                            ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                           >  1:1811039/1‑150 (MQ=255)
                                            ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                           >  2:639332/1‑150 (MQ=255)
                                            ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                           >  2:479061/1‑150 (MQ=255)
                                             gCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGc                                                                                                          <  1:1373715/150‑1 (MQ=255)
                                              cACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCa                                                                                                         <  1:29067/150‑1 (MQ=255)
                                                cTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAAc                                                                                                       <  1:241566/150‑1 (MQ=255)
                                                cTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAAc                                                                                                       <  1:568067/150‑1 (MQ=255)
                                                cTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAAc                                                                                                       >  2:888881/1‑150 (MQ=255)
                                                    aaCTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTg                                                                                                   <  1:276995/150‑1 (MQ=255)
                                                      cTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGaa                                                                                                 >  1:1736865/1‑150 (MQ=255)
                                                      attCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGaa                                                                                                 <  1:1951348/149‑1 (MQ=255)
                                                       ttCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAgggg                                                                                                                       <  2:1362905/127‑1 (MQ=255)
                                                       ttCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAgggg                                                                                                                       >  1:1362905/1‑127 (MQ=255)
                                                       ttCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAg                                                                                                <  2:368821/150‑1 (MQ=255)
                                                       ttCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAg                                                                                                <  1:950430/150‑1 (MQ=255)
                                                            cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgt                                                                                           >  1:1687417/1‑149 (MQ=255)
                                                              tttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTa                                                                                         <  2:868882/150‑1 (MQ=255)
                                                               ttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAg                                                                                        <  2:900646/150‑1 (MQ=255)
                                                                tGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGt                                                                                       <  1:869699/150‑1 (MQ=255)
                                                                  gCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCa                                                                                     >  2:1424068/1‑150 (MQ=255)
                                                                       taaCGGCCTATCATCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTa                                                                                <  2:114334/149‑1 (MQ=255)
                                                                         aCGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGc                                                                              <  2:1736865/150‑1 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCAACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                            <  2:1526164/150‑1 (MQ=255)
                                                                              cTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACg                                                                         >  2:1586125/1‑150 (MQ=255)
                                                                               tATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGc                                                                        >  2:145519/1‑150 (MQ=255)
                                                                                   aGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCTGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTc                                                                    <  2:704857/150‑1 (MQ=255)
                                                                                   aGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCtttg                                                                    >  2:1991491/1‑149 (MQ=255)
                                                                                       atatAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGt                                                                <  1:1497859/150‑1 (MQ=255)
                                                                                         ataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCt                                                              >  2:1060925/1‑150 (MQ=255)
                                                                                         ataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCt                                                              >  1:1920207/1‑150 (MQ=255)
                                                                                          taAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTg                                                             >  2:761693/1‑150 (MQ=255)
                                                                                          taAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTg                                                             >  1:1370929/1‑150 (MQ=255)
                                                                                          taAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTg                                                             >  2:1648814/1‑150 (MQ=255)
                                                                                           aaaaGCATTCAGATATATATCATTTCTACAAACATACCCGATAAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGc                                                            <  2:1665283/150‑1 (MQ=255)
                                                                                            aaaGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGcc                                                           >  2:230959/1‑150 (MQ=255)
                                                                                               gCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCGGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCa                                                        <  1:544172/150‑1 (MQ=255)
                                                                                                  ttCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTg                                                     >  1:658010/1‑150 (MQ=255)
                                                                                                   tCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTgg                                                    <  1:761693/150‑1 (MQ=255)
                                                                                                       atatatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCAt                                                <  1:853865/150‑1 (MQ=255)
                                                                                                          tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                             >  2:1251801/1‑150 (MQ=255)
                                                                                                             atCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTg                                          >  2:1868394/1‑150 (MQ=255)
                                                                                                               cATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAt                                        >  1:1100369/1‑150 (MQ=255)
                                                                                                                aTTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAtt                                       >  2:1103922/1‑150 (MQ=255)
                                                                                                                 tttCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTg                                      <  1:1348636/150‑1 (MQ=255)
                                                                                                                   tCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                    <  2:1920207/150‑1 (MQ=255)
                                                                                                                    cTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGt                                   <  1:1668005/150‑1 (MQ=255)
                                                                                                                      aCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAg                                 <  2:1100369/150‑1 (MQ=255)
                                                                                                                       cAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                                >  2:618236/1‑150 (MQ=255)
                                                                                                                        aaaCATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTACGTAGCa                               <  1:1081474/150‑1 (MQ=255)
                                                                                                                           cATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAggg                            <  1:479061/150‑1 (MQ=255)
                                                                                                                           cATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAggg                            <  1:1622606/150‑1 (MQ=255)
                                                                                                                              aCCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCa                         >  2:1166152/1‑150 (MQ=255)
                                                                                                                               cccGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCAc                        >  2:1379865/1‑150 (MQ=255)
                                                                                                                                 cGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTg                      <  2:392986/150‑1 (MQ=255)
                                                                                                                                    ttAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTa                   <  1:496744/150‑1 (MQ=255)
                                                                                                                                      aaCGGCAGCATCTAAGGATTCTGCCGCATTTTTATTATTTTCCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCCCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATc                 >  2:1786807/1‑150 (MQ=255)
                                                                                                                                      aaCGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATc                 <  1:230959/150‑1 (MQ=255)
                                                                                                                                      aaCGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTCCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATc                 >  2:1786802/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg               <  1:639332/150‑1 (MQ=255)
                                                                                                                                        cGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg               >  1:368774/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTgg                     <  2:538328/144‑1 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTgg                     >  1:538328/1‑144 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg               >  2:635968/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg               >  1:477984/1‑150 (MQ=255)
                                                                                                                                         ggCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGATTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc              >  2:1962749/1‑150 (MQ=255)
                                                                                                                                          gCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGcc             <  1:1586125/150‑1 (MQ=255)
                                                                                                                                          gCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGcc             <  1:1131734/150‑1 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg           >  1:1699151/1‑150 (MQ=255)
                                                                                                                                               aTCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTAt        >  1:925942/1‑150 (MQ=255)
                                                                                                                                                    aGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCg   <  2:1370929/150‑1 (MQ=255)
                                                                                                                                                     ggATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCTTAGCAGGGCCACTGGTATCTGCCGGTATCACCGa  <  2:678954/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGA  >  NC_000913/2775356‑2775654

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