Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,152 G→A 6.0% intergenic (+40/‑168) ihfB → / → ycaI integration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1520GA 6.0% 382.4 / 8.3 116intergenic (+40/‑168)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (60/49);  new base A (2/5);  total (62/54)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 119 total aligned reads displayed.

ACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACTTTC  >  NC_000913/964005‑964300
                                                                                                                                                   |                                                                                                                                                    
aCCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                                                                                                     <  4:625368/149‑1 (MQ=255)
   ggACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTg                                                                                                                                                  >  2:644830/1‑149 (MQ=255)
     aCGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                                                                                <  4:411373/149‑1 (MQ=255)
        tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCCCCCGACAGGTGCtttt                                                                                                                                             <  4:247955/149‑1 (MQ=255)
        tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtttt                                                                                                                                             <  1:239456/149‑1 (MQ=255)
          aTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAAAAAAGCACCCGACAGGTGCttttct                                                                                                                                           <  4:927179/149‑3 (MQ=255)
               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                                                                      <  2:710986/149‑1 (MQ=255)
                      gCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAAt                                                                                                                               <  4:809835/149‑1 (MQ=255)
                        gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                                                                              >  1:949925/1‑148 (MQ=255)
                        gATAAAGTAGAACTCGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                                                                             >  4:260355/1‑149 (MQ=255)
                            aaGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaa                                                                                                                         <  4:333243/149‑1 (MQ=255)
                             aGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                                                                        <  1:359210/149‑1 (MQ=255)
                             aGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                                                                        >  4:604376/1‑149 (MQ=255)
                                aGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgg                                                                                                                     >  3:111034/1‑149 (MQ=255)
                                     tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                                                                                >  3:265044/1‑117 (MQ=255)
                                     tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                                                                                <  4:265044/117‑1 (MQ=255)
                                     tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCt                                                                                                                >  1:767145/1‑149 (MQ=255)
                                     tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCt                                                                                                                <  4:141815/149‑1 (MQ=255)
                                      ggaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTg                                                                                                               >  4:101536/1‑149 (MQ=255)
                                      ggaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTg                                                                                                               <  3:820750/149‑1 (MQ=255)
                                       gaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGTGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtct                                                                                                                                    <  4:660082/127‑1 (MQ=255)
                                       gaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtct                                                                                                                                    >  3:660082/1‑127 (MQ=255)
                                       gaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGc                                                                                                              >  2:977666/1‑149 (MQ=255)
                                         aggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCt                                                                                                            >  3:649586/1‑149 (MQ=255)
                                         aggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACATGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCt                                                                                                            >  2:143179/1‑149 (MQ=255)
                                            aaaaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCCCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGc                                                                                                         >  4:863153/1‑149 (MQ=255)
                                              aaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCag                                                                                                       >  4:446617/1‑149 (MQ=255)
                                               aTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCaga                                                                                                      <  4:423207/149‑1 (MQ=255)
                                                   gTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                                  >  2:252313/1‑149 (MQ=255)
                                                    ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                                 <  2:966051/149‑1 (MQ=255)
                                                    ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                                 <  1:1011882/149‑1 (MQ=255)
                                                    ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                                 <  1:878506/149‑1 (MQ=255)
                                                    gtCCTCACTTTAAACCTGGTAAATAACTGCGCGATCGCGCCAATATTTACGGTTAAGATTTTTACTCAAACTCGAACGACAGAAAAGCACCTGTCGGGTGCTTTTTTAATTACTCTAATCTGGAACCGGAAGCTGGCTCGCAGAGttta                                                                                                 <  4:908813/148‑2 (MQ=255)
                                                     tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                                                                                <  1:424559/149‑1 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGaa                                                                                              <  3:609107/148‑1 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                             >  3:1041326/1‑149 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                             >  4:854975/1‑149 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                             >  2:620639/1‑149 (MQ=255)
                                                          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAg                                                                                           >  1:357491/1‑149 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGcc                                                                                                             >  2:1047602/1‑129 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGcc                                                                                                             <  1:1047602/129‑1 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtt                                                                                         >  2:274299/1‑149 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtt                                                                                         >  2:927931/1‑149 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctc                                                                                                                                          <  4:225210/100‑1 (MQ=39)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctc                                                                                                                                          >  3:225210/1‑100 (MQ=39)
                                                              tAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtttt                                                                                       <  1:274299/149‑1 (MQ=255)
                                                               aaaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTc                                                                                      >  3:61570/1‑149 (MQ=255)
                                                                   cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAccc                                                                                  >  4:908323/1‑149 (MQ=255)
                                                                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACa                                                                                            <  3:597823/138‑1 (MQ=255)
                                                                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCt                                                                                 >  2:464008/1‑149 (MQ=255)
                                                                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTAAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACa                                                                                            >  4:597823/1‑138 (MQ=255)
                                                                       tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                                                                              >  4:935147/1‑149 (MQ=255)
                                                                       tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                                                                              >  1:559218/1‑149 (MQ=255)
                                                                        aaaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                                                                             >  3:1039623/1‑149 (MQ=255)
                                                                          aGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGtt                                                                           >  2:753566/1‑149 (MQ=255)
                                                                                         cgcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAata                                                            <  1:329181/149‑1 (MQ=255)
                                                                                          gcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat                                                           >  1:10649/1‑149 (MQ=255)
                                                                                            cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc                                                         <  1:252313/149‑1 (MQ=255)
                                                                                             aaTATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGcc                                                        >  1:1058437/1‑149 (MQ=255)
                                                                                                 tttACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCgg                                                    >  3:513485/1‑149 (MQ=255)
                                                                                                     cGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGaa                                                >  2:378830/1‑149 (MQ=255)
                                                                                                      ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAAc                                               <  1:378830/149‑1 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  4:548400/1‑149 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  2:247532/1‑149 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  2:500936/1‑149 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  4:359170/1‑149 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  3:572082/1‑149 (MQ=255)
                                                                                                        ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTGTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                             >  1:1032389/1‑149 (MQ=255)
                                                                                                            gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCaaaa                                         <  2:559218/149‑1 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                                 <  2:816449/92‑1 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                                 >  1:816449/1‑92 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAAt                                        <  2:360147/149‑1 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAAt                                        >  4:944818/1‑149 (MQ=255)
                                                                                                              tttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAAtt                                       <  2:870925/149‑1 (MQ=255)
                                                                                                                tttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTcc                                     <  1:644830/149‑1 (MQ=255)
                                                                                                                    cTCAAACTTGAACGAGAGAAGAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACGGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                                 <  1:143179/149‑1 (MQ=255)
                                                                                                                    cTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                                 >  1:322098/1‑149 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGtt                                                                           >  4:263724/1‑106 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGtt                                                                           <  3:263724/106‑1 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacca                                <  1:642085/149‑1 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                                                                              >  1:70052/1‑53 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                              >  2:204783/1‑149 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                              <  4:137840/149‑1 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                              <  3:854975/149‑1 (MQ=255)
                                                                                                                           ttGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAAc                          >  3:353252/1‑149 (MQ=255)
                                                                                                                             gAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCg                        >  2:490327/1‑149 (MQ=255)
                                                                                                                               aCGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCt                      <  1:435016/149‑1 (MQ=255)
                                                                                                                                 gagagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTca                    >  1:307138/1‑149 (MQ=255)
                                                                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac                  >  1:710643/1‑149 (MQ=255)
                                                                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac                  >  3:367136/1‑149 (MQ=255)
                                                                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac                  >  3:960437/1‑149 (MQ=255)
                                                                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac                  <  4:513485/149‑1 (MQ=255)
                                                                                                                                      aaaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACaaa               <  1:464008/149‑1 (MQ=255)
                                                                                                                                      aaaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACaaa               <  1:1009836/149‑1 (MQ=255)
                                                                                                                                        aaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGt             >  4:43596/1‑149 (MQ=255)
                                                                                                                                         aGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGtt            <  3:68043/149‑1 (MQ=255)
                                                                                                                                            aCCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGca         <  4:58959/149‑1 (MQ=255)
                                                                                                                                            aCCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGca         >  3:31943/1‑149 (MQ=255)
                                                                                                                                                   gggTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACTTTc  <  4:1039623/149‑1 (MQ=255)
                                                                                                                                                   gggTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACTTTc  <  2:194929/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACACTTTC  >  NC_000913/964005‑964300

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