Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,105,384 G→T 6.9% intergenic (+91/‑30) csgC → / → ymdA curli assembly protein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,105,3840GT 6.9% 253.2 / 7.5 72intergenic (+91/‑30)csgC/ymdAcurli assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base T (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGG  >  NC_000913/1105239‑1105528
                                                                                                                                                 |                                                                                                                                                
tGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGAt                                                                                                                                               <  2:126955/149‑1 (MQ=255)
tGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGAt                                                                                                                                               >  1:442367/1‑149 (MQ=255)
 ggACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATg                                                                                                                                              >  4:70001/1‑149 (MQ=255)
  gACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGa                                                                                                                                             <  4:110405/149‑1 (MQ=255)
   aCAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGAc                                                                                                                                            >  1:639429/1‑149 (MQ=255)
      gTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGAttagtca                                                                                                                                         >  3:277026/1‑143 (MQ=255)
         aCTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCa                                                                                                                                      <  3:37072/149‑1 (MQ=255)
           ttCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCAtt                                                                                                                                    <  3:102636/149‑1 (MQ=255)
           ttCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCAtt                                                                                                                                    <  4:528934/149‑1 (MQ=255)
             cATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGa                                                                                                                                  <  4:338386/149‑1 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACGCAACCCTCAAGGATGACTAATCATTGAg                                                                                                                                 <  3:425727/149‑1 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAg                                                                                                                                 >  4:628880/1‑149 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAg                                                                                                                                 <  2:325800/149‑1 (MQ=255)
               tttATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAgg                                                                                                                                <  4:277026/149‑1 (MQ=255)
                    cacaACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATa                                                                                                                           <  4:634127/149‑1 (MQ=255)
                     acaacaATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAg                                                                                                                          >  4:218047/1‑149 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        >  2:34692/1‑149 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        <  4:343611/149‑1 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        <  3:218047/149‑1 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        <  2:442367/149‑1 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        >  2:113304/1‑149 (MQ=255)
                       aacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                        >  1:623299/1‑149 (MQ=255)
                          aaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAg                                                                                                                                 >  4:103726/1‑137 (MQ=255)
                          aaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAg                                                                                                                                 <  3:103726/137‑1 (MQ=255)
                           aTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATaaa                                                                                                                    <  3:612436/149‑1 (MQ=255)
                               ccgccCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGtt                                                                                                                <  2:639429/149‑1 (MQ=255)
                                   ccTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGa                                                                                                            <  4:37488/149‑1 (MQ=255)
                                   ccTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGa                                                                                                            >  1:641892/1‑149 (MQ=255)
                                    ctctTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGAc                                                                                                           <  3:70001/149‑1 (MQ=255)
                                       ttCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCa                                                                                                                                      <  4:348021/119‑1 (MQ=255)
                                       ttCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCa                                                                                                                                      >  3:348021/1‑119 (MQ=255)
                                       ttCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGa                                                                                                                                  <  2:433302/123‑1 (MQ=255)
                                       ttCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGa                                                                                                                                  >  1:433302/1‑123 (MQ=255)
                                       ttCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCtt                                                                                                        <  4:503660/149‑1 (MQ=255)
                                          aGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCttttt                                                                                                     <  2:682529/149‑1 (MQ=255)
                                                  cTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAActct                                                                                             >  4:553147/1‑149 (MQ=255)
                                                   ttAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAActctc                                                                                            >  3:144256/1‑149 (MQ=255)
                                                   ttAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAActctc                                                                                            >  1:58260/1‑149 (MQ=255)
                                                             gAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCgg                                                                                  >  1:177553/1‑149 (MQ=255)
                                                              aaaTATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGc                                                                                 >  2:407368/1‑149 (MQ=255)
                                                                  atCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAgtt                                                                             <  3:553147/149‑1 (MQ=255)
                                                                      aGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAgtttgtt                                                                         <  1:476039/149‑1 (MQ=255)
                                                                        cATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGTAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGtttt                                                                       <  1:34692/149‑1 (MQ=255)
                                                                               tGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCttt                                                                >  3:74902/1‑149 (MQ=255)
                                                                                gAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCtttt                                                               >  4:248536/1‑149 (MQ=255)
                                                                                 aaTCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTa                                                              >  3:480169/1‑149 (MQ=255)
                                                                                        agagaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCAt                                                       >  1:270930/1‑149 (MQ=255)
                                                                                         gagaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCAtt                                                      >  3:271025/1‑149 (MQ=255)
                                                                                            aaTGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCt                                                   >  2:622417/1‑149 (MQ=255)
                                                                                                     ctgcGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCAcc                                          >  3:530427/1‑149 (MQ=255)
                                                                                                       gcGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGt                                        <  1:407368/149‑1 (MQ=255)
                                                                                                        cGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTg                                       <  2:641892/149‑1 (MQ=255)
                                                                                                           aTCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCaa                                    <  2:489959/149‑1 (MQ=255)
                                                                                                            tCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAAACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGaa                                               <  3:424832/137‑1 (MQ=255)
                                                                                                            tCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAAACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGaa                                               >  4:424832/1‑137 (MQ=255)
                                                                                                              aGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAggg                                 <  3:248536/149‑1 (MQ=255)
                                                                                                               gCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAgggg                                >  3:188530/1‑149 (MQ=255)
                                                                                                                 cAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGc                              <  2:240497/149‑1 (MQ=255)
                                                                                                                      tGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCTTGAt                         <  4:458188/149‑1 (MQ=255)
                                                                                                                        aaaGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATcc                       <  4:144256/149‑1 (MQ=255)
                                                                                                                            tATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAgggg                                >  1:148059/1‑136 (MQ=255)
                                                                                                                            tATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAgggg                                <  2:148059/136‑1 (MQ=255)
                                                                                                                            tATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCAttt                   >  4:278118/1‑149 (MQ=255)
                                                                                                                               ggATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTa                <  2:623299/149‑1 (MQ=255)
                                                                                                                                gATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCt                                                   >  3:548682/1‑113 (MQ=255)
                                                                                                                                gATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCt                                                   <  4:548682/113‑1 (MQ=255)
                                                                                                                                   aaCACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGc            <  2:177553/149‑1 (MQ=255)
                                                                                                                                        aaCCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAAt       >  3:363642/1‑149 (MQ=255)
                                                                                                                                           ccTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATtgt    <  3:628880/149‑1 (MQ=255)
                                                                                                                                             tCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTgg  >  2:6888/1‑149 (MQ=255)
                                                                                                                                             tCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTgg  >  3:654158/1‑149 (MQ=255)
                                                                                                                                             tCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTgg  >  4:211628/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGG  >  NC_000913/1105239‑1105528

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