Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,783,857 A→T 6.2% I283I (ATT→ATA tdh ← L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,783,8570AT 6.2% 604.6 / 23.8 192I283I (ATT→ATAtdhL‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (101/79);  new base T (6/6);  total (107/85)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 201 total aligned reads displayed.

ACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACAT  >  CP009273/3783709‑3784005
                                                                                                                                                    |                                                                                                                                                    
acgggCCCGAACGCATAGCGTCAAACCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTa                                                                                                                                                      <  1:2278075/146‑1 (MQ=255)
aCTGGCCCGAACGCATCGCGTCAAAGCCCTTCTGTAAATCATCGATAGAGAAACTATGGGTAATGATCGGCTAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATATTGTACCAGGTTTAAAACATCTCACGACCGTAAATACCTTTa                                                                                                                                                      <  1:1866272/149‑1 (MQ=255)
 cgggCCCGAACGCATCGCTTCAAACCCCTTCTGGAAATCATCTATAGATAAACGTTGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTaa                                                                                                                                                     <  1:2104289/147‑1 (MQ=255)
 cTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTaa                                                                                                                                                     >  1:1107377/1‑149 (MQ=255)
   ggCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATg                                                                                                                                                   >  1:1228552/1‑149 (MQ=255)
    gCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGa                                                                                                                                                  <  2:3719205/149‑1 (MQ=255)
    ccccGAACGAATCGCTACAACGCCATTCTGTAAATAATCTATACATAAACGTTGGGTAATGATCGGCTAAAGATCGTGGCCAGAATGAATCAGCGCCCCCATCTTGTACCCGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGa                                                                                                                                                  <  1:1752083/148‑1 (MQ=255)
       ctaaCGCATAGCGTCAAAGCCCTTCTGCAAATCATCGATAGAGAAACTTTGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACa                                                                                                                                               <  1:2528474/147‑1 (MQ=255)
       cGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACa                                                                                                                                               <  2:765522/149‑1 (MQ=255)
       cGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACa                                                                                                                                               >  1:1882394/1‑149 (MQ=255)
         aaCGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAg                                                                                                                                             <  2:2834561/149‑1 (MQ=255)
            gCATAGCGTCAAAGCCCTTCTGGAAATAATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCt                                                                                                                                          <  2:1357534/149‑1 (MQ=255)
              aTAGCGTCAAAGCCCTTCTGGAGATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCttt                                                                                                                                        >  1:2582866/1‑149 (MQ=255)
                aGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTaa                                                                                                                                      <  2:2405826/149‑1 (MQ=255)
                aGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTaa                                                                                                                                      >  2:3157030/1‑149 (MQ=255)
                aGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTaa                                                                                                                                      >  1:3013752/1‑149 (MQ=255)
                  cGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAg                                                                                                                                    <  2:1583986/149‑1 (MQ=255)
                    tCAACGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGAt                                                                                                                                  <  2:1076126/149‑1 (MQ=255)
                    tCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGAt                                                                                                                                  >  1:3957555/1‑149 (MQ=255)
                     cAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATc                                                                                                                                 <  2:677571/149‑1 (MQ=255)
                     cAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATc                                                                                                                                 >  1:1086766/1‑149 (MQ=255)
                     cAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATc                                                                                                                                 >  1:3210793/1‑149 (MQ=255)
                      aaaGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCa                                                                                                                                <  1:2708626/149‑1 (MQ=255)
                       aaGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCAc                                                                                                                               >  2:4058428/1‑149 (MQ=255)
                       aaGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCAc                                                                                                                               >  1:2842184/1‑149 (MQ=255)
                         gCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACtt                                                                                                                             <  2:2901448/149‑1 (MQ=255)
                          cccTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACttt                                                                                                                            <  2:1423005/149‑1 (MQ=255)
                          cccTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACttt                                                                                                                            <  2:1979114/149‑1 (MQ=255)
                            cTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTgg                                                                                                                          <  2:880104/149‑1 (MQ=255)
                             ttCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGt                                                                                                                         >  1:992598/1‑149 (MQ=255)
                                        atcGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAg                                                                                                              >  2:770048/1‑149 (MQ=255)
                                        atcGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAg                                                                                                              >  2:2120972/1‑149 (MQ=255)
                                          cGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGAc                                                                                                            >  2:1659879/1‑149 (MQ=255)
                                          cGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGAc                                                                                                            <  2:1294646/149‑1 (MQ=255)
                                          cGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGAc                                                                                                            <  2:3545430/149‑1 (MQ=255)
                                           gATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACa                                                                                                           <  2:3814734/149‑1 (MQ=255)
                                           gATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACa                                                                                                           <  2:806440/149‑1 (MQ=255)
                                            aTAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACat                                                                                                          >  2:76592/1‑149 (MQ=255)
                                              agagAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACatat                                                                                                        >  1:3509154/1‑149 (MQ=255)
                                              agagAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACatat                                                                                                        <  1:2164072/149‑1 (MQ=255)
                                               gagaAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATc                                                                                                       >  2:2010660/1‑149 (MQ=255)
                                               gagaAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATc                                                                                                       >  2:113933/1‑149 (MQ=255)
                                                 gaAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACt                                                                                                                              >  1:761151/1‑124 (MQ=255)
                                                 gaAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAg                                                                                                     >  2:2900403/1‑149 (MQ=255)
                                                   aaCGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTCAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGAc                                                                                                   >  2:3876011/1‑149 (MQ=255)
                                                    aCGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGAcg                                                                                                  >  2:2631094/1‑149 (MQ=255)
                                                    aCGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGAcg                                                                                                  >  2:913222/1‑149 (MQ=255)
                                                         gggTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGa                                                                                             >  2:939861/1‑149 (MQ=255)
                                                         gggTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGa                                                                                             >  2:1253857/1‑149 (MQ=255)
                                                          ggTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGaa                                                                                            <  2:2577075/149‑1 (MQ=255)
                                                           gTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTTATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAAt                                                                                           >  1:2345213/1‑149 (MQ=255)
                                                           gTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGact                                                                                           >  2:3740600/1‑147 (MQ=255)
                                                            tAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATa                                                                                          <  2:3270755/149‑1 (MQ=255)
                                                                gATCGGCTAAATATCGAGGCCATACTGAATCAGCGCCGCAATCTTGTACCAGGTTACAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCa                                                                                      <  1:723264/149‑1 (MQ=255)
                                                                gATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCa                                                                                      <  2:3201212/149‑1 (MQ=255)
                                                                gATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCa                                                                                      <  1:723205/149‑1 (MQ=255)
                                                                  tCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGc                                                                                    >  1:433665/1‑149 (MQ=255)
                                                                        aaaGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCa                                                                              >  1:2213946/1‑149 (MQ=255)
                                                                         aaGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCaa                                                                             >  2:2515990/1‑149 (MQ=255)
                                                                                aGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGgcc                                                                      >  1:777202/1‑149 (MQ=255)
                                                                                 ggCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGgccg                                                                     >  1:3861246/1‑149 (MQ=255)
                                                                                     aGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACACGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGt                                                                 <  2:3001781/149‑1 (MQ=255)
                                                                                      gACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTg                                                                >  2:6177/1‑149 (MQ=255)
                                                                                        cTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGAt                                                              >  2:2363438/1‑149 (MQ=255)
                                                                                          gAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTc                                                            <  2:4084887/149‑1 (MQ=255)
                                                                                               cgcgcCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGgt                                                       <  2:775/148‑1 (MQ=255)
                                                                                                   ccgccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCa                                                   >  1:2899869/1‑149 (MQ=255)
                                                                                                      ccATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGc                                                >  1:605117/1‑149 (MQ=255)
                                                                                                         tCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATg                                             >  2:599100/1‑149 (MQ=255)
                                                                                                           ttGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGt                                           >  1:2381714/1‑149 (MQ=255)
                                                                                                           ttGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGt                                           <  2:153093/149‑1 (MQ=255)
                                                                                                             gTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTAc                                         <  1:1550709/149‑1 (MQ=255)
                                                                                                             gTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTAc                                         >  1:2632870/1‑149 (MQ=255)
                                                                                                              tACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACg                                        >  1:2427057/1‑149 (MQ=255)
                                                                                                                 cAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGaaa                                     <  2:3382312/149‑1 (MQ=255)
                                                                                                                  aGGTTTCAAACATCTTACGACCGTAAATACCTTTTATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAAc                                    <  2:583337/149‑1 (MQ=255)
                                                                                                                       tCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTgg                               <  2:2528013/149‑1 (MQ=255)
                                                                                                                              tctcACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCa                        >  1:3826096/1‑149 (MQ=255)
                                                                                                                                    gACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGAc                  <  1:1789612/149‑1 (MQ=255)
                                                                                                                                     aCCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACa                 <  1:1340524/149‑1 (MQ=255)
                                                                                                                                       cGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACAtt               >  1:2652573/1‑149 (MQ=255)
                                                                                                                                       cGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACAtt               <  2:3901138/149‑1 (MQ=255)
                                                                                                                                         tAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTc             >  1:2342025/1‑149 (MQ=255)
                                                                                                                                         tAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTc             >  2:88001/1‑149 (MQ=255)
                                                                                                                                         tAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTc             <  1:96907/149‑1 (MQ=255)
                                                                                                                                          aaaTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTcc            >  1:502075/1‑149 (MQ=255)
                                                                                                                                          aaaTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTcc            >  1:1456931/1‑149 (MQ=255)
                                                                                                                                           aaTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCa           >  2:3574523/1‑149 (MQ=255)
                                                                                                                                           aaTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCa           >  1:1156880/1‑149 (MQ=255)
                                                                                                                                            aTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAg          >  1:786192/1‑149 (MQ=255)
                                                                                                                                            aTACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAg          >  2:3174612/1‑149 (MQ=255)
                                                                                                                                             tACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGa         >  1:2423624/1‑149 (MQ=255)
                                                                                                                                              aCCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGAc        >  2:3782227/1‑149 (MQ=255)
                                                                                                                                               ccTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGAcc       >  2:3614474/1‑149 (MQ=255)
                                                                                                                                               ccTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGAcc       <  1:867655/149‑1 (MQ=255)
                                                                                                                                                cTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCg      >  1:3888252/1‑149 (MQ=255)
                                                                                                                                                  ttAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGAc    <  2:2842184/149‑1 (MQ=255)
                                                                                                                                                   tAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACa   <  2:2980860/149‑1 (MQ=255)
                                                                                                                                                   tAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACa   <  2:1228552/149‑1 (MQ=255)
                                                                                                                                                    aaTGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACAt  <  2:2124779/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGGCCGCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACAT  >  CP009273/3783709‑3784005

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