Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 667,899 C→A 100% D781Y (GAT→TAT)  leuS ← leucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273667,8990CA97.4% 600.1 / 6.2 193D781Y (GAT→TAT) leuSleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (74/114);  total (77/116)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 211 total aligned reads displayed.

GCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCAT  >  CP009273/667754‑668043
                                                                                                                                                 |                                                                                                                                                
gCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagc                                                                                                                                                <  1:915209/148‑1 (MQ=255)
 cACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcc                                                                                                                                               <  1:1132938/148‑1 (MQ=255)
   cGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcctt                                                                                                                                             >  2:1418813/1‑148 (MQ=255)
   cGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcctt                                                                                                                                             >  1:482760/1‑148 (MQ=255)
   cGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcctt                                                                                                                                             <  2:566605/148‑1 (MQ=255)
    gTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttc                                                                                                                                            <  2:2857425/148‑1 (MQ=255)
          gAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCttt                                                                                                                                      <  1:602803/148‑1 (MQ=255)
            aCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCa                                                                                                                                    >  2:2740162/1‑148 (MQ=255)
            aCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCa                                                                                                                                    >  2:2688048/1‑148 (MQ=255)
            aCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCa                                                                                                                                    >  2:1699396/1‑148 (MQ=255)
             ccTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAg                                                                                                                                   <  2:2752662/148‑1 (MQ=255)
             ccTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAg                                                                                                                                   >  2:2841709/1‑148 (MQ=255)
              cTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGt                                                                                                                                  >  1:1940436/1‑148 (MQ=255)
               tGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGtt                                                                                                                                 >  1:1579520/1‑148 (MQ=255)
                     tcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGcc                                                                                                                           <  2:2844159/148‑1 (MQ=255)
                       cGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCcac                                                                                                                         <  2:42771/148‑1 (MQ=255)
                       cGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCcac                                                                                                                         <  2:1822059/148‑1 (MQ=255)
                          tGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGc                                                                                                                      >  2:1239570/1‑148 (MQ=255)
                             gTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTg                                                                                                                   <  1:2981428/148‑1 (MQ=255)
                              tCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGa                                                                                                                  >  1:1698753/1‑148 (MQ=255)
                              tCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGa                                                                                                                  >  2:3035443/1‑148 (MQ=255)
                                cACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAg                                                                                                                <  2:526295/148‑1 (MQ=255)
                                cACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAg                                                                                                                <  2:2018499/148‑1 (MQ=255)
                                   cGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAg                                                                                                             <  1:2252543/148‑1 (MQ=255)
                                   cGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAg                                                                                                             >  1:1465954/1‑148 (MQ=255)
                                      aaCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtg                                                                                                          <  2:303415/148‑1 (MQ=255)
                                      aaCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtg                                                                                                          <  1:394001/148‑1 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                         <  2:2664491/148‑1 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                         >  2:3315878/1‑148 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                         <  2:2433309/148‑1 (MQ=255)
                                        cGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgtg                                                                                                        <  2:1063722/148‑1 (MQ=255)
                                        cGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgtg                                                                                                        <  2:727668/148‑1 (MQ=255)
                                        cGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgtg                                                                                                        >  1:1323245/1‑148 (MQ=255)
                                         ggTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGc                                                                                                       <  2:2198708/148‑1 (MQ=255)
                                            gATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCggg                                                                                                    <  2:1452354/148‑1 (MQ=255)
                                            gATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCggg                                                                                                    >  2:224946/1‑148 (MQ=255)
                                             aTTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCgggg                                                                                                   >  2:627452/1‑148 (MQ=255)
                                                 tGGCACGGACTTGACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGaa                                                                                               <  2:2474853/148‑1 (MQ=255)
                                                  ggCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAAc                                                                                              <  1:1498946/148‑1 (MQ=255)
                                                  ggCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAAc                                                                                              <  2:2806517/148‑1 (MQ=255)
                                                   gCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACg                                                                                             >  1:2137526/1‑148 (MQ=255)
                                                   gCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACg                                                                                             >  1:1259444/1‑148 (MQ=255)
                                                   gCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACg                                                                                             >  1:3008008/1‑148 (MQ=255)
                                                    cACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACgg                                                                                            <  1:2448859/148‑1 (MQ=255)
                                                    cACGGACTTTACCGTGAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACgg                                                                                            <  2:2288412/148‑1 (MQ=255)
                                                        gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTa                                                                                        <  1:253586/148‑1 (MQ=255)
                                                        gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTa                                                                                        <  2:3280391/148‑1 (MQ=255)
                                                        gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTa                                                                                        <  1:3248436/148‑1 (MQ=255)
                                                        gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTa                                                                                        <  2:3280402/148‑1 (MQ=255)
                                                              aCCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATa                                                                                  <  2:330522/148‑1 (MQ=255)
                                                               ccGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATAc                                                                                 >  1:1585491/1‑148 (MQ=255)
                                                               ccGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATAc                                                                                 >  1:3389487/1‑148 (MQ=255)
                                                                  ttAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACCGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                                                                              >  1:1005354/1‑148 (MQ=255)
                                                                  ttAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                                                                              >  1:290461/1‑148 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                            <  2:1543887/148‑1 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                            <  1:3415248/148‑1 (MQ=255)
                                                                            ccacGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcag                                                                    >  2:374525/1‑148 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGTATACGGACAACGGCcagcag                                                                 <  1:3041163/148‑1 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                 >  2:3421877/1‑148 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                 <  2:2114385/148‑1 (MQ=255)
                                                                                gACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGt                                                                <  1:1213741/148‑1 (MQ=255)
                                                                                 aCCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTg                                                               <  1:2692744/148‑1 (MQ=255)
                                                                                  ccAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGc                                                              <  2:937395/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                           <  1:3323962/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                           <  2:1931176/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                           <  2:3348874/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                           <  1:2872020/148‑1 (MQ=255)
                                                                                            gTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATc                                                    >  1:632050/1‑148 (MQ=255)
                                                                                             tCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCa                                                   <  2:1469452/148‑1 (MQ=255)
                                                                                              cTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCag                                                  <  2:1640652/148‑1 (MQ=255)
                                                                                               ttCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCaga                                                 <  2:496126/148‑1 (MQ=255)
                                                                                                tCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCagag                                                <  1:3489955/148‑1 (MQ=255)
                                                                                                tCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCagag                                                <  1:3202131/148‑1 (MQ=255)
                                                                                                  caccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAgcg                                              <  1:858859/148‑1 (MQ=255)
                                                                                                     caTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGa                                           >  2:301007/1‑148 (MQ=255)
                                                                                                      aTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGAt                                          <  2:2488645/148‑1 (MQ=255)
                                                                                                        cGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATcc                                        <  2:1906589/148‑1 (MQ=255)
                                                                                                           tttttCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGc                                     >  1:1826496/1‑148 (MQ=255)
                                                                                                            ttttCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCt                                    <  1:989204/148‑1 (MQ=255)
                                                                                                                cGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGcc                                <  2:1495722/148‑1 (MQ=255)
                                                                                                                  tCAGCAACCGGCCACTGCGCGTTGTAGATATAGCGTTCGCCTTTCATTTACTGCCACAGCGTTAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCAt                              <  1:749004/148‑1 (MQ=255)
                                                                                                                  tCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCAt                              <  1:1418813/148‑1 (MQ=255)
                                                                                                                   cAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATc                             >  2:2372159/1‑148 (MQ=255)
                                                                                                                   cAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATc                             <  2:825762/148‑1 (MQ=255)
                                                                                                                         ccGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTgg                       <  2:857584/148‑1 (MQ=255)
                                                                                                                         ccGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTgg                       <  1:1789977/148‑1 (MQ=255)
                                                                                                                         ccGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTgg                       <  2:2429020/148‑1 (MQ=255)
                                                                                                                            gCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGc                    <  1:1816481/148‑1 (MQ=255)
                                                                                                                            gCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGc                    <  1:789194/148‑1 (MQ=255)
                                                                                                                              cACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCtt                  <  1:1992831/148‑1 (MQ=255)
                                                                                                                                cGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGAGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTc                >  1:419722/1‑148 (MQ=255)
                                                                                                                                  gcgcGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGc              >  1:1911882/1‑148 (MQ=255)
                                                                                                                                  gcgcGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAAAGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGc              >  1:3196328/1‑148 (MQ=255)
                                                                                                                                     cgTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAg           >  1:725529/1‑148 (MQ=255)
                                                                                                                                      gTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGt          <  1:248182/148‑1 (MQ=255)
                                                                                                                                      gTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGt          <  2:1971651/148‑1 (MQ=255)
                                                                                                                                       ttGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGtt         >  2:3201227/1‑148 (MQ=255)
                                                                                                                                         gTCGATATAGCCTTCGCCTGTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTg       <  2:2970459/148‑1 (MQ=255)
                                                                                                                                             atatAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtca   <  2:3456247/148‑1 (MQ=255)
                                                                                                                                              tataGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtcat  <  2:90185/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCAT  >  CP009273/667754‑668043

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: