Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,273,229 A→T 7.6% intergenic (+27/+17) chaC → / ← ychN glutathione‑specific gamma‑glutamylcyclotransferase/DsrE/F sulfur relay family protein YchN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,273,2290AT 7.6% 243.6 / 7.8 79intergenic (+27/+17)chaC/ychNglutathione‑specific gamma‑glutamylcyclotransferase/DsrE/F sulfur relay family protein YchN
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base T (4/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACC  >  NC_000913/1273082‑1273378
                                                                                                                                                   |                                                                                                                                                     
ggAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCCGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGa                                                                                                                                                     >  1:1221513/1‑150 (MQ=255)
     aGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGccc                                                                                                                                                <  2:46760/150‑1 (MQ=255)
        aGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccac                                                                                                                                             <  1:127674/150‑3 (MQ=255)
           ttATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactat                                                                                                                                          <  1:413294/150‑6 (MQ=255)
           ttATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCaa                                                                                                                                          <  1:1397491/150‑1 (MQ=255)
             aTCAAACTGGTAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCTGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAAc                                                                                                                                        <  1:1206146/150‑1 (MQ=255)
             aTCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAAc                                                                                                                                        >  1:326796/1‑150 (MQ=255)
             aTCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAAc                                                                                                                                        <  1:1414414/150‑1 (MQ=255)
              tCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACa                                                                                                                                       <  2:335355/150‑1 (MQ=255)
                aaaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACAtt                                                                                                                                     >  1:1409943/1‑150 (MQ=255)
                 aaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTa                                                                                                                                    >  2:950745/1‑150 (MQ=255)
                     gggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTTCGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAt                                                                                                                                <  2:850005/150‑1 (MQ=255)
                     gggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAt                                                                                                                                >  2:608982/1‑150 (MQ=255)
                        aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTTGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgtg                                                                                                                             <  2:1196036/150‑1 (MQ=255)
                        aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgtg                                                                                                                             <  1:1144869/150‑1 (MQ=255)
                          tGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAg                                                                                                                           >  2:1526764/1‑150 (MQ=255)
                           gCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTTGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGc                                                                                                                          <  2:994939/150‑1 (MQ=255)
                            cAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCa                                                                                                                         >  1:830401/1‑150 (MQ=255)
                              ggACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGCCGCCCCCTCAAACATTAAAATGTGAGCACt                                                                                                                       >  1:904735/1‑150 (MQ=255)
                                     gggCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCgg                                                                                                                <  2:724834/150‑1 (MQ=255)
                                      ggCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAATAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                                                                               >  2:1325772/1‑150 (MQ=255)
                                       gCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                              <  2:1221513/150‑1 (MQ=255)
                                        cTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTCTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                                                                             <  1:993942/150‑1 (MQ=255)
                                          gAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGAc                                                                                                           >  2:283320/1‑150 (MQ=255)
                                           aaTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACa                                                                                                          <  1:206021/150‑1 (MQ=255)
                                             tGAGTTGCTGGTATCGGTAAAAAAACTGCAGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGCCGCCCCCTCAAACATTAAAAGGTGAGGACTTTATCGGCTGAccgg                                                                                                        >  2:681639/1‑147 (MQ=255)
                                                        tATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcgc                                                                                             >  1:1143755/1‑150 (MQ=255)
                                                            ggTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAACCATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGt                                                                                         >  2:528308/1‑150 (MQ=255)
                                                            ggTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGt                                                                                         >  2:644626/1‑150 (MQ=255)
                                                                 aaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCcac                                                                                    <  1:1086579/150‑1 (MQ=255)
                                                                  aaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCcaca                                                                                   <  1:1518602/150‑1 (MQ=255)
                                                                      ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGCCGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAagag                                                                               >  2:677771/1‑150 (MQ=255)
                                                                      ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAagag                                                                               <  2:1409943/150‑1 (MQ=255)
                                                                        gctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTa                                                                             >  1:1495768/1‑150 (MQ=255)
                                                                          tgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTAcc                                                                           >  2:880944/1‑150 (MQ=255)
                                                                          tgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTAcc                                                                           >  1:872103/1‑150 (MQ=255)
                                                                                 gaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTcc                                                                    <  2:485174/150‑1 (MQ=255)
                                                                                 gaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCAAGTTCCACAAGAGTACCGATTTcc                                                                    <  2:830401/150‑1 (MQ=255)
                                                                                  aaTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCa                                                                   <  2:1143755/150‑1 (MQ=255)
                                                                                    ttttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAATCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                 >  1:597538/1‑150 (MQ=255)
                                                                                     tttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAccc                                                                <  1:752750/150‑1 (MQ=255)
                                                                                         cGGATGGTGTGTTACGTCCGGGATTCGCCTTAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatc                                                            <  2:71400/150‑1 (MQ=255)
                                                                                          ggATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatca                                                           <  1:950745/150‑1 (MQ=255)
                                                                                              ggTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaatca                                                       <  1:644626/150‑1 (MQ=255)
                                                                                               gtgtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCag                                                      <  2:1462675/150‑1 (MQ=255)
                                                                                                tgtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCaga                                                     <  1:829365/150‑1 (MQ=255)
                                                                                                 gtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCagag                                                    >  1:1066152/1‑150 (MQ=255)
                                                                                                 gtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCagag                                                    <  1:1606938/150‑1 (MQ=255)
                                                                                                  tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                 >  2:744529/1‑136 (MQ=255)
                                                                                                  tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                 <  1:744529/136‑1 (MQ=255)
                                                                                                  tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAgg                                                   <  2:1495768/150‑1 (MQ=255)
                                                                                                  tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAgg                                                   >  1:1190795/1‑150 (MQ=255)
                                                                                                  tgtTACATCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                 >  2:1273249/1‑136 (MQ=255)
                                                                                                  tgtTACATCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                 <  1:1273249/136‑1 (MQ=255)
                                                                                                     tACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAg                                                >  2:265207/1‑150 (MQ=255)
                                                                                                     tACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAg                                                >  2:971009/1‑150 (MQ=255)
                                                                                                      aCGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgt                                               <  2:1429781/150‑1 (MQ=255)
                                                                                                         tacGGCATTCGCCTGAGTAAACATCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACAATAAAATGTGAGCACTTTATCGGCTGACAGCGTCCAGTGCGCGAGTTCCACAAGACTACCGATTTCGCCCCCATCAATCAGAGGAAGTGTa                                            >  2:1066088/3‑150 (MQ=255)
                                                                                                         tCCGTGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTa                                            <  1:1052010/150‑1 (MQ=255)
                                                                                                            gggATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGAGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTa                                         >  2:1340393/1‑150 (MQ=255)
                                                                                                                ttCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATcc                                     <  2:1066152/150‑1 (MQ=255)
                                                                                                                  cGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcg                                   >  2:49326/1‑150 (MQ=255)
                                                                                                                  cGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcg                                   >  1:293610/1‑150 (MQ=255)
                                                                                                                   gCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcgc                                  >  2:918182/1‑150 (MQ=255)
                                                                                                                       gAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCg                                                                                                                 <  1:607199/67‑1 (MQ=255)
                                                                                                                       gAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCg                                                                                                                 >  2:607199/1‑67 (MQ=255)
                                                                                                                         gTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGt                            >  2:1072308/1‑150 (MQ=255)
                                                                                                                            aaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCgg                         <  2:796793/150‑1 (MQ=255)
                                                                                                                              cTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTa                       >  2:1480360/1‑150 (MQ=255)
                                                                                                                                  ccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAgg                   <  1:265207/150‑1 (MQ=255)
                                                                                                                                  ccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAgg                   >  2:120227/1‑150 (MQ=255)
                                                                                                                                     tttgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGttt                >  2:397277/7‑150 (MQ=255)
                                                                                                                                       tgagggggCTTCTGACGCCCCCTCAAACATTAAAATGTGAGCACATTATCGTCTGAAAGCCACCATTGCGACAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGATGAAGTGTACTAATACAGAGCCCGTCCGTACAGCttttt              >  2:86241/5‑149 (MQ=255)
                                                                                                                                       tgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTg              >  1:1380287/5‑150 (MQ=255)
                                                                                                                                       tAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTg              >  1:1284634/1‑150 (MQ=255)
                                                                                                                                        gagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGc             >  2:1507350/4‑150 (MQ=255)
                                                                                                                                           ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGcaca          >  1:739571/1‑150 (MQ=255)
                                                                                                                                                  cAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAc   <  1:1325772/150‑1 (MQ=255)
                                                                                                                                                  cAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAc   >  2:1271107/1‑150 (MQ=255)
                                                                                                                                                   aGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAcc  <  1:880944/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACC  >  NC_000913/1273082‑1273378

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: