Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,681 T→A 5.0% V381D (GTT→GAT)  waaA → KDO transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,6810TA 5.0% 539.6 / 10.2 174V381D (GTT→GAT) waaAKDO transferase
Reads supporting (aligned to +/- strand):  ref base T (74/90);  new base A (5/4);  total (79/95)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 180 total aligned reads displayed.

TGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAA  >  NC_000913/3809532‑3809830
                                                                                                                                                     |                                                                                                                                                     
tggtggGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGt                                                                                                                                                       >  2:193207/1‑150 (MQ=255)
 ggtggGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATTCCACTACGCTTGCAAAAGAGGtt                                                                                                                                                      <  1:890830/150‑1 (MQ=255)
  gtggGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGttt                                                                                                                                                     >  2:832895/1‑150 (MQ=255)
      gCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                                                                                                                                 >  1:1033253/1‑150 (MQ=255)
      gCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                                                                                                                                 >  2:459039/1‑150 (MQ=255)
       cATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctct                                                                                                                                                >  2:506777/1‑150 (MQ=255)
          aaTCCGCTGGAAGCTGCCGCCCACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                             >  2:745505/1‑150 (MQ=255)
          aaTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                             >  1:1102380/1‑150 (MQ=255)
          aaTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                             <  1:1192379/150‑1 (MQ=255)
           aTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAc                                                                                                                                            >  2:959384/1‑150 (MQ=255)
           aTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAc                                                                                                                                            >  2:354395/1‑150 (MQ=255)
            tCCGCTGGAAGCTGCCGCACACGCTATTCCGGTCTTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACt                                                                                                                                           <  1:969159/150‑1 (MQ=255)
             ccGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgca                                                                                                                                          <  2:447626/150‑5 (MQ=255)
               gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAc                                                                                                                                        <  2:180197/150‑1 (MQ=255)
                    aaGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGTCg                                                                                                                                   <  1:1161638/150‑1 (MQ=255)
                      gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGcc                                                                                                                                 <  2:134615/150‑1 (MQ=255)
                      gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGcc                                                                                                                                 <  2:783500/150‑1 (MQ=255)
                      gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTCCCTCTTTACTCACCGACGcc                                                                                                                                 <  2:846522/150‑1 (MQ=255)
                           cGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                                                                                                                            >  1:378955/1‑150 (MQ=255)
                            gCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAc                                                                                                                           >  1:326700/1‑150 (MQ=255)
                               cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGt                                                                                                                        >  2:85428/1‑150 (MQ=255)
                                acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTa                                                                                                                       <  2:378955/150‑1 (MQ=255)
                                acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTa                                                                                                                       <  2:669289/150‑1 (MQ=255)
                                 cGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAg                                                                                                                      >  1:191397/1‑150 (MQ=255)
                                     aTTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc                                                                                                                  >  1:1015498/1‑150 (MQ=255)
                                       tCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa                                                                                                                >  2:1421100/1‑150 (MQ=255)
                                        ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                               >  1:1252760/1‑150 (MQ=255)
                                        ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                               >  1:457684/1‑150 (MQ=255)
                                        ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                               >  2:102175/1‑150 (MQ=255)
                                        ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATg                                                                                                              >  1:1301507/1‑150 (MQ=255)
                                         cGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATg                                                                                                              >  1:1292317/1‑150 (MQ=255)
                                           gTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGc                                                                                                            <  1:53420/150‑1 (MQ=255)
                                               tGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTc                                                                                                        >  2:1226864/1‑150 (MQ=255)
                                               tGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTc                                                                                                        >  2:1036784/1‑150 (MQ=255)
                                                   ggggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGc                                                                                                    >  2:1111232/1‑150 (MQ=255)
                                                        cGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTg                                                                                               <  2:177507/150‑1 (MQ=255)
                                                         gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGa                                                                                              <  1:1000279/150‑1 (MQ=255)
                                                           aTACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAg                                                                                            <  2:667413/150‑1 (MQ=255)
                                                           aTACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAg                                                                                            <  2:18210/150‑1 (MQ=255)
                                                            tACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGt                                                                                           <  2:1039979/150‑1 (MQ=255)
                                                              cTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTAc                                                                                         >  1:824710/1‑150 (MQ=255)
                                                              cTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTAc                                                                                         >  2:256358/1‑150 (MQ=255)
                                                               ttttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACt                                                                                        >  1:803755/1‑150 (MQ=255)
                                                                 ttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGt                                                                                      <  2:767536/150‑1 (MQ=255)
                                                                     cTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCa                                                                                  >  2:5379/1‑150 (MQ=255)
                                                                      tttAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCaa                                                                                 >  1:1055337/1‑150 (MQ=255)
                                                                      tttAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCaa                                                                                 >  1:816724/1‑150 (MQ=255)
                                                                        tAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCaaaa                                                                               >  1:229611/1‑150 (MQ=255)
                                                                        tAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCaaaa                                                                               <  1:506777/150‑1 (MQ=255)
                                                                         aaaGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAAc                                                                              <  1:699242/150‑1 (MQ=255)
                                                                           aGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCa                                                                            >  1:782918/1‑150 (MQ=255)
                                                                            gACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa                                                                                                                <  1:164784/113‑1 (MQ=255)
                                                                            gACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa                                                                                                                >  2:164784/1‑113 (MQ=255)
                                                                                tttGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTGACCGACGCCGATTACCTTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCATAACCAGGgcg                                                                       <  2:1138276/150‑1 (MQ=255)
                                                                                   gcgcgcGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                    >  2:820979/1‑150 (MQ=255)
                                                                                   gcgcgcGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                    >  2:846265/1‑150 (MQ=255)
                                                                                   gcgcgcGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                    <  1:598713/150‑1 (MQ=255)
                                                                                                gCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCt                                                       >  1:462370/1‑150 (MQ=255)
                                                                                                gCAGGCAAGCGGGCTGATTACCCTTACCGATGCCACTACGCTTGCAACAGAGGTGTCCTCTTTACTCCCCGACGCCGAGTACCGTACTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCt                                                       <  2:598009/150‑1 (MQ=255)
                                                                                                   ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAGACCAGGGCGCGCTACAGCGTCTGCTTCa                                                    >  1:952511/1‑150 (MQ=255)
                                                                                                   ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCa                                                    >  1:46769/1‑150 (MQ=255)
                                                                                                   ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCa                                                    >  1:139523/1‑150 (MQ=255)
                                                                                                    gCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCaa                                                   <  1:256358/150‑1 (MQ=255)
                                                                                                       ctccggCTGATTACCAAAAACGATGCCACTACGCTTGCAAAACAGGTTTCCTCTTAACTCACCGACCCCGAATACCCTAGTTTCTATGGCCGTAATGCCGTTGAAGGACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActg                                                >  2:1267785/5‑150 (MQ=255)
                                                                                                        gCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgc                                               <  2:1403479/150‑1 (MQ=255)
                                                                                                         cGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgct                                              <  1:354395/150‑1 (MQ=255)
                                                                                                          gggCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgctg                                             <  2:1292317/150‑1 (MQ=255)
                                                                                                              tGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATAAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAAc                                         >  1:590073/1‑150 (MQ=255)
                                                                                                                   accgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGCACCTTAc                                    >  1:1080801/1‑150 (MQ=255)
                                                                                                                    ccgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAcc                                   >  2:763281/1‑150 (MQ=255)
                                                                                                                    ccgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTAACAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAcc                                   <  1:1186053/150‑1 (MQ=255)
                                                                                                                      gttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTg                                 <  1:786869/150‑1 (MQ=255)
                                                                                                                          ccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCAc                             >  2:149801/1‑150 (MQ=255)
                                                                                                                           cgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCAcc                            <  1:1111232/150‑1 (MQ=255)
                                                                                                                            gATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAGCCTTACCTGCCACCg                           >  2:316435/1‑150 (MQ=255)
                                                                                                                            gATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCg                           >  2:584917/1‑150 (MQ=255)
                                                                                                                            gATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCg                           >  1:194679/1‑150 (MQ=255)
                                                                                                                             aTGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGa                          <  2:409878/150‑1 (MQ=255)
                                                                                                                                   cTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGc                    <  1:97859/150‑1 (MQ=255)
                                                                                                                                    tACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCa                   <  1:85428/150‑1 (MQ=255)
                                                                                                                                      cGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCAcc                            <  1:750636/139‑1 (MQ=255)
                                                                                                                                      cGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCAtt                 <  1:960949/150‑1 (MQ=255)
                                                                                                                                        gtgagtaaaGAGGAAACCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa               >  2:749738/7‑150 (MQ=255)
                                                                                                                                        cTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCTCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa               <  1:1039789/150‑1 (MQ=255)
                                                                                                                                        cTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa               <  1:871810/150‑1 (MQ=255)
                                                                                                                                        cTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa               >  1:393452/1‑150 (MQ=255)
                                                                                                                                         tgagtaaaGAGGAAACCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAg              >  1:1087561/6‑150 (MQ=255)
                                                                                                                                          gagtaaaGAGGAAACCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAgg             >  2:92223/5‑150 (MQ=255)
                                                                                                                                            cAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGc                                                        >  2:725307/1‑105 (MQ=255)
                                                                                                                                            cAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGc                                                        <  1:725307/105‑1 (MQ=255)
                                                                                                                                             aaaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttg          <  1:292521/150‑1 (MQ=255)
                                                                                                                                             aaaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttg          >  2:213875/1‑150 (MQ=255)
                                                                                                                                             aaaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttg          <  2:1103511/150‑1 (MQ=255)
                                                                                                                                              aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttgt         <  1:1226864/150‑1 (MQ=255)
                                                                                                                                                agagGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttgtt        <  2:388745/149‑1 (MQ=255)
                                                                                                                                                  agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATg     <  2:462370/150‑1 (MQ=255)
                                                                                                                                                  agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATg     <  1:1421100/150‑1 (MQ=255)
                                                                                                                                                   ggTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGc    >  1:111357/1‑150 (MQ=255)
                                                                                                                                                     tttCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCaa  <  2:194679/150‑1 (MQ=255)
                                                                                                                                                     tttCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGGTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCaa  <  1:745505/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAA  >  NC_000913/3809532‑3809830

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: