Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 115,331 C→T G132S (GGT→AGT)  hofC ← assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913115,3310CT98.7% 308.4 / ‑3.9 77G132S (GGT→AGT) hofCassembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (42/34);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

AGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACCCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTCCGCCA  >  NC_000913/115183‑115474
                                                                                                                                                    |                                                                                                                                               
agcaTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCg                                                                                                                                               <  1:770399/151‑1 (MQ=255)
  caTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTg                                                                                                                                             <  2:310579/151‑1 (MQ=255)
   aTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGc                                                                                                                                            >  1:1277263/1‑151 (MQ=255)
      gCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGa                                                                                                                                         <  1:470124/151‑1 (MQ=255)
          caACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGatcat                                                                                                                                     <  2:1323152/151‑1 (MQ=255)
                    gATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGAt                                                                                                                             >  1:1101169/1‑149 (MQ=255)
                     aTTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGAt                                                                                                                             >  1:692584/1‑148 (MQ=255)
                        gCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAgc                                                                                                                        <  1:992059/150‑1 (MQ=255)
                          cATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAgcggc                                                                                                                     <  1:510486/151‑1 (MQ=255)
                            tCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGcggcgg                                                                                                                   >  2:151102/1‑151 (MQ=255)
                            tCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGcggcgg                                                                                                                   >  1:560710/1‑151 (MQ=255)
                              gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGaa                                                                                                                 <  2:361514/151‑1 (MQ=255)
                              gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGa                                                                                                                  >  1:1150560/1‑150 (MQ=255)
                                 aaaaTGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAAt                                                                                                               <  2:652930/150‑1 (MQ=255)
                                 aaaaTGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAAt                                                                                                               >  1:555303/1‑150 (MQ=255)
                                  aaaTGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATAc                                                                                                             >  1:1380517/1‑151 (MQ=255)
                                      gatgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACctct                                                                                                         >  1:794929/1‑151 (MQ=255)
                                      gatgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACctc                                                                                                          <  2:439502/150‑1 (MQ=255)
                                       atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTg                                                                                                        <  1:846503/151‑1 (MQ=255)
                                              gATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAggg                                                                                                 >  2:1437203/1‑151 (MQ=255)
                                                 aacgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGtaa                                                                                              <  1:144397/151‑1 (MQ=255)
                                                  acgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATAc                                                                                                             <  2:1215532/135‑1 (MQ=255)
                                                  acgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATAc                                                                                                             >  1:1215631/1‑135 (MQ=255)
                                                  acgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGtaat                                                                                             <  1:399062/151‑1 (MQ=255)
                                                   cgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGtaat                                                                                             >  1:466359/1‑150 (MQ=255)
                                                   cgtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGta                                                                                               >  1:136367/1‑148 (MQ=255)
                                                    gtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGtaataa                                                                                           <  2:903876/151‑1 (MQ=255)
                                                    gtaacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGtaata                                                                                            >  2:1079962/1‑150 (MQ=255)
                                                      aacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAgg                                                                                         <  2:405403/151‑1 (MQ=255)
                                                      aacgCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAg                                                                                          <  2:503277/150‑1 (MQ=255)
                                                           tGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCCGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTg                                                                                    >  1:212638/1‑151 (MQ=255)
                                                           tGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTg                                                                                    >  2:446846/1‑151 (MQ=255)
                                                            gATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTgg                                                                                   <  2:275997/151‑1 (MQ=255)
                                                            gATTTCACTTTGCCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTg                                                                                    <  2:762793/150‑1 (MQ=255)
                                                              tttCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGa                                                                                  >  1:495459/1‑150 (MQ=255)
                                                               ttCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGaa                                                                                 >  2:1114866/1‑150 (MQ=255)
                                                                   cTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGc                                                                            >  2:1464918/1‑151 (MQ=255)
                                                                   cTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGc                                                                            >  2:1282242/1‑151 (MQ=255)
                                                                      tGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCaa                                                                          <  2:984064/150‑1 (MQ=255)
                                                                      tGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAAt                                                                         >  1:1341307/1‑151 (MQ=255)
                                                                       gTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATg                                                                        <  2:804725/151‑1 (MQ=255)
                                                                         cGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACCCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGcc                                                                      >  2:1449040/1‑151 (MQ=255)
                                                                              aaCTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGtt                                                                 >  2:1025113/1‑151 (MQ=255)
                                                                                  gACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCgag                                                             <  1:1244349/151‑1 (MQ=255)
                                                                                   aCGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCgaga                                                            <  1:843237/151‑1 (MQ=255)
                                                                                   aCGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCg                                                               >  1:1184648/1‑148 (MQ=255)
                                                                                     gCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCg                                                               >  2:1097436/1‑146 (MQ=255)
                                                                                     gCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATc                                                          <  2:166295/151‑1 (MQ=255)
                                                                                       tGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCg                                                         >  2:352620/1‑150 (MQ=255)
                                                                                        gggCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGt                                                        <  1:854191/150‑1 (MQ=255)
                                                                                           cTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTgcgc                                                    >  1:310747/1‑151 (MQ=255)
                                                                                           cTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTgcgc                                                    >  2:1353169/1‑151 (MQ=255)
                                                                                            tttttGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGcc                                                   >  2:29742/1‑151 (MQ=255)
                                                                                                 gCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGa                                              <  1:529491/151‑1 (MQ=255)
                                                                                                      cgcgCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgca                                         <  1:352666/151‑1 (MQ=255)
                                                                                                      cgcgCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgca                                         <  2:1000224/151‑1 (MQ=255)
                                                                                                       gcgcCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgc                                          >  2:887557/1‑149 (MQ=255)
                                                                                                       gcgcCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgc                                          >  1:1260456/1‑149 (MQ=255)
                                                                                                       gcgcCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgc                                          >  1:267452/1‑149 (MQ=255)
                                                                                                        cgcCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgcagc                                       >  1:505381/1‑151 (MQ=255)
                                                                                                        cgcCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTgcag                                        >  2:846287/1‑150 (MQ=255)
                                                                                                               tCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCtt                                >  1:1190861/1‑151 (MQ=255)
                                                                                                                cAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCtt                                <  1:619379/150‑1 (MQ=255)
                                                                                                                cAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCtt                                >  2:557818/1‑150 (MQ=255)
                                                                                                                 aaaGCAGCATTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGc                              >  2:1024393/1‑151 (MQ=255)
                                                                                                                     cagcaTTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCAtt                          <  2:843399/151‑1 (MQ=255)
                                                                                                                     cagcaTTCATCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCAtt                          >  2:155248/1‑151 (MQ=255)
                                                                                                                             aTCCAGCTTACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTg                  >  1:801880/1‑151 (MQ=255)
                                                                                                                                    ttACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCt           >  1:331205/1‑151 (MQ=255)
                                                                                                                                     tACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCt           <  2:495401/150‑1 (MQ=255)
                                                                                                                                     tACCGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCt           <  2:1293379/150‑1 (MQ=255)
                                                                                                                                       ccGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGtt        >  1:212258/1‑151 (MQ=255)
                                                                                                                                        cGGTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTc       <  2:1087283/151‑1 (MQ=255)
                                                                                                                                         ggTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTcc      >  2:542225/1‑151 (MQ=255)
                                                                                                                                          gTCAGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTCCg     >  1:731600/1‑151 (MQ=255)
                                                                                                                                             aGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTCCGCca  <  2:1313922/151‑1 (MQ=255)
                                                                                                                                             aGTTCACTCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTCCGCca  <  2:1392758/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCCAGTTCAAAGCAGCATTCATCCAGCTTACCGGTCAGTTCACCCGTGCGGATCATCGCCTGATAGAGCGGCGGAAATACCTCTGACCAGGGTAATAAGGCATTGGAAAAAGCAATGCCCTGTTCGAGATCGTGCGCCAGCGATTGCAGCAACGCTTGCCATTGCTTACTGGGATGCTGTTCCGCCA  >  NC_000913/115183‑115474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: