Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,191,651 G→T 5.1% E340D (GAG→GAT)  yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,191,6510GT 5.1% 282.3 / 4.9 79E340D (GAG→GAT) yqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (35/40);  new base T (2/2);  total (37/42)
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.30e-01

TGCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCAT  >  NC_000913/3191503‑3191794
                                                                                                                                                    |                                                                                                                                               
tGCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAg                                                                                                                                                 <  3:750595/149‑1 (MQ=255)
tGCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAg                                                                                                                                                 <  4:716953/149‑1 (MQ=255)
 gCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAgg                                                                                                                                                >  2:796715/1‑149 (MQ=255)
 gCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAgg                                                                                                                                                >  4:172293/1‑149 (MQ=255)
 gCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTAGAgg                                                                                                                                                <  2:104895/149‑1 (MQ=255)
       gAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCACTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGt                                                                                                                                          >  4:536890/1‑149 (MQ=255)
        aaGATATTTCATAACGATGGAACAGCAACGGAATTTAACCACTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGtt                                                                                                                                         <  3:536890/149‑1 (MQ=255)
          gATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTa                                                                                                                                       <  1:468183/149‑1 (MQ=255)
          gATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTa                                                                                                                                       <  2:80835/149‑1 (MQ=255)
              tttCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTaa                                                                                                                                   >  2:480268/1‑149 (MQ=255)
              tttCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTaa                                                                                                                                   >  1:580318/1‑149 (MQ=255)
                          ggAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCt                                                                                                                       >  4:777479/1‑149 (MQ=255)
                             aCAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGaa                                                                                                                    <  2:340948/149‑1 (MQ=255)
                             aCAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGaa                                                                                                                    <  2:116944/149‑1 (MQ=255)
                                  aaCGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAg                                                                                                                                                 <  3:97554/115‑1 (MQ=255)
                                  aaCGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAg                                                                                                                                                 >  4:97554/1‑115 (MQ=255)
                                  aaCGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCa                                                                                                               >  4:42661/1‑149 (MQ=255)
                                   aCGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCAc                                                                                                              >  2:225916/1‑149 (MQ=255)
                                      gAATTTAACTGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCta                                                                                                           >  4:511615/1‑149 (MQ=255)
                                      gAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATc                                                                                                                >  1:553453/1‑144 (MQ=255)
                                      gAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATc                                                                                                                <  2:553453/144‑1 (MQ=255)
                                       aaTTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATcc                                                                                                                        >  3:667452/1‑135 (MQ=255)
                                       aaTTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATcc                                                                                                                        <  4:667452/135‑1 (MQ=255)
                                          ttAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAAc                                                                                                       <  1:590866/149‑1 (MQ=255)
                                            aaCCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTa                                                                                                     >  3:316762/1‑149 (MQ=255)
                                                 cTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATc                                                                                                <  2:580318/149‑1 (MQ=255)
                                                  ttCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTa                                                                                                     <  1:632167/143‑1 (MQ=255)
                                                  ttCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTa                                                                                                     >  2:632167/1‑143 (MQ=255)
                                                  ttCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCt                                                                                               <  2:638769/149‑1 (MQ=255)
                                                    ctcACGATTTGCTTCTTTTTCATCGTCTAGTGCACCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTAGAGGGACCGTTTAGTAAAGACGTCATCCTGAAAATCACCTATAACTAGTATCTaa                                                                                             >  1:231529/1‑149 (MQ=255)
                                                    ctcACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTaa                                                                                             <  1:71886/149‑1 (MQ=255)
                                                     tcACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAg                                                                                                    >  4:242775/1‑141 (MQ=255)
                                                     tcACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAg                                                                                                    <  3:242775/141‑1 (MQ=255)
                                                     tcACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAAt                                                                                            >  3:800095/1‑149 (MQ=255)
                                                      cacaATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATa                                                                                           <  4:781087/149‑1 (MQ=255)
                                                       acaATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATcc                                                                                                                        <  3:515580/119‑1 (MQ=255)
                                                       acaATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATcc                                                                                                                        >  4:515580/1‑119 (MQ=255)
                                                       acaATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATAc                                                                                          >  4:145268/1‑149 (MQ=255)
                                                         aaTTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCtat                                                                                                          >  4:24012/1‑131 (MQ=255)
                                                         aaTTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCtat                                                                                                          <  3:24012/131‑1 (MQ=255)
                                                          aTTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACaaa                                                                                       <  1:505768/149‑1 (MQ=255)
                                                          aTTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACaaa                                                                                       >  3:801514/1‑149 (MQ=255)
                                                           tttGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAac                                                                                      <  3:503221/149‑1 (MQ=255)
                                                              gCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACt                                                                                   <  3:169748/149‑1 (MQ=255)
                                                                ttcttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACaaa                                                                                       <  3:723148/143‑1 (MQ=255)
                                                                ttcttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACaaa                                                                                       >  4:723148/1‑143 (MQ=255)
                                                                ttcttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTat                                                                                 >  3:117592/1‑148 (MQ=255)
                                                                ttcttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTaa                                                                                 >  4:458226/1‑149 (MQ=255)
                                                                   tttttCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGATGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAAc                                                                              >  1:513970/1‑149 (MQ=255)
                                                                   tttttCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGATGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAAc                                                                              >  1:70200/1‑149 (MQ=255)
                                                                          tCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCAt                                                                       <  1:586506/149‑1 (MQ=255)
                                                                           cGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATc                                                                      <  1:225916/149‑1 (MQ=255)
                                                                                 gTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCAcc                                                                                                             <  2:704171/104‑1 (MQ=255)
                                                                                 gTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCAcc                                                                                                             >  1:704171/1‑104 (MQ=255)
                                                                                    cgcCTTCGGTCACCGCACCATTCAGTGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGtt                                                             <  3:172293/149‑1 (MQ=255)
                                                                                    cgcCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGtt                                                             <  1:796715/149‑1 (MQ=255)
                                                                                       cTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTaa                                                                                             <  2:77830/114‑1 (MQ=255)
                                                                                       cTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTaa                                                                                             >  1:77830/1‑114 (MQ=255)
                                                                                        ttCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTatat                                                         <  3:42661/149‑1 (MQ=255)
                                                                                           ggTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTATAACGGGCCTGATGGCCCGTtttagtgt                                                      <  4:117592/149‑5 (MQ=255)
                                                                                            gTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGt                                                     <  4:801514/149‑1 (MQ=255)
                                                                                               aCCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTAGCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGct                                                  <  1:480268/149‑1 (MQ=255)
                                                                                                 cGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                          <  1:232208/123‑1 (MQ=255)
                                                                                                 cGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                          >  2:232208/1‑123 (MQ=255)
                                                                                                   cACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTaa                                              >  3:763249/1‑149 (MQ=255)
                                                                                                    aCCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAAc                                             <  3:458226/149‑1 (MQ=255)
                                                                                                        ttCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                         <  3:777479/149‑1 (MQ=255)
                                                                                                        ttCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                         >  2:299492/1‑149 (MQ=255)
                                                                                                              gCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctctc                                   >  4:126683/1‑149 (MQ=255)
                                                                                                              gCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctctc                                   >  4:274849/1‑149 (MQ=255)
                                                                                                                aGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCtt                                 <  4:97924/149‑1 (MQ=255)
                                                                                                                  aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGc                               >  3:119462/1‑149 (MQ=255)
                                                                                                                         tGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatc                        >  2:635476/1‑149 (MQ=255)
                                                                                                                            gAACCCGGCAGAAACGGTTGTTGATGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCatcatc                     <  2:513970/149‑1 (MQ=255)
                                                                                                                            gAACCCGGCAGAAACGGTTGTTGATGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCatcatc                     <  2:70200/149‑1 (MQ=255)
                                                                                                                            gAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAg                                                                                                    <  4:696967/70‑1 (MQ=255)
                                                                                                                            gAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAg                                                                                                    >  3:696967/1‑70 (MQ=255)
                                                                                                                                ccGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                 >  4:91597/1‑149 (MQ=255)
                                                                                                                                   gCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGAc              <  4:800095/149‑1 (MQ=255)
                                                                                                                                        aaCGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc         >  2:695127/1‑149 (MQ=255)
                                                                                                                                           ggTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTAt                                                                                                 <  3:602616/58‑1 (MQ=255)
                                                                                                                                           ggTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTAt                                                                                                 >  4:602616/1‑58 (MQ=255)
                                                                                                                                           ggTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAg      <  1:299492/149‑1 (MQ=255)
                                                                                                                                               gttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                 <  3:480722/134‑1 (MQ=255)
                                                                                                                                               gttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                 >  4:480722/1‑134 (MQ=255)
                                                                                                                                               gttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCAt  <  3:91597/149‑1 (MQ=255)
                                                                                                                                                    gggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                 <  1:703668/129‑1 (MQ=255)
                                                                                                                                                    gggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGtt                 >  2:703668/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGCTATTGAAGATATTTCATAACGATGGAACAGCAACGGAATTTAACCGCTTCTCACAATTTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCAT  >  NC_000913/3191503‑3191794

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