Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,682 T→A 7.3% V381V (GTT→GTA) ‡ waaA → KDO transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,6820TA 7.3% 317.6 / 9.0 106V381V (GTT→GTA) ‡waaAKDO transferase
Reads supporting (aligned to +/- strand):  ref base T (44/54);  new base A (2/6);  total (46/60)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 106 total aligned reads displayed.

TGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCA  >  NC_000913/3809535‑3809829
                                                                                                                                                   |                                                                                                                                                   
tggGCATAATCCGCTGGGAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                                                                                                                                   <  1:1468825/150‑1 (MQ=255)
 gggCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTcc                                                                                                                                                  >  2:1342841/1‑150 (MQ=255)
 gggCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTcc                                                                                                                                                  <  2:577679/150‑1 (MQ=255)
  ggCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCct                                                                                                                                                 <  1:467470/150‑1 (MQ=255)
   gCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                                                                                                                                <  2:851540/150‑1 (MQ=255)
   gCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctc                                                                                                                                                <  1:1478801/150‑1 (MQ=255)
    cATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctct                                                                                                                                               <  1:813670/150‑1 (MQ=255)
     aTAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCtt                                                                                                                                              <  2:590527/150‑1 (MQ=255)
     aTAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCtt                                                                                                                                              <  1:1393358/150‑1 (MQ=255)
      tAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTcc                                                                                                                                                  >  2:1535303/1‑145 (MQ=255)
      tAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTcc                                                                                                                                                  <  1:1535303/145‑1 (MQ=255)
      tAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCttt                                                                                                                                             <  1:1069484/150‑1 (MQ=255)
       aaTCCGCTGGAAGTTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                            <  1:46827/150‑1 (MQ=255)
       aaTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                            <  1:1347173/150‑1 (MQ=255)
        aTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttg                                                                                                                                           <  2:195787/150‑3 (MQ=255)
         tCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgc                                                                                                                                          <  1:1528607/150‑4 (MQ=255)
         tCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCGGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgc                                                                                                                                          <  2:1356824/150‑4 (MQ=255)
         tCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACGTGTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACt                                                                                                                                          >  2:1537810/1‑150 (MQ=255)
         cccGCTGCAAGCTGCCGCACACGCTATTCCGGTCTTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgc                                                                                                                                          <  2:464907/149‑4 (MQ=255)
          ccGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTc                                                                                                                                         >  2:1089609/1‑150 (MQ=255)
          ccGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTc                                                                                                                                         <  2:908795/150‑1 (MQ=255)
            gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCGGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgcaag                                                                                                                                       <  2:1384370/150‑7 (MQ=255)
                  aGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGc                                                                                                                                 <  1:909463/150‑1 (MQ=255)
                   gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGcc                                                                                                                                <  2:460972/150‑1 (MQ=255)
                    cTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTGCTCTTTACTCACCGACGCCg                                                                                                                               <  2:1290171/150‑1 (MQ=255)
                        cGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                                                                                                                           >  1:387383/1‑150 (MQ=255)
                        cGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                                                                                                                           >  1:710465/1‑150 (MQ=255)
                        cGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                                                                                                                           <  1:778981/150‑1 (MQ=255)
                           acacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCg                                                                                                                        <  2:723042/150‑1 (MQ=255)
                            cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTCCGCTTGCAAAAGAGGTTTCCTCGTTACTCACCGCCGCCGATTGCCgg                                                                                                                       >  2:298930/1‑149 (MQ=255)
                            cacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGt                                                                                                                       >  2:1116139/1‑150 (MQ=255)
                             acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAc                                                                                                                                           >  2:1361347/1‑129 (MQ=255)
                             acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAc                                                                                                                                           <  1:1361347/129‑1 (MQ=255)
                              cGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAg                                                                                                                     >  2:1336488/1‑150 (MQ=255)
                                  aTTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc                                                                                                                 <  1:32749/150‑1 (MQ=255)
                                  aTTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc                                                                                                                 <  1:935735/150‑1 (MQ=255)
                                    tCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa                                                                                                               >  1:831303/1‑150 (MQ=255)
                                     ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                              <  1:494356/150‑1 (MQ=255)
                                     ccGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                              >  2:1232197/1‑150 (MQ=255)
                                     ccGATATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                              <  2:912097/150‑1 (MQ=255)
                                        gTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGc                                                                                                           <  2:1035470/150‑1 (MQ=255)
                                         tATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGcc                                                                                                          >  1:930771/1‑150 (MQ=255)
                                                ggggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTAATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGc                                                                                                   >  2:1203760/1‑150 (MQ=255)
                                                 gggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGcc                                                                                                  >  1:790465/1‑150 (MQ=255)
                                                      gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGa                                                                                             <  1:1446337/150‑1 (MQ=255)
                                                      gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGa                                                                                             <  1:1165317/150‑1 (MQ=255)
                                                         tACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGt                                                                                          >  2:491891/1‑150 (MQ=255)
                                                         tACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGt                                                                                          <  1:1502257/150‑1 (MQ=255)
                                                           cTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTAc                                                                                        >  1:194923/1‑150 (MQ=255)
                                                           cTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTAc                                                                                        <  1:298930/150‑1 (MQ=255)
                                                              ttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGt                                                                                     <  2:916126/150‑1 (MQ=255)
                                                               tAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTa                                                                                    >  2:554943/1‑150 (MQ=255)
                                                                aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTAt                                                                                   <  1:187528/150‑1 (MQ=255)
                                                                  cTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCa                                                                                 >  2:432241/1‑150 (MQ=255)
                                                                     tAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCaaaa                                                                              <  1:69995/150‑1 (MQ=255)
                                                                      aaaGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGCCGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAAc                                                                             >  2:10990/1‑150 (MQ=255)
                                                                      aaaGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAAc                                                                             >  2:528127/1‑150 (MQ=255)
                                                                        aGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCa                                                                           <  1:1089609/150‑1 (MQ=255)
                                                                         gACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAg                                                                          <  2:930771/150‑1 (MQ=255)
                                                                         gACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAg                                                                          >  2:548896/1‑150 (MQ=255)
                                                                           cttttGCGCGCTTCTGGAGGAGGGAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTAGTCAGCGACGCCGATTACCGTAGTTTGTATGGGGGTCATGCCGTTGATTTACTGTATGAAAACCAggg                                                                        <  2:710465/148‑1 (MQ=255)
                                                                           cATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAggg                                                                        <  2:678572/150‑1 (MQ=255)
                                                                           cATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAggg                                                                        >  1:615239/1‑150 (MQ=255)
                                                                                 cgcgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACACCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCt                                                                  >  2:913294/1‑150 (MQ=255)
                                                                                   cgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTAc                                                                >  1:665674/1‑150 (MQ=255)
                                                                                   cgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTAc                                                                <  1:173414/150‑1 (MQ=255)
                                                                                     cgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAg                                                              >  1:579857/1‑150 (MQ=255)
                                                                                       gCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCg                                                            <  2:790465/150‑1 (MQ=255)
                                                                                         tGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTc                                                          <  1:593138/150‑1 (MQ=255)
                                                                                            aGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGc                                                       <  1:888266/150‑1 (MQ=255)
                                                                                             gCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCt                                                      <  2:579857/150‑1 (MQ=255)
                                                                                              cAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCtt                                                     <  2:86210/150‑1 (MQ=255)
                                                                                                ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCa                                                   >  2:1371247/1‑150 (MQ=255)
                                                                                                ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCa                                                   >  2:116835/1‑150 (MQ=255)
                                                                                                          cTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGaa                                         >  2:995102/1‑150 (MQ=255)
                                                                                                            gATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAAcc                                       >  1:541307/1‑150 (MQ=255)
                                                                                                            gATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAAcc                                       <  1:812798/150‑1 (MQ=255)
                                                                                                                accgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAc                                   <  1:466859/150‑1 (MQ=255)
                                                                                                                accgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAc                                   <  1:461714/150‑1 (MQ=255)
                                                                                                                   gttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTg                                >  2:253988/1‑150 (MQ=255)
                                                                                                                      accgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCa                             >  1:1001217/1‑150 (MQ=255)
                                                                                                                       ccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCAc                            >  1:127324/1‑150 (MQ=255)
                                                                                                                         gATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCg                          >  1:242355/1‑150 (MQ=255)
                                                                                                                         gATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCg                          <  2:935822/150‑1 (MQ=255)
                                                                                                                          aTGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGa                         <  2:1205047/150‑1 (MQ=255)
                                                                                                                              cACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAAc                     <  1:1342841/150‑1 (MQ=255)
                                                                                                                                      ttGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAg             <  1:1203760/150‑1 (MQ=255)
                                                                                                                                       tGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAgg            >  1:1131961/1‑150 (MQ=255)
                                                                                                                                        gCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgt           >  2:1094682/1‑150 (MQ=255)
                                                                                                                                        agtaaaGAGGAAACCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgt           >  1:297421/4‑150 (MQ=255)
                                                                                                                                         gtaaaGAGGAAACCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgtt          >  1:670035/3‑150 (MQ=255)
                                                                                                                                         cAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgtt          >  2:916063/1‑150 (MQ=255)
                                                                                                                                          aaaaGAGGATTCCTCTTTACTGAACGACCCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTGGAAGTACTGTATCAAAACCAGGCCGCGCTACAGAGTCTGCTTCAACTGCTGGACCCTTACCTGCCACCGAAAACGCATTGAGgttt         >  2:597570/1‑149 (MQ=255)
                                                                                                                                           aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttgt        >  2:78127/1‑150 (MQ=255)
                                                                                                                                           aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttgt        >  2:821039/1‑150 (MQ=255)
                                                                                                                                           aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttgt        <  1:253988/150‑1 (MQ=255)
                                                                                                                                             agagGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTa      <  1:1116139/150‑1 (MQ=255)
                                                                                                                                              gagGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTAt     <  1:6026/150‑1 (MQ=255)
                                                                                                                                               agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATTGCCGGCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCGACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAACACGCATTGAGGTTGTTATg    <  2:1533407/150‑1 (MQ=255)
                                                                                                                                               agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATg    <  2:937382/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCA  >  NC_000913/3809535‑3809829

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: