Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,068,408 T→A 15.9% intergenic (‑1087/‑227) insH‑6 ← / → yoeA CP4‑44 prophage; IS5 transposase and trans‑activator/CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,068,4080TA15.9% 200.1 / 25.4 81intergenic (‑1087/‑227)insH‑6/yoeACP4‑44 prophage; IS5 transposase and trans‑activator/CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base A (6/7);  total (37/44)
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.99e-01

ACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAT  >  NC_000913/2068260‑2068557
                                                                                                                                                    |                                                                                                                                                     
aCCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATg                                                                                                                                                      <  1:850661/150‑1 (MQ=255)
   tgatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCa                                                                                                                                                   <  1:103537/150‑1 (MQ=255)
      tactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAaga                                                                                                                                                <  2:863053/150‑1 (MQ=255)
        ctgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGt                                                                                                                                              >  1:382100/1‑150 (MQ=255)
         tgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTg                                                                                                                                             >  1:688344/1‑150 (MQ=255)
          tatatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                            <  1:339633/149‑1 (MQ=255)
          gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                            <  1:1049707/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                          <  1:258118/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                          <  1:1253236/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  1:321034/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  1:699782/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  1:1252017/150‑1 (MQ=255)
            tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  2:1121856/150‑1 (MQ=255)
              tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTc                                                                                                                                        <  2:946647/150‑1 (MQ=255)
              tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTc                                                                                                                                        <  1:1207641/150‑1 (MQ=255)
              tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                        <  1:737927/150‑1 (MQ=255)
               gAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTcc                                                                                                                                       <  1:490357/150‑1 (MQ=39)
                  tCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTg                                                                                                                                    <  2:853550/150‑1 (MQ=255)
                   cTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCGCTCTTATGTCAAGAGTGCGGTCCGTgg                                                                                                                                   >  1:1207794/1‑150 (MQ=255)
                    ttACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGTGGTCCGTGGa                                                                                                                                  >  2:768787/1‑150 (MQ=255)
                      aCCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATa                                                                                                                                >  1:351703/1‑150 (MQ=255)
                        cGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  2:71857/138‑1 (MQ=255)
                        cGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          >  1:71857/1‑138 (MQ=255)
                        cGCATATACGGTCAATTCCAGCAATTAATTACATTATCCACGCCAAAGTATTTATCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          >  1:144606/1‑138 (MQ=255)
                        cGCATATACGGTCAATTCCAGCAATTAATTACATTATCCACGCCAAAGTATTTATCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                          <  2:144606/138‑1 (MQ=255)
                         gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGaaaacg                                                                                                                             <  1:547502/150‑6 (MQ=38)
                           atatACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAAcc                                                                                                                           <  2:495433/150‑1 (MQ=255)
                           atatACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAAcc                                                                                                                           <  1:1350092/150‑1 (MQ=255)
                             ataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCag                                                                                                                         <  1:1197771/150‑1 (MQ=255)
                              taCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGACTGCGGTCCGTGGATACAACCaga                                                                                                                        >  2:624693/1‑150 (MQ=255)
                               aCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCagag                                                                                                                       >  2:1419951/1‑150 (MQ=255)
                                 ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCTGTCCGTGGATACAACCAGAGAc                                                                                                                     <  2:1097439/150‑1 (MQ=255)
                                  gTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGAcc                                                                                                                    >  1:382310/1‑150 (MQ=255)
                                     aaTTTCAGCAATTAATTACTTTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGAc                                                                                                                 <  2:208880/150‑1 (MQ=255)
                                            gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                          <  2:1431498/150‑1 (MQ=255)
                                            gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                          <  1:1104768/150‑1 (MQ=255)
                                                  aattaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCag                                                                                                    >  2:1008296/1‑150 (MQ=255)
                                                     taCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAgg                                                                                                 <  2:754692/150‑1 (MQ=255)
                                                        aTTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                   <  2:1277393/145‑1 (MQ=255)
                                                        aTTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                   >  1:1277393/1‑145 (MQ=255)
                                                        aTTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGaaa                                                                                              >  1:550364/1‑150 (MQ=255)
                                                          aaTCCACCCCAAACTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCATCTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAacg                                                                                            <  2:243802/149‑1 (MQ=255)
                                                              cACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACgg                                                                                        <  2:1120164/150‑1 (MQ=255)
                                                                 gCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGata                                                                                     >  1:1063396/1‑150 (MQ=255)
                                                                  ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACAGCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGatat                                                                                    >  1:390445/1‑150 (MQ=255)
                                                                  ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGatat                                                                                    >  1:7643/1‑150 (MQ=255)
                                                                         aTTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                             >  1:241321/1‑150 (MQ=255)
                                                                         aTTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                             <  1:424924/150‑1 (MQ=255)
                                                                                  cacaATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAata                                                                    >  2:1105924/1‑150 (MQ=255)
                                                                                  cacaATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAata                                                                    >  1:18248/1‑150 (MQ=255)
                                                                                    caATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAatata                                                                  >  2:770334/1‑149 (MQ=255)
                                                                                               ccTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAAc                                                       <  1:1414051/150‑1 (MQ=255)
                                                                                                   cttTCAGACACCATTTTTTCAACTCCGTTTACCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGa                                                   <  1:1349197/150‑1 (MQ=255)
                                                                                                    tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGAt                                                  <  2:166305/150‑1 (MQ=255)
                                                                                                    tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGAt                                                  <  2:161226/150‑1 (MQ=255)
                                                                                                         gACACCATTTTTTCAACTCCGTTTACCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                             >  2:820561/1‑150 (MQ=255)
                                                                                                           cacCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTc                                           >  1:1009262/1‑150 (MQ=255)
                                                                                                             ccATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaa                                         >  2:421950/1‑150 (MQ=255)
                                                                                                              cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                        >  1:647588/1‑150 (MQ=255)
                                                                                                                     tCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGttt                                 <  1:624693/150‑1 (MQ=255)
                                                                                                                               tttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACa                       >  1:1249716/1‑150 (MQ=255)
                                                                                                                               tttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACa                       <  2:93395/150‑1 (MQ=255)
                                                                                                                                  ccACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGCGACCGACTGACACGAGTCAGCGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTTAGGTCAAAAATGTTTGATATTTACAAtt                    >  1:1292409/1‑150 (MQ=38)
                                                                                                                                  ccACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                                >  2:732121/1‑138 (MQ=38)
                                                                                                                                  ccACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                                <  1:732121/138‑1 (MQ=38)
                                                                                                                                   cACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAttt                   >  1:244461/1‑150 (MQ=39)
                                                                                                                                   cACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAttt                   <  2:1063396/150‑1 (MQ=255)
                                                                                                                                   cACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAttt                   >  1:1251458/1‑150 (MQ=255)
                                                                                                                                    aCGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTa                  >  2:858472/1‑150 (MQ=39)
                                                                                                                                     cGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt                 >  2:1048198/1‑150 (MQ=39)
                                                                                                                                      ggACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTTAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATg                <  2:382100/150‑1 (MQ=255)
                                                                                                                                      ggACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATg                >  1:748991/1‑150 (MQ=255)
                                                                                                                                       gACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa               >  2:353266/1‑150 (MQ=255)
                                                                                                                                         ccGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAg             <  1:1419951/150‑1 (MQ=255)
                                                                                                                                         ccGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAg             >  1:1306143/1‑150 (MQ=255)
                                                                                                                                          cGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGa            >  1:1147675/1‑150 (MQ=255)
                                                                                                                                            cACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATg          <  1:1105924/150‑1 (MQ=255)
                                                                                                                                               tctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACa       <  2:390445/150‑1 (MQ=255)
                                                                                                                                                ctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaa      >  1:111009/1‑150 (MQ=255)
                                                                                                                                                 ttATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaat     >  1:486440/1‑150 (MQ=255)
                                                                                                                                                 ttATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTCTGATATTTACAATTTATGAAGATGACaat     <  1:770334/150‑1 (MQ=255)
                                                                                                                                                    tGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaataat  >  1:158323/1‑150 (MQ=255)
                                                                                                                                                    tGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaataat  >  1:982338/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAT  >  NC_000913/2068260‑2068557

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: