Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,184,692 +G intergenic (‑226/‑148) ribB ← / → ubiK 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/ubiquinone biosynthesis accessory factor UbiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,184,6881.G96.7% 205.5 / ‑0.3 60intergenic (‑222/‑152)ribB/ubiK3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/ubiquinone biosynthesis accessory factor UbiK
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (20/38);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTAT  >  NC_000913/3184552‑3184821
                                                                                                                                         |                                                                                                                                     
gTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTggg                                                                                                                                    <  1:728439/140‑1 (MQ=255)
      gCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTgggggaatt                                                                                                                              <  1:789614/140‑6 (MQ=255)
          gACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgc                                                                                                                           >  2:89644/1‑140 (MQ=255)
          gACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgc                                                                                                                           >  1:318908/1‑140 (MQ=255)
          gACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgc                                                                                                                           >  1:745595/1‑140 (MQ=255)
           aCAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgcc                                                                                                                          <  2:318908/140‑1 (MQ=255)
           aCAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgcc                                                                                                                          <  2:158055/140‑1 (MQ=255)
           aCAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgcc                                                                                                                          <  1:265977/140‑1 (MQ=255)
           aCAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgcc                                                                                                                          <  1:794812/140‑1 (MQ=255)
            cAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgccc                                                                                                                         <  2:975026/140‑1 (MQ=255)
            cAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTcgccc                                                                                                                         <  1:636543/140‑1 (MQ=255)
                  cGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCt                                                                                                                   <  2:662276/140‑1 (MQ=255)
                  cGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCt                                                                                                                   <  1:976706/140‑1 (MQ=255)
                           ctcCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAg                                                                                                          >  2:1141341/1‑140 (MQ=255)
                           ctcCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAg                                                                                                          <  1:832177/140‑1 (MQ=255)
                           ctcCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCACACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAg                                                                                                          <  1:1093671/140‑1 (MQ=255)
                            tcCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGc                                                                                                         <  1:1048026/140‑1 (MQ=255)
                                aTCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAt                                                                                                     <  2:764733/140‑1 (MQ=255)
                                 tCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                    <  2:75243/140‑1 (MQ=255)
                                 tCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                    >  2:396412/1‑140 (MQ=255)
                                 tCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTCCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                    <  2:635217/140‑1 (MQ=255)
                                 tCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCgg                                                                                                       >  1:1216056/1‑136 (MQ=255)
                                 tCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCgg                                                                                                       <  2:1216056/136‑1 (MQ=255)
                                      aCTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCTGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACta                                                                                               >  2:918056/1‑140 (MQ=255)
                                      aCTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACta                                                                                               >  2:623335/1‑140 (MQ=255)
                                      aCTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACta                                                                                               >  1:962483/1‑140 (MQ=255)
                                       ctctAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACtat                                                                                              <  1:60221/140‑1 (MQ=255)
                                       ctctAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACtat                                                                                              <  1:1115320/140‑1 (MQ=255)
                                         ctAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATaa                                                                                            <  2:807483/140‑1 (MQ=255)
                                         ctAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATaa                                                                                            <  2:745595/140‑1 (MQ=255)
                                          tAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAAc                                                                                           >  2:424046/1‑140 (MQ=255)
                                           aaCCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACg                                                                                          <  1:971832/140‑1 (MQ=255)
                                               gTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTaa                                                                                      >  1:139593/1‑140 (MQ=255)
                                                    ccccGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGAtt                                                                                 >  2:698201/1‑140 (MQ=255)
                                                    ccccGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGAtt                                                                                 >  2:202274/1‑140 (MQ=255)
                                                         gAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCgg                                                                            <  2:1093560/140‑1 (MQ=255)
                                                             tACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTATTGATTAGCGtcacca                                                                       <  2:152819/140‑1 (MQ=255)
                                                                acCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAAc                                                                     <  2:139593/140‑1 (MQ=255)
                                                                    gATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCAt                                                                 >  1:155890/1‑140 (MQ=255)
                                                                        tgctgTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCTGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCt                                                             <  1:918056/140‑1 (MQ=255)
                                                                               cTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATa                                                      >  1:592678/1‑140 (MQ=255)
                                                                                  tGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAAtt                                                   >  2:761403/1‑140 (MQ=255)
                                                                                   gAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTc                                                  <  2:1154547/140‑1 (MQ=255)
                                                                                   gAGTTCGCACCCAAAACGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTc                                                  <  1:530279/140‑1 (MQ=255)
                                                                                       tCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAAt                                                    <  2:475502/134‑1 (MQ=255)
                                                                                       tCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAAt                                                    >  1:475502/1‑134 (MQ=255)
                                                                                       tCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGt                                              >  2:1120491/1‑140 (MQ=255)
                                                                                             ccAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGa                                        <  1:89644/140‑1 (MQ=255)
                                                                                              cAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGAc                                       <  2:610929/140‑1 (MQ=255)
                                                                                               aaaGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACt                                      <  1:761403/140‑1 (MQ=255)
                                                                                                  gcgcTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTAc                                   <  1:1141341/140‑1 (MQ=255)
                                                                                                     cTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGa                                        >  2:426768/1‑132 (MQ=255)
                                                                                                     cTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGa                                        <  1:426768/132‑1 (MQ=255)
                                                                                                     cTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCt                                >  2:1040969/1‑140 (MQ=255)
                                                                                                       cgcgGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTa                              <  1:698201/140‑1 (MQ=255)
                                                                                                       cgcgGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTa                              <  1:992284/140‑1 (MQ=255)
                                                                                                          gggCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCg                           <  2:963339/140‑1 (MQ=255)
                                                                                                           ggCTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGc                          <  2:589780/140‑1 (MQ=255)
                                                                                                             cTTTCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCAc                        <  1:202274/140‑1 (MQ=255)
                                                                                                              tttCAACTGAGTTGATTTACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACt                       <  2:962483/140‑1 (MQ=255)
                                                                                                                              ttACCGCCGGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAaca       <  1:623335/140‑1 (MQ=255)
                                                                                                                                   gccgGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTAt  >  2:397259/1‑140 (MQ=255)
                                                                                                                                   gccgGTGGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTAt  <  1:396412/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTAT  >  NC_000913/3184552‑3184821

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