Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 380,091 G→T 5.7% A251D (GCC→GAC)  yaiO ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0910GT 5.7% 310.4 / 7.7 88A251D (GCC→GAC) yaiOouter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (45/38);  new base T (3/2);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCA  >  NC_000913/379947‑380239
                                                                                                                                                 |                                                                                                                                                    
ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                                                                                                                   >  3:590489/1‑149 (MQ=255)
      gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGGCAGACGGGCt                                                                                                                                                    >  2:298222/1‑141 (MQ=255)
      gTATATAACACCGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                                                                                                                    <  1:298222/141‑1 (MQ=255)
           taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAAAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                                                                                                        <  4:590489/149‑1 (MQ=255)
              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                                                                                                           >  2:568063/1‑142 (MQ=255)
              cacTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                                                                                                           <  1:568063/142‑1 (MQ=255)
                                                 ggAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACtt                                                                                                   >  4:740570/1‑147 (MQ=255)
                                                 ggAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACtt                                                                                                   <  3:740570/147‑1 (MQ=255)
                                                  gAAAGCAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGTCTGCaac                                                                                                                                               <  3:1617336/102‑2 (MQ=255)
                                                  gAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGTCTGCaac                                                                                                                                               >  4:1617336/1‑101 (MQ=255)
                                                  gAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAc                                                                                                >  4:1764387/1‑149 (MQ=255)
                                                        aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                                                                                                                    <  1:307601/123‑1 (MQ=255)
                                                        aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                                                                                                                    >  2:307601/1‑123 (MQ=255)
                                                        aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg                                                                                          >  2:330903/1‑149 (MQ=255)
                                                         aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGGCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt                                                                                         >  4:1790496/1‑149 (MQ=255)
                                                         aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACttt                                                                                                  >  3:374800/1‑140 (MQ=255)
                                                         aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACttt                                                                                                  <  4:374800/140‑1 (MQ=255)
                                                                        gggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTaa                                                                           <  3:1764387/149‑1 (MQ=255)
                                                                          gggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                                                                                                                   <  3:1654543/107‑1 (MQ=255)
                                                                          gggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                                                                                                                   >  4:1654543/1‑107 (MQ=255)
                                                                           ggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGtt                                                                        <  1:330903/149‑1 (MQ=255)
                                                                            gCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTg                                                                       >  2:1389692/1‑149 (MQ=255)
                                                                             cTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                                                                                                                  <  4:807924/73‑1 (MQ=255)
                                                                             cTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                                                                                                                  >  3:807924/1‑73 (MQ=255)
                                                                                aTCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCtca                                                                   <  3:1790496/149‑1 (MQ=255)
                                                                                 tCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                                                                                                                          <  3:141689/77‑1 (MQ=255)
                                                                                 tCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                                                                                                                          >  4:141689/1‑77 (MQ=255)
                                                                                      aGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCg                                                                                  >  1:1168827/1‑128 (MQ=255)
                                                                                      aGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCg                                                                                  <  2:1168827/128‑1 (MQ=255)
                                                                                        ggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                 <  4:1579556/127‑1 (MQ=255)
                                                                                        ggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                 >  3:1579556/1‑127 (MQ=255)
                                                                                                 ccccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                 >  3:380165/1‑118 (MQ=255)
                                                                                                 ccccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                 <  4:380165/118‑1 (MQ=255)
                                                                                                  cccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGt                                                 <  1:1389692/149‑1 (MQ=255)
                                                                                                   ccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGTGATTAAGTTGCTCAGTCAGCGGTTGAATACGt                                                 <  4:1620599/148‑1 (MQ=255)
                                                                                                   ccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGt                                                 >  3:1620599/1‑148 (MQ=255)
                                                                                                    cGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTg                                               >  2:1043223/1‑149 (MQ=255)
                                                                                                              gcggcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                                                              <  3:378305/92‑1 (MQ=255)
                                                                                                              gcggcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                                                              >  4:378305/1‑92 (MQ=255)
                                                                                                              gcggcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTaa                                                                           >  4:92238/1‑111 (MQ=255)
                                                                                                              gcggcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTaa                                                                           <  3:92238/111‑1 (MQ=255)
                                                                                                                 gcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGctc                                  >  3:942984/1‑149 (MQ=255)
                                                                                                                 gcAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGctc                                  <  1:1043223/149‑1 (MQ=255)
                                                                                                                  cAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGctct                                 >  3:1306060/1‑149 (MQ=255)
                                                                                                                   aGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                                                                                                          >  1:505167/1‑75 (MQ=255)
                                                                                                                   aGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCtt                                >  2:737657/1‑149 (MQ=255)
                                                                                                                   aGAAGGAATCAGAATTTCCAAGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                                                                                                          <  2:505167/75‑1 (MQ=255)
                                                                                                                    gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCtca                                                                   <  1:66075/113‑1 (MQ=255)
                                                                                                                    gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCtca                                                                   <  1:314079/113‑1 (MQ=255)
                                                                                                                    gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCtca                                                                   >  2:314079/1‑113 (MQ=255)
                                                                                                                    gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCtca                                                                   >  2:66075/1‑113 (MQ=255)
                                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                                                                                                                   <  1:1504820/64‑1 (MQ=255)
                                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                                                                                                                   >  2:1504820/1‑64 (MQ=255)
                                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCa                              <  4:942984/149‑1 (MQ=255)
                                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCa                              <  4:1306060/149‑1 (MQ=255)
                                                                                                                          aTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGAc                                          >  4:219661/1‑132 (MQ=255)
                                                                                                                          aTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGAc                                          <  3:219661/132‑1 (MQ=255)
                                                                                                                          aTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                             <  4:94482/145‑1 (MQ=255)
                                                                                                                          aTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                             >  3:94482/1‑145 (MQ=255)
                                                                                                                             aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCgaga                                                                               <  3:1241522/92‑1 (MQ=255)
                                                                                                                             aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCgaga                                                                               >  4:1241522/1‑92 (MQ=255)
                                                                                                                             aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGa                                       <  3:708912/132‑1 (MQ=255)
                                                                                                                             aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGa                                       >  4:708912/1‑132 (MQ=255)
                                                                                                                              gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGt                     >  4:240747/1‑149 (MQ=255)
                                                                                                                               aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTa                    <  3:240747/149‑1 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATa                                                    >  1:543286/1‑115 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATa                                                    <  2:543286/115‑1 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGtt                                                >  1:237487/1‑119 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGtt                                                <  2:237487/119‑1 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGc                  >  3:961766/1‑149 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGc                  >  3:1253449/1‑149 (MQ=255)
                                                                                                                                 tttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGc                  >  3:889084/1‑149 (MQ=255)
                                                                                                                                   tCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCgt                >  3:1572966/1‑149 (MQ=255)
                                                                                                                                    ccAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCg                                                            <  4:771968/104‑1 (MQ=255)
                                                                                                                                    ccAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCg                                                            >  3:771968/1‑104 (MQ=255)
                                                                                                                                      aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGtt             >  3:700727/1‑149 (MQ=255)
                                                                                                                                      aGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGtt             >  3:1516392/1‑149 (MQ=255)
                                                                                                                                       ggTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACCCTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAtgc                           <  2:742073/134‑1 (MQ=255)
                                                                                                                                       ggTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAtgc                           >  1:742073/1‑134 (MQ=255)
                                                                                                                                        gTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTc           >  1:713545/1‑149 (MQ=255)
                                                                                                                                          ttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGa                                       >  2:1802436/4‑119 (MQ=255)
                                                                                                                                          ttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGa                                       <  1:1802436/116‑1 (MQ=255)
                                                                                                                                          cAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTg         <  4:1572966/149‑1 (MQ=255)
                                                                                                                                          cAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTg         <  4:961766/149‑1 (MQ=255)
                                                                                                                                           aGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                                                                                                                     >  4:642948/1‑40 (MQ=255)
                                                                                                                                           aGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                                                                                                                     <  3:642948/40‑1 (MQ=255)
                                                                                                                                            gACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGc       >  3:351138/1‑149 (MQ=255)
                                                                                                                                              cGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTc                               >  3:62873/1‑123 (MQ=255)
                                                                                                                                                 gCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCa  <  1:737657/149‑1 (MQ=255)
                                                                                                                                                 gCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCa  >  4:764943/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCA  >  NC_000913/379947‑380239

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