Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,473,8820TG20.9% 129.6 / 10.4 67L94R (CTG→CGG) yfdPCPS‑53 (KpLE1) prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base G (14/0);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGA  >  NC_000913/2473741‑2474023
                                                                                                                                             |                                                                                                                                              
tcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTg                                                                                                                                               >  2:131369/1‑143 (MQ=255)
tcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTg                                                                                                                                               >  1:763629/1‑143 (MQ=255)
 caACCTGGACATCGAACCATTTGCTGCGAACCGTCTGAAAATCAAACACAAACTCCCATATGTTCGCCCGCCTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGa                                                                                                                                              <  2:384250/143‑1 (MQ=255)
     cTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACTAt                                                                                                                                          >  1:690993/1‑143 (MQ=255)
       ggAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTc                                                                                                                                        >  2:647905/1‑143 (MQ=255)
        gAAATCGCACCCTTTGATGAGAACCGGCGGAAAACCAAACACAAGCGCGTAGATGTTCGCCCGACTACCCGTGGGAAAATTTCCGAAGAAGACACCACCTAAACGCCTATGTATGTCAATCGCGGTGGTCGCCTGACTATTCt                                                                                                                                       <  2:170429/143‑1 (MQ=255)
        gAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTAttat                                                                                                                                       >  2:693774/1‑141 (MQ=255)
          aaTCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCCTTCTga                                                                                                                                     >  1:64951/1‑142 (MQ=255)
          aaTCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTGc                                                                                                                                     >  1:392873/1‑143 (MQ=255)
           aTCGAACCAGTTGATGAGAACCGGCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAAACGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCa                                                                                                                                    <  2:468390/143‑1 (MQ=255)
            tCGAACCGTTGGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGGTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAg                                                                                                                                   <  2:515814/143‑1 (MQ=255)
                              aCCGTCTGAAAAGCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTa                                                                                                                 <  2:1206199/143‑1 (MQ=255)
                                     gAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTGAAGGAAGACCAGGGACACTCACTGActc                                                                                                          >  1:1223822/1‑143 (MQ=255)
                                        aaTCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTg                                                                                                      >  2:356875/1‑143 (MQ=255)
                                             aaCACACACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCTGGt                                                                                                 <  2:305805/143‑1 (MQ=255)
                                               cacaAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGa                                                                                               >  2:872253/1‑143 (MQ=255)
                                                acaAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAg                                                                                              >  2:890893/1‑143 (MQ=255)
                                                 caAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGCCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGCCAATTACTGACTCTTGCCGGTGAGc                                                                                             <  2:1221842/143‑1 (MQ=255)
                                                  aaaCTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGcc                                                                                            >  2:843102/1‑143 (MQ=255)
                                                  aaaCTCTCATATGTCCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGcc                                                                                            <  2:616836/143‑1 (MQ=255)
                                                    acgcgcATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTg                                                                                          <  2:1222278/138‑1 (MQ=255)
                                                      tctcATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTGCAGGAAGACCAGGGACAGTTACTGACTCTTGCCGGTGAGCCTGAc                                                                                        >  1:32544/1‑143 (MQ=255)
                                                       tcgacatGTCCGCCCGACTGACCGTGCAACATTTACCGAAGAAGCCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGAACAATTACTGACTCTTGCCGGTGAGCCTGACg                                                                                       <  1:37446/138‑1 (MQ=255)
                                                         cATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTACAAAAAAAACAGGGACACTTACTGACTCTCGCCGGAGAGCCTGACGGa                                                                                     >  1:434867/1‑143 (MQ=255)
                                                          atatGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaa                                                                                    >  2:154150/1‑143 (MQ=255)
                                                           tatGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaaa                                                                                   <  2:1009753/143‑1 (MQ=255)
                                                             tGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAAc                                                                                 >  2:855092/1‑143 (MQ=255)
                                                              gTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACt                                                                                >  2:39455/1‑143 (MQ=255)
                                                                 cGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcg                                                                             <  2:323819/143‑1 (MQ=255)
                                                                  gCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcgc                                                                            <  2:784267/143‑1 (MQ=255)
                                                                    ccGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTACAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACCCCgcgc                                                                          >  1:649222/1‑143 (MQ=255)
                                                                       aCTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTGCAAAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCgcacc                                                                       >  1:405832/1‑141 (MQ=255)
                                                                                 ggAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTat                                                             >  2:741843/1‑142 (MQ=255)
                                                                                 ggAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTaa                                                             >  2:1138310/1‑143 (MQ=255)
                                                                                  gAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGATCATTATAAAGAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGAAAAACTCCGCGCAGCAGGTCATTAat                                                            >  1:32425/1‑143 (MQ=255)
                                                                                       tttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAt                                                       >  2:901657/1‑143 (MQ=255)
                                                                                         tCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGCAAGTTTTAAATCTTCTCCAGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTc                                                     >  2:1240158/1‑143 (MQ=255)
                                                                                          ccGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCt                                                    >  2:920451/1‑143 (MQ=255)
                                                                                           cGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCtt                                                   >  2:1172352/1‑143 (MQ=255)
                                                                                             aagaagACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaa                                                 >  2:742043/1‑143 (MQ=255)
                                                                                             aagaagACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGTACCATTCTAAAGAAAAACCCGGGACCATTACTGACTCTTGCCGGTGAACCTCACGGAAAACTCCGCGCAGCAGGTCATTAATCCCATTTttac                                                 >  1:865353/1‑142 (MQ=255)
                                                                                              agaagaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaat                                                <  1:855092/143‑1 (MQ=255)
                                                                                                  gaCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAAc                                            <  2:223088/143‑1 (MQ=255)
                                                                                                     accaccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGAAAGAAAACAAAGCGCAGCAGGTCATTAATATCATTCTTAATTAACTaa                                         >  2:35747/1‑143 (MQ=255)
                                                                                                        accGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGACCATTCAACAAAAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAatta                                      >  1:286617/1‑143 (MQ=255)
                                                                                                        accGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCCGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGGGCAGCAGGTCGTTAATATCATTCTTAATTAACtcgtta                                      >  2:75624/1‑138 (MQ=255)
                                                                                                        accGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAatta                                      >  2:332059/1‑143 (MQ=255)
                                                                                                         ccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGTCTCTTGCCGGTGAGCCTGACGGAAAACTCAGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAattat                                     >  2:114775/1‑143 (MQ=255)
                                                                                                         ccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAattat                                     <  1:356875/143‑1 (MQ=255)
                                                                                                          cGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGATCATTAAAAAGAAAAAACACGGTTAATTACTGACTCTTGCCGGTGAGACTGACGAAAAACTCAGCGAAGCAGGTCATAAATCACAATCATAattaattaattaga                                    >  1:1156862/1‑141 (MQ=255)
                                                                                                             aaCCCCGATGTATGTCAACCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAAttatttat                                 <  2:136642/143‑1 (MQ=255)
                                                                                                                gCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATCATTTAtctc                              >  1:1175717/1‑143 (MQ=255)
                                                                                                                  cGATTTATGCCAATCGCGGTGCTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtcat                            <  2:339197/143‑1 (MQ=255)
                                                                                                                  cGATGTATGTCAATCTCGGTGGTCGCCTGATCATTCAGCAAAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGACAACTCCGCGCCGCAGGTCACTAATATCATGATCACTTAACTAATGATTCGTCtcat                            >  1:217018/1‑143 (MQ=255)
                                                                                                                  cGATGTATGTCAATCGCGGTGGTCGCCTGACGATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtcat                            <  2:907887/143‑1 (MQ=255)
                                                                                                                      gtatgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCCGCTGGTCCTTACTATCATTCTTAATTAACTAATTATTTAAATCATCACt                        >  2:352214/1‑143 (MQ=255)
                                                                                                                      gtatgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTAATCACt                        >  2:897683/1‑143 (MQ=255)
                                                                                                                      gtatgtCAATCGCGGTGGTCGCCGGACTATTAAAAAAAAAGACCAGGGACAATTACTGACACTTGCCGGGGAGCCTGCCGGAAAACTCAGCGCAGCAGGTCATAAAAAGCAATCTTAATTAACTAAATAATTAACATATCAcg                        >  1:775967/1‑142 (MQ=255)
                                                                                                                      gtacgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACt                        <  2:376374/143‑1 (MQ=255)
                                                                                                                         tgcaaaTCGCCGGGGGTGCCTCACCATTCTCCAGGAAGACCAGAGACAATTGCTGACTCTAGCCGATGTGCCTGACTGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGaa                     <  1:1013491/139‑1 (MQ=255)
                                                                                                                            caaaagccaagatCGCCTGGCTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGACTAACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAatat                  <  1:133357/131‑1 (MQ=255)
                                                                                                                            cAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTCATTAACTAATTATTTATCTCATCACTGCatat                  >  2:706173/1‑143 (MQ=255)
                                                                                                                             aaTCTCGGTGGTCGCCGTATCATTAAACAAAAACACCATGGACAATTACTGACTCTTGCCGGTGACCCTGCCGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTCATTACCTAATTCTTCATCTCCTCACTGAatacc                 >  1:907068/1‑141 (MQ=255)
                                                                                                                             aaTCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATc                 <  1:154150/143‑1 (MQ=255)
                                                                                                                             aaTCGCGGTGGTCGCCGGACTATTCTGCAGGAAGACCAGGGACAATTACTGCCTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAAGTATTTATCTCATCACTGAATATc                 >  1:476089/1‑143 (MQ=255)
                                                                                                                                  cgGTGGTCGCCTGACTATTCTGCAGGAAGACCAGCGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAat            <  2:112844/143‑1 (MQ=255)
                                                                                                                                            cTGACTGCTCGGCAGGAAGACCAGGGAAAATACATTTCGCTTTCCGTTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGa  <  1:1209688/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGA  >  NC_000913/2473741‑2474023

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