Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,535,948 G→C 100% D39H (GAT→CAT)  crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,535,9480GC99.7% 1197.7 / ‑5.2 319D39H (GAT→CAT) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (147/171);  total (148/171)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 319 total aligned reads displayed.

ATGCGGCCCAATTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTT  >  NC_000913/2535802‑2536091
                                                                                                                                                  |                                                                                                                                               
     gCCCAATTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGt                                                                                                                                          <  3:211483/149‑1 (MQ=255)
       ccAATTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGttt                                                                                                                                        >  3:2202890/1‑149 (MQ=255)
         aaTTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGttttt                                                                                                                                      <  4:1987989/149‑1 (MQ=255)
         aaTTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGttttt                                                                                                                                      <  3:8755/149‑1 (MQ=255)
           tttACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGc                                                                                                                                    <  1:2010822/149‑1 (MQ=255)
             tACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTCCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCACTCTCTGGCGAGATCGTCAATATCGAAGACGTACCGCATGTCGTCTTTGCgg                                                                                                                                  <  4:2075803/149‑1 (MQ=255)
              aCTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGttt                                                                                                                                        >  1:1633062/1‑142 (MQ=255)
                tGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGTTCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGaaa                                                                                                                               >  2:882307/1‑149 (MQ=255)
                   ttAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGaaaaaa                                                                                                                            <  2:2316234/149‑1 (MQ=255)
                   ttAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGaaaaaa                                                                                                                            >  1:663999/1‑149 (MQ=255)
                       gagaAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGt                                                                                                                        >  3:403467/1‑149 (MQ=255)
                        agaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                       <  3:2039935/149‑1 (MQ=255)
                        agaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                       <  1:1755143/149‑1 (MQ=255)
                        agaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                       <  1:1821996/149‑1 (MQ=255)
                        agaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                       <  1:463491/149‑1 (MQ=255)
                         gaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTg                                                                                                                      >  4:1088894/1‑149 (MQ=255)
                         gaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTg                                                                                                                      >  1:1265311/1‑149 (MQ=255)
                            gaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCAGCATGTCGTTTTTGCGGAAAAAATCGTTGGTg                                                                                                                   >  4:1466141/1‑149 (MQ=255)
                              tCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGAt                                                                                                                 >  2:2011356/1‑149 (MQ=255)
                               cATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGGAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                >  2:1893495/1‑149 (MQ=255)
                               cATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                <  1:762602/149‑1 (MQ=255)
                               cATGGGTTTGTTCGATAAACTGAAATCTCTGGCTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                <  2:293996/149‑1 (MQ=255)
                                aTGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATgg                                                                                                               >  1:963325/1‑149 (MQ=255)
                                 tGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGt                                                                                                              >  1:2099678/1‑149 (MQ=255)
                                 tGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGATGGTGATgg                                                                                                               >  4:528964/1‑148 (MQ=255)
                                        gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTa                                                                                                       >  4:1750718/1‑149 (MQ=255)
                                         ttCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTAt                                                                                                      <  4:619506/149‑1 (MQ=255)
                                         ttCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTAt                                                                                                      >  3:1109208/1‑149 (MQ=255)
                                          tCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATc                                                                                                     <  1:2490398/149‑1 (MQ=255)
                                            gATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCaa                                                                                                   >  3:764096/1‑149 (MQ=255)
                                             aTAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCaaa                                                                                                  <  1:1077367/149‑1 (MQ=255)
                                                  cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                             <  4:159449/149‑1 (MQ=255)
                                                     aaaTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACggg                                                                                          >  2:384228/1‑149 (MQ=255)
                                                      aaTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGt                                                                                         <  3:612657/149‑1 (MQ=255)
                                                       aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                        <  4:2127650/149‑1 (MQ=255)
                                                         ctctGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGa                                                                                                                  >  3:1313604/1‑121 (MQ=255)
                                                          tctGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACa                                                                                     >  1:2158144/1‑149 (MQ=255)
                                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                 <  2:1826762/149‑1 (MQ=255)
                                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                 <  2:1377654/149‑1 (MQ=255)
                                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                 <  3:1577191/149‑1 (MQ=255)
                                                               tttCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATg                                                                                <  4:1790258/149‑1 (MQ=255)
                                                                 tCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATGTCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGt                                                                              <  4:1289901/149‑1 (MQ=255)
                                                                  ccGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                             <  4:1508007/149‑1 (MQ=255)
                                                                  ccGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                             <  3:1855372/149‑1 (MQ=255)
                                                                    gacgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATTGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTTCTATCAAACCAACGGGTAACAAAATGGTcgc                                                                           <  2:2189461/149‑1 (MQ=255)
                                                                    gacgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTcgc                                                                           >  3:2325697/1‑149 (MQ=255)
                                                                       gacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                 <  3:670585/140‑1 (MQ=255)
                                                                          aagaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAGTGGTCGCGCCAGt                                                                     >  2:921586/1‑149 (MQ=255)
                                                                           agaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTa                                                                    >  3:1440896/1‑149 (MQ=255)
                                                                           agaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTa                                                                    >  2:2145922/1‑149 (MQ=255)
                                                                            gaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTAt                                                                                                      >  1:2280727/1‑114 (MQ=255)
                                                                            gaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAg                                                                   >  3:2237266/1‑149 (MQ=255)
                                                                              agGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGAc                                                                 <  1:2432512/149‑1 (MQ=255)
                                                                               ggATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGcc                                                                        >  2:1589103/1‑141 (MQ=255)
                                                                                    ccGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCAcc                                                           <  3:494280/149‑1 (MQ=255)
                                                                                     cGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCa                                                          >  1:1155211/1‑149 (MQ=255)
                                                                                     cGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCa                                                          >  1:924205/1‑149 (MQ=255)
                                                                                      ggAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCAt                                                         <  4:747423/149‑1 (MQ=255)
                                                                                       gAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCAtt                                                        <  2:1502349/149‑1 (MQ=255)
                                                                                        aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                       <  1:2393703/149‑1 (MQ=255)
                                                                                        aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                       <  2:2052256/149‑1 (MQ=255)
                                                                                        aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                       <  2:1639715/149‑1 (MQ=255)
                                                                                           tATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCaa                                                                                              >  2:568728/1‑107 (MQ=255)
                                                                                           tATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTa                                                    >  1:1861970/1‑149 (MQ=255)
                                                                                            aTTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTAt                                                                                                            <  1:572635/92‑1 (MQ=255)
                                                                                            aTTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaa                                                   >  2:2382225/1‑149 (MQ=255)
                                                                                             ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGt                                                                                                                        <  3:1053000/79‑1 (MQ=255)
                                                                                             ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCTAGTAGACGGCACCATTGGTaaa                                                  >  4:2015790/1‑149 (MQ=255)
                                                                                             ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaaa                                                  >  4:2383818/1‑149 (MQ=255)
                                                                                             ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaaa                                                  <  3:1742973/149‑1 (MQ=255)
                                                                                               gagaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAAt                                                <  2:1699921/149‑1 (MQ=255)
                                                                                               gagaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAAt                                                >  1:1699921/1‑149 (MQ=255)
                                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                             >  3:1657921/1‑149 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGCGATCGTCAATATCGAACACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACCAAATGTTCGCGCCAGTAGACGGCACCATTGGTAATATCttt                                            <  2:254644/149‑1 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            <  1:431420/149‑1 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            <  1:1941097/149‑1 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            <  4:1657921/149‑1 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            <  4:1348999/149‑1 (MQ=255)
                                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGAAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            >  3:339129/1‑149 (MQ=255)
                                                                                                         cTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                             >  4:560188/1‑142 (MQ=255)
                                                                                                         cTCCGCTCTCTGGCGAGATCGTCAATAACGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                             <  3:560188/142‑1 (MQ=255)
                                                                                                            cGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaac                                   <  4:2202890/149‑1 (MQ=255)
                                                                                                                 tctGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGc                              >  2:1125424/1‑149 (MQ=255)
                                                                                                                 tctGGCGAGATCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGc                              >  3:746097/1‑149 (MQ=255)
                                                                                                                         gaTCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGG‑AAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                     <  1:1499107/149‑1 (MQ=255)
                                                                                                                          aTCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACGAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                     >  3:2390117/1‑149 (MQ=255)
                                                                                                                          aTCGTCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCa                                                          <  4:704703/112‑1 (MQ=255)
                                                                                                                              tCAATATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGaa                 >  2:68271/1‑149 (MQ=255)
                                                                                                                                aaTATCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATc               <  3:2010355/149‑1 (MQ=255)
                                                                                                                                 atatCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCt              >  1:807683/1‑149 (MQ=255)
                                                                                                                                  tatCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTg             >  2:512228/1‑149 (MQ=255)
                                                                                                                                    tCGAAGTCGTGCCGCATGTCGTGTTTGCGGAAAAAATCGTTGGTAATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCCCCATTGGTAAACTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGAt           <  4:1892966/149‑1 (MQ=255)
                                                                                                                                    tCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGAt           <  1:1043342/149‑1 (MQ=255)
                                                                                                                                    tCGAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGAt           >  1:1297186/1‑149 (MQ=255)
                                                                                                                                      gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAg         >  3:2135686/1‑149 (MQ=255)
                                                                                                                                      gAAGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                            >  2:1167409/1‑114 (MQ=255)
                                                                                                                                       aaGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgc        >  3:297526/1‑149 (MQ=255)
                                                                                                                                       aaGACGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgc        >  2:1344292/1‑149 (MQ=255)
                                                                                                                                          aCGTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcggc     <  4:2237266/149‑1 (MQ=255)
                                                                                                                                            gTGCCGCATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGt   >  4:1045371/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATGCGGCCCAATTTACTGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTT  >  NC_000913/2535802‑2536091

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