Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,812,983 C→T 6.4% Q796Q (CAG→CAA)  alaS ← alanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,812,9830CT 6.4% 1154.5 / 44.5 329Q796Q (CAG→CAA) alaSalanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (134/172);  new base T (10/11);  total (146/183)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 334 total aligned reads displayed.

CGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCTTTAAGCTGTTG  >  CP009273/2812834‑2813132
                                                                                                                                                     |                                                                                                                                                     
       caccCTTGCCGGCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCAAGCACGATAATTTTCTACCCCAGCTGAtttt                                                                                                                                                <  1:1440901/150‑1 (MQ=255)
         cccTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGAtttttt                                                                                                                                              <  1:329995/150‑1 (MQ=255)
          ccTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTa                                                                                                                                             >  2:2040971/1‑150 (MQ=255)
          ccTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTa                                                                                                                                             <  1:4045304/150‑1 (MQ=255)
           cTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTaa                                                                                                                                            <  2:906910/150‑1 (MQ=255)
           cTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGGTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTaa                                                                                                                                            <  1:2336589/150‑1 (MQ=255)
             tGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAAt                                                                                                                                          >  2:3584033/1‑150 (MQ=255)
              gccgccCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATc                                                                                                                                         <  1:3381606/150‑1 (MQ=255)
               ccgccCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCg                                                                                                                                        <  2:4749528/150‑1 (MQ=255)
                  cccACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCa                                                                                                                                     <  2:2306651/150‑1 (MQ=255)
                  cccACCAGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCa                                                                                                                                     <  2:2743993/150‑1 (MQ=255)
                   ccACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCaa                                                                                                                                    <  2:3404295/150‑1 (MQ=255)
                       ctgctgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCAt                                                                                                                                <  1:3673338/150‑1 (MQ=255)
                        tgctgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATg                                                                                                                               <  2:781087/150‑1 (MQ=255)
                        tgctgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATg                                                                                                                               <  2:2277658/150‑1 (MQ=255)
                         gctgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATgg                                                                                                                              >  2:810282/1‑150 (MQ=255)
                         gctgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATgg                                                                                                                              <  2:311476/150‑1 (MQ=255)
                           tgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTa                                                                                                                            <  1:2084398/150‑1 (MQ=255)
                           tgAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTa                                                                                                                            <  1:4687690/150‑1 (MQ=255)
                             aGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACg                                                                                                                          >  2:1306114/1‑150 (MQ=255)
                               cGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCAACCATGGTACGCa                                                                                                                        <  1:978553/150‑1 (MQ=255)
                               cGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCa                                                                                                                        <  2:3758683/150‑1 (MQ=255)
                                gACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCaa                                                                                                                       <  2:2593167/150‑1 (MQ=255)
                                  ccATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCAACCATGGTACGCAACa                                                                                                                     >  1:4655409/1‑150 (MQ=255)
                                  ccATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACa                                                                                                                     >  1:503478/1‑150 (MQ=255)
                                   cATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACAt                                                                                                                    >  2:4892642/1‑150 (MQ=255)
                                    aTACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCAACCATGGTACGCAACAtt                                                                                                                   <  1:2027608/150‑1 (MQ=255)
                                    aTACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACAtt                                                                                                                   <  2:3927842/150‑1 (MQ=255)
                                      aCCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACAtttt                                                                                                                 <  1:1606408/150‑1 (MQ=255)
                                       ccACTCAGTTCCCCTTCTTTCACACGATCTGTGACGTACTTTGATACGCCTGCAATCAGAGTTACCTTACCTTAGACTTGCTTTGCCAGCACGATATTTGTCGACCCCAGCTGATTTTTTATTTCGTGATCCATGGGACTCATAATTTTc                                                                                                                <  2:2323917/150‑1 (MQ=255)
                                          aTCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGc                                                                                                             >  1:3752407/1‑150 (MQ=255)
                                           tCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCt                                                                                                            <  2:3688158/150‑1 (MQ=255)
                                                      gCTTTCACACGATCTGAGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTa                                                                                                 <  2:545293/150‑1 (MQ=255)
                                                          tCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCt                                                                                             <  1:914643/150‑1 (MQ=255)
                                                            acacGATCTGTGACGTCCTTAGATACGCCTGCAATCAGCGAACCCTTACCTTCGACTACCGTTGCCAGCACGATAATGGTCGACCCCAGCTGATTTTTTAAATCGTCAACCCTGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCg                                                                                           <  2:4316862/150‑1 (MQ=255)
                                                            acacGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCg                                                                                           >  2:3352121/1‑150 (MQ=255)
                                                             cacGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGATCAACACCGCTAAGCTCGc                                                                                          >  2:1383969/1‑150 (MQ=255)
                                                               cGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTGACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTa                                                                                        <  1:2811225/150‑1 (MQ=255)
                                                                 aTCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAAc                                                                                      >  1:1113957/1‑150 (MQ=255)
                                                                   cTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCa                                                                                    <  2:3600339/150‑1 (MQ=255)
                                                                   cTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCa                                                                                    <  1:1924693/150‑1 (MQ=255)
                                                                   cTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCa                                                                                    <  2:3834931/150‑1 (MQ=255)
                                                                      tgACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACa                                                                                 >  1:5122574/1‑150 (MQ=255)
                                                                      tgACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACa                                                                                 >  1:4410980/1‑150 (MQ=255)
                                                                      tgACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACa                                                                                 <  2:4797412/150‑1 (MQ=255)
                                                                       gACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAg                                                                                <  2:1787674/150‑1 (MQ=255)
                                                                         cGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCt                                                                              >  1:1972211/1‑150 (MQ=255)
                                                                            ccTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTaa                                                                           <  2:1128735/150‑1 (MQ=255)
                                                                                  aTACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCAt                                                                     <  2:895449/150‑1 (MQ=255)
                                                                                  aTACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCAt                                                                     <  1:2780151/150‑1 (MQ=255)
                                                                                    aCGCCTTCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTa                                                                   <  2:3207986/150‑1 (MQ=255)
                                                                                     cGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTaa                                                                  <  2:707469/150‑1 (MQ=255)
                                                                                     cGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTaa                                                                  >  1:3510080/1‑150 (MQ=255)
                                                                                      gCCTTCAATCAGAGAAACCTTACCTTCGACTACCTTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACAATGGTACGCAACATTATCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAAc                                                                 <  2:439308/150‑1 (MQ=255)
                                                                                       ccTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACa                                                                >  1:185731/1‑150 (MQ=255)
                                                                                         tGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATc                                                              <  2:4710229/150‑1 (MQ=255)
                                                                                           cAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCaa                                                            >  1:1274507/1‑150 (MQ=255)
                                                                                            aaTCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGTTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAAt                                                           >  1:3074549/1‑150 (MQ=255)
                                                                                            aaTCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAAt                                                           >  2:4131672/1‑150 (MQ=255)
                                                                                             aTCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAAtt                                                          >  1:1680944/1‑150 (MQ=255)
                                                                                               cAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCCTGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGc                                                        <  2:1438338/150‑1 (MQ=255)
                                                                                               cAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGc                                                        <  2:3635661/150‑1 (MQ=255)
                                                                                                 gagaAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCt                                                      <  1:1545865/150‑1 (MQ=255)
                                                                                                 gagaAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCt                                                      <  2:5170109/150‑1 (MQ=255)
                                                                                                  agaAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCtt                                                     <  1:3329015/150‑1 (MQ=255)
                                                                                                   gaAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTa                                                    <  1:171938/150‑1 (MQ=255)
                                                                                                    aaaCCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTAc                                                   <  1:924605/150‑1 (MQ=255)
                                                                                                     aaCCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACt                                                  >  2:5002842/1‑150 (MQ=255)
                                                                                                     aaCCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACt                                                  <  1:2960254/150‑1 (MQ=255)
                                                                                                     aaCCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACt                                                  <  2:427387/150‑1 (MQ=255)
                                                                                                      accttaccttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTg                                                 <  1:2480502/150‑1 (MQ=255)
                                                                                                         ttaccttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGaa                                              >  2:2101476/1‑150 (MQ=255)
                                                                                                         ttaccttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGaa                                              >  2:105950/1‑150 (MQ=255)
                                                                                                          taccttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGaaa                                             <  1:319500/150‑1 (MQ=255)
                                                                                                           accttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAg                                            >  2:525510/1‑150 (MQ=255)
                                                                                                            ccttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCAGGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGa                                           <  2:1134778/150‑1 (MQ=255)
                                                                                                             cttCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGAt                                          >  1:1654439/1‑150 (MQ=255)
                                                                                                                cGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTg                                       >  2:1345804/1‑150 (MQ=255)
                                                                                                                     aCCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGctc                                  <  2:1390958/150‑1 (MQ=255)
                                                                                                                         ttGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACCGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCt                              <  2:161122/150‑1 (MQ=255)
                                                                                                                         ttGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTTCGCTCTCCt                              >  2:457695/1‑150 (MQ=255)
                                                                                                                         ttGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCt                              >  1:1797532/1‑150 (MQ=255)
                                                                                                                          tGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCtg                             >  1:3990160/1‑150 (MQ=255)
                                                                                                                              aGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCg                         <  2:3875134/150‑1 (MQ=255)
                                                                                                                                  cGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAg                     >  1:133750/1‑150 (MQ=255)
                                                                                                                                   gATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCCCCCCTGGTACGCAACTTTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATGAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGc                    <  1:784118/150‑1 (MQ=255)
                                                                                                                                   gATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGc                    >  2:1976451/1‑150 (MQ=255)
                                                                                                                                     tAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCtt                  <  1:5201505/150‑1 (MQ=255)
                                                                                                                                     tAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAAGTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCtt                  <  2:3020086/150‑1 (MQ=255)
                                                                                                                                      aaTTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTg                 >  1:4514463/1‑150 (MQ=255)
                                                                                                                                      aaTTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTg                 <  1:852685/150‑1 (MQ=255)
                                                                                                                                      aaTTGTCGACCCAAGTAGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAAAAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGAGTTGCTCTCACCAGTGCGGCAGCtag                 >  2:3441056/1‑148 (MQ=255)
                                                                                                                                         tGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTc              >  2:1629271/1‑150 (MQ=255)
                                                                                                                                          gacgcccccAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGATAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCt             <  2:5225402/145‑1 (MQ=255)
                                                                                                                                          gTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCt             <  2:2926312/150‑1 (MQ=255)
                                                                                                                                          gTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCt             <  1:4300750/150‑1 (MQ=255)
                                                                                                                                            cGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCttt           >  1:2658897/1‑150 (MQ=255)
                                                                                                                                             gACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCTTTa          >  1:933030/1‑150 (MQ=255)
                                                                                                                                               ccccAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCTTTAAg        <  2:5238936/150‑1 (MQ=255)
                                                                                                                                                     cTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCTTTAAGCtgttg  >  2:51768/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCCACCACCCTTGCCGCCCACCTGCTGAGCGACCATACCAATCAGTTCCCCTGCTTTCACACGATCTGTGACGTCCTTAGATACGCCTGCAATCAGAGAAACCTTACCTTCGACTACCGTTGCCAGCACGATAATTGTCGACCCCAGCTGATTTTTTAAATCGTCAACCATGGTACGCAACATTTTCGGCTCAACACCGCTAAGCTCGCTAACCAACAGCTTAACACCATTAACATCAATTGCCTTACTGGAAAGATTTGCGCTCTCCTGTGCGGCAGCTTGTTCTTTAAGCTGTTG  >  CP009273/2812834‑2813132

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