Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,504 A→G 8.2% intergenic (+483/‑41) 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,5040AG 8.2% 284.3 / 11.5 97intergenic (+483/‑41)yfjW/ypjICP4‑57 prophage; uncharacterized protein YfjW/putative UPF0401 protein YpjI
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base G (4/4);  total (48/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 = 9.71e-01

GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCG  >  NC_000913/2775356‑2775653
                                                                                                                                                    |                                                                                                                                                     
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:368774/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg              <  1:639332/150‑1 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTgg                    <  2:538328/144‑1 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTgg                    >  1:538328/1‑144 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg              >  1:477984/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg              >  2:635968/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)
                                                                                                                                                    |                                                                                                                                                     
GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCG  >  NC_000913/2775356‑2775653

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: