Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,291 G→A 5.3% intergenic (+174/‑1) pagP → / → cspE phospholipid:lipid A palmitoyltransferase/constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,2910GA 5.3% 253.5 / 2.8 75intergenic (+174/‑1)pagP/cspEphospholipid:lipid A palmitoyltransferase/constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (25/46);  new base A (2/2);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGC  >  NC_000913/657145‑657438
                                                                                                                                                  |                                                                                                                                                   
aaaaTCAATATCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTTCTGGTAACCGACACAGCTTTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  4:118345/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGCCACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  5:66736/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  4:140946/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  6:167156/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  6:86305/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  1:302863/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  2:185825/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  7:242165/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  8:214802/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  1:257665/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  8:302352/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  1:219572/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  3:277688/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  8:318126/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  3:238413/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  3:128238/147‑1 (MQ=255)
aaaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATCGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  1:7060/147‑1 (MQ=255)
 aaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCTTGTAAAGGTAATTTTGa                                                                                                                                                    <  5:238455/147‑1 (MQ=255)
 aaaTCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGa                                                                                                                                                    <  1:11460/147‑1 (MQ=255)
       aTAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTaa                                                                                                                                             <  3:170558/148‑1 (MQ=255)
        tAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAg                                                                                                                                            <  8:23744/148‑1 (MQ=255)
           cTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGAtt                                                                                                                                         <  7:170277/148‑1 (MQ=255)
           cTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGAtt                                                                                                                                         <  1:13848/148‑1 (MQ=255)
            tAAAATTAGTTACATCATAATTAACGCGTCTTTTATCACTTTTTAGTAAAGTTTGACTGGATAAAGCGTTGCAAAATTGTTGTACTGGTAACCGACTCAGGATTTGTGTCTATTTTTCATGTAAAGGTTATTTTg                                                                                                                                                     <  1:20995/135‑1 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  2:120096/134‑1 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     >  1:120096/1‑134 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAg                                                                                                                                            >  1:36295/1‑143 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAg                                                                                                                                            <  2:36295/143‑1 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTa                                                                                                                                        >  5:14359/1‑147 (MQ=255)
             aaaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTa                                                                                                                                        >  1:19490/1‑147 (MQ=255)
              aactAAGTAACATCAATAATAACGCGACTTTTATCTCTTTTTAGTAAAGTTACACTGGACAATGCGTACCACAATTGGTGTACTGTTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTaa                                                                                                                                       <  2:21803/144‑1 (MQ=255)
              aaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTaa                                                                                                                                       <  4:214720/147‑1 (MQ=255)
              aaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTaa                                                                                                                                       <  4:140086/147‑1 (MQ=255)
              aaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTaa                                                                                                                                       <  1:166416/147‑1 (MQ=255)
              aaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTaa                                                                                                                                       <  2:207397/147‑1 (MQ=255)
              aaaTAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTa                                                                                                                                        >  1:254329/1‑146 (MQ=255)
                   aGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTa                                                                                                                                  <  2:261692/147‑1 (MQ=255)
                   aGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTa                                                                                                                                  <  2:158230/147‑1 (MQ=255)
                        cATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTa                                                                                                                            <  7:5731/148‑1 (MQ=255)
                        cATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTa                                                                                                                            <  2:115967/148‑1 (MQ=255)
                          tCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAg                                                                                                                          >  6:125109/1‑148 (MQ=255)
                            aaaaaTAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGt                                                                                                                         <  1:67542/147‑1 (MQ=255)
                              aaaTAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGt                                                                                                                      <  6:1550/148‑1 (MQ=255)
                                 tAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   <  2:254329/148‑1 (MQ=255)
                                 tAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   <  7:325112/148‑1 (MQ=255)
                                 tAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   <  2:19490/148‑1 (MQ=255)
                                 tAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   <  1:28829/148‑1 (MQ=255)
                                   aCGCGACTTTTATCACTTTTTAGTTAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGGATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTaa                                                                                                                  <  8:212784/147‑1 (MQ=255)
                                   aCGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTaa                                                                                                                  <  7:46390/147‑1 (MQ=255)
                                    cgcgACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATg                                                                                                                <  8:124177/148‑1 (MQ=255)
                                         cTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAAc                                                                                                                                >  8:153869/1‑127 (MQ=255)
                                         cTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAAc                                                                                                                                <  7:153869/127‑1 (MQ=255)
                                           tttATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTAATGTCTAAGATTAAAGGTAACTTTAAGTGGTTTAATGAGTCCa                                                                                                          >  3:244087/1‑147 (MQ=255)
                                           tttATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTAATGTCTAAGATTAAAGGTAACTTTAAGTGGTTTAATGAGTCCa                                                                                                          >  5:135963/1‑147 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGCCACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCgg                                                                                                >  6:88895/1‑148 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     <  4:2496/95‑1 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTg                                                                                                                                                     >  3:2496/1‑95 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTa                                                                                                                                  <  1:220579/114‑1 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTa                                                                                                                                  >  2:220579/1‑114 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCgg                                                                                                >  5:33882/1‑148 (MQ=255)
                                                    ttttAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCgg                                                                                                >  7:120319/1‑148 (MQ=255)
                                                        aGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGttt                                                                                             <  6:33882/147‑1 (MQ=255)
                                                         gTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCa                                                                                           <  6:179939/148‑1 (MQ=255)
                                                            aaGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCAtt                                                                                         <  4:224696/147‑1 (MQ=255)
                                                              gTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   >  7:41421/1‑119 (MQ=255)
                                                              gTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTa                                                                                                                   <  8:41421/119‑1 (MQ=255)
                                                                tACAATGGAAAAAGCGTACCACAATTTGTGTAATGTTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTc                                                                                     <  5:48193/147‑1 (MQ=255)
                                                                  cacTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCgg                                                                                  <  8:308099/148‑1 (MQ=255)
                                                                   acTGGACAAAGCGTACCACAATTTGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTATCATGTAAATGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGa                                                                                 <  4:12689/148‑1 (MQ=255)
                                                                   acTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGa                                                                                 <  3:16109/148‑1 (MQ=255)
                                                                          aaaGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACgg                                                                           >  8:48724/1‑147 (MQ=255)
                                                                           aaGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACgg                                                                           >  3:111123/1‑146 (MQ=255)
                                                                                  ccACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTAATGTCTAAGATTAAAGGTAACTTTAAGTGGTTTAATGAGTCCAAATGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAg                                                                   <  6:135963/147‑1 (MQ=255)
                                                                                  ccACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTAATGTCTAAGATTAAAGGTAACTTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAg                                                                   <  4:244087/147‑1 (MQ=255)
                                                                                   cacaATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAc                                                                 >  1:299235/1‑148 (MQ=255)
                                                                                    acaATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAg                                                                               >  6:54175/1‑133 (MQ=255)
                                                                                    acaATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAg                                                                               <  5:54175/133‑1 (MQ=255)
                                                                                        ttGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGtt                                                            <  4:111123/148‑1 (MQ=255)
                                                                                           gtgtACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                          <  7:48724/147‑1 (MQ=255)
                                                                                           gtgtACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGt                                                         <  4:29588/148‑1 (MQ=255)
                                                                                                  ggTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACg                                                                            <  7:176804/122‑1 (MQ=255)
                                                                                                  ggTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACg                                                                            >  8:176804/1‑122 (MQ=255)
                                                                                                   gTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtc                                                  <  5:125109/147‑1 (MQ=255)
                                                                                                    tAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                >  1:280405/1‑148 (MQ=255)
                                                                                                                   ttGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                       <  3:155083/126‑1 (MQ=255)
                                                                                                                   ttGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                       >  4:155083/1‑126 (MQ=255)
                                                                                                                      tgtCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATgg                              >  1:88436/1‑148 (MQ=255)
                                                                                                                        tCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                             <  2:299235/147‑1 (MQ=255)
                                                                                                                                   tGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCtt                 >  2:190424/1‑148 (MQ=255)
                                                                                                                                      aaaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt              >  2:67840/1‑148 (MQ=255)
                                                                                                                                       aaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt              >  4:4943/1‑147 (MQ=255)
                                                                                                                                        aGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg             >  5:174985/1‑147 (MQ=255)
                                                                                                                                        aGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa            <  8:120319/148‑1 (MQ=255)
                                                                                                                                        aGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa            <  5:88895/148‑1 (MQ=255)
                                                                                                                                           tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAgg         >  4:400/1‑148 (MQ=255)
                                                                                                                                           tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAgg         >  3:32151/1‑148 (MQ=255)
                                                                                                                                             aTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc       <  6:174985/148‑1 (MQ=255)
                                                                                                                                                  gATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcgc  >  3:262764/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGC  >  NC_000913/657145‑657438

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