Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 2,209,361 A→G H98R (CAC→CGC)  suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,209,3610AG99.0% 1212.6 / ‑0.9 302H98R (CAC→CGC) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (194/105);  total (196/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 323 total aligned reads displayed.

GCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCACTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGCTCGC  >  AP012306/2209213‑2209506
                                                                                                                                                     |                                                                                                                                                  
 cGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCt                                                                                                                                                  <  1:4137510/150‑1 (MQ=255)
   gTGATTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTc                                                                                                                                                <  1:1479197/150‑1 (MQ=255)
      aTTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTcgcg                                                                                                                                             >  2:3579539/1‑150 (MQ=255)
      aTTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTcgcg                                                                                                                                             <  2:2732352/150‑1 (MQ=255)
       ttATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCgg                                                                                                                                            <  1:945381/150‑1 (MQ=255)
       ttATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCgg                                                                                                                                            >  1:860494/1‑150 (MQ=255)
          tCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGtat                                                                                                                                         <  1:98925/150‑1 (MQ=255)
           cGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGtatc                                                                                                                                        >  2:1258059/1‑150 (MQ=255)
            gACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGtatct                                                                                                                                       >  2:3720635/1‑150 (MQ=255)
               acGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGtatctatc                                                                                                                                    <  1:1324731/150‑1 (MQ=255)
               acGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGtatctatc                                                                                                                                    >  2:164346/1‑150 (MQ=255)
                cGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCg                                                                                                                                   <  2:2698910/150‑1 (MQ=255)
                cGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCg                                                                                                                                   >  1:366632/1‑150 (MQ=255)
                 gATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGc                                                                                                                                  >  2:1775243/1‑150 (MQ=255)
                  aTTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCt                                                                                                                                 >  1:3221615/1‑150 (MQ=255)
                  aTTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCt                                                                                                                                 >  1:3392001/1‑150 (MQ=255)
                    tCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGt                                                                                                                               >  1:1550233/1‑150 (MQ=255)
                      gTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTc                                                                                                                             <  1:2921972/150‑1 (MQ=255)
                       tAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCg                                                                                                                            <  2:2914335/150‑1 (MQ=255)
                        aaaTCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGt                                                                                                                           >  2:2296126/1‑150 (MQ=255)
                        aaaTCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGt                                                                                                                           <  1:1054714/150‑1 (MQ=255)
                         aaTCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTa                                                                                                                          >  2:2767916/1‑150 (MQ=255)
                         aaTCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTa                                                                                                                          >  2:2517594/1‑150 (MQ=255)
                          atattACCCACAGCGCACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTGCAATGGGTTATCGATCCACTGTATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTTTCTCTCGCTGTTCGTAt                                                                                                                         <  2:3195737/147‑1 (MQ=255)
                           tCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGGTGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATc                                                                                                                        <  1:311898/150‑1 (MQ=255)
                           tCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATc                                                                                                                        <  1:564730/150‑1 (MQ=255)
                               aCCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAg                                                                                                                    >  1:2304359/1‑150 (MQ=255)
                                cccACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAgg                                                                                                                   >  1:2517788/1‑150 (MQ=255)
                                 ccACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGc                                                                                                                 >  1:2060294/1‑150 (MQ=255)
                                 ccACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGc                                                                                                                 >  1:3920952/1‑150 (MQ=255)
                                     aGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCa                                                                                                             >  1:1123091/1‑150 (MQ=255)
                                      gCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAACCGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCAc                                                                                                            <  2:1205585/150‑1 (MQ=255)
                                       cacaCCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCAcc                                                                                                           >  1:2390400/1‑150 (MQ=255)
                                        acacCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCg                                                                                                          <  1:515490/150‑1 (MQ=255)
                                            catcatCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGt                                                                                                      >  1:2337350/1‑150 (MQ=255)
                                             atcatcACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGtt                                                                                                     <  2:4105004/150‑1 (MQ=255)
                                              tcatcaCCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTg                                                                                                    >  1:2731555/1‑150 (MQ=255)
                                               catcaCCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGc                                                                                                   >  2:2646133/1‑150 (MQ=255)
                                                  caCCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCtgt                                                                                                <  2:3165742/150‑1 (MQ=255)
                                                   aCCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCtgtg                                                                                               >  2:3816357/1‑150 (MQ=255)
                                                   aCCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCACTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCtgtg                                                                                               <  2:3891607/150‑1 (MQ=255)
                                                    ccGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTgg                                                                                              <  2:1823106/150‑1 (MQ=255)
                                                      gaagaaAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGta                                                                                            >  2:3102457/1‑150 (MQ=255)
                                                        agaaAGCGATGAACTTGAAGGTACTGATCAGGATCTTCAATTGGTTAGCGATCGACTGGATGTCGCCACCAACCTGATCAAACGTCTGCCCAGGTTCGCGGTATCAAGCGCTGCTCGTCTCAAAGACCGCACCGAAGTAGCTGCGGtata                                                                                          <  2:1302053/150‑1 (MQ=255)
                                                         gaaAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATAc                                                                                         <  1:1638816/150‑1 (MQ=255)
                                                           aaGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGa                                                                                       >  1:3650148/1‑150 (MQ=255)
                                                             gCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATc                                                                                     <  2:1496550/150‑1 (MQ=255)
                                                             gCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATc                                                                                     <  1:1774832/150‑1 (MQ=255)
                                                               ggTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCt                                                                                   <  2:2815463/150‑1 (MQ=255)
                                                                gTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTa                                                                                  <  2:4276899/150‑1 (MQ=255)
                                                                 tGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTAt                                                                                 >  2:894098/1‑150 (MQ=255)
                                                                  gAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATg                                                                                <  2:1716574/150‑1 (MQ=255)
                                                                   aaCTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGc                                                                               <  2:3286121/150‑1 (MQ=255)
                                                                    aCTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTTCCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCg                                                                              <  1:3997215/150‑1 (MQ=255)
                                                                     cTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGt                                                                             <  1:100269/150‑1 (MQ=255)
                                                                     cTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGt                                                                             >  1:1609614/1‑150 (MQ=255)
                                                                     cTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGt                                                                             >  1:1675008/1‑150 (MQ=255)
                                                                     cTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGt                                                                             >  2:1583308/1‑150 (MQ=255)
                                                                     cTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGt                                                                             <  2:844895/150‑1 (MQ=255)
                                                                       tGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTaa                                                                           >  2:3538888/1‑150 (MQ=255)
                                                                          aGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGa                                                                        >  1:3274616/1‑150 (MQ=255)
                                                                             tACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACt                                                                     <  1:3576918/150‑1 (MQ=255)
                                                                              aCTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTg                                                                    <  2:4237899/150‑1 (MQ=255)
                                                                              aCTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTg                                                                    >  1:4221503/1‑150 (MQ=255)
                                                                                  aTCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCa                                                                >  1:3371184/1‑150 (MQ=255)
                                                                                    cAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCAcc                                                              >  1:814170/1‑150 (MQ=255)
                                                                                        atcttCACTGGGTTATCGATCCACTGGATGGCACTCCCACCTTTATCAACCGTCTGCCGCGCTTCGCGGTATCTATCCCTGTTCGTATCAAAGCCCGCCCCGAAGTTGCTGTGGTATCCGATCCTATGCGTAACGAACTGTTCACCGCCa                                                          <  2:7638/147‑1 (MQ=255)
                                                                                        aTGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCa                                                          >  1:3030798/1‑150 (MQ=255)
                                                                                        aTGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCa                                                          <  2:544544/150‑1 (MQ=255)
                                                                                         tGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCAc                                                         >  2:2125144/1‑150 (MQ=255)
                                                                                          gTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACt                                                        >  1:258965/1‑150 (MQ=255)
                                                                                           ttCAATGGGTTTCCGCTCCACTGGATGGCACTCCCACCTTTATCAACCGCCTGCCGCGCTCCGCGGTTTCTCTCGCTGTTCGTATCAAGGGCCCCCCCGAAGTTGGGGTGGTATACGTTCCTATGCGTAACGAACTGTCCACCGCCACTc                                                       <  2:3211622/150‑1 (MQ=255)
                                                                                            tCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTcg                                                      >  1:2219572/1‑150 (MQ=255)
                                                                                             taaTGGGTTTTCCGACCCACTGGATGGCACTCCCACCTTTACAAACCGTCTGCCGCCCTTCGCGTTTTCTACCGCTGTTCGTATAAAAGCCCCCCCCGAAGTTGCTGGGGTATCCGTTCCTATGCGAAACAAACTGTTCACCGCCACTcg                                                      <  2:2306789/149‑1 (MQ=255)
                                                                                              aaTGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTcgcg                                                    >  2:1032853/1‑150 (MQ=255)
                                                                                              aaTGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTcgcg                                                    >  1:1376753/1‑150 (MQ=255)
                                                                                                   gTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAg                                               >  2:1925703/1‑150 (MQ=255)
                                                                                                   gTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAg                                               >  2:3960607/1‑150 (MQ=255)
                                                                                                          gctCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGca                                         >  2:4000261/3‑150 (MQ=255)
                                                                                                           tcccaggggATGCCACTCCCAATTTTAAAAAACGTCTCCCGCGCTCCGGGGTTTCTCCCGCTGTTCGTACAAAGGCCCCCCCCGAGGTGGCTGGGGTATCCGATCCTATGCGAAACGACCTTTTCACCGCCACTCGCGGTCAGGGCGcac                                        <  2:1069392/143‑1 (MQ=255)
                                                                                                           aTCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGcac                                        >  1:3777399/1‑150 (MQ=255)
                                                                                                              cACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGc                                     >  1:2201510/1‑150 (MQ=255)
                                                                                                                cTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTg                                   >  2:3080463/1‑150 (MQ=255)
                                                                                                                 tGGATGGCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGa                                  >  1:1930882/1‑150 (MQ=255)
                                                                                                                      gccACTACCATCTTAAAAAACCGTCGGCCGCGCTTCGCGGTATCTACCGTTGTTCGTATCAAAGCCCGCCCCGAAGTTGCTGGGGTATACGACCCTATGCGAAACGACCTGTCCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGc                             <  2:1390919/148‑1 (MQ=255)
                                                                                                                      gccACTACCAATTTTACCAACCGTCGGCCGCGCTTCGCGGTTTCTATCGCTGTTCGTATCAAGGGCCGCCCCGAAGTTGCTGGGGTATACGATCTTATGCGTAACGAACTGTCCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGc                             <  2:3444153/148‑1 (MQ=255)
                                                                                                                       gCACTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCt                            >  1:1017385/1‑150 (MQ=255)
                                                                                                                         aCTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTAc                          >  1:1412824/1‑150 (MQ=255)
                                                                                                                          cTACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTAcc                         >  1:3490647/1‑150 (MQ=255)
                                                                                                                           tACCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCg                        >  2:3631908/1‑150 (MQ=255)
                                                                                                                            aCCAACTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGa                       >  1:3535863/1‑150 (MQ=255)
                                                                                                                             ccAACTTTATCAAACGTCTGCCGCACTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGAc                      >  1:661489/1‑150 (MQ=255)
                                                                                                                                 cTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTgcg                  >  2:637197/1‑150 (MQ=255)
                                                                                                                                 cTTTATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTgcg                  >  2:2408891/1‑150 (MQ=255)
                                                                                                                                   ttATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGcgcg                >  2:2491764/1‑150 (MQ=255)
                                                                                                                                    tATCAAACGTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCgg               >  1:3101250/1‑150 (MQ=255)
                                                                                                                                            gTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCg       >  1:2274115/1‑150 (MQ=255)
                                                                                                                                            gTCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCg       >  2:302291/1‑150 (MQ=255)
                                                                                                                                             tCTGCCGCGCTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGc      >  1:1217037/1‑150 (MQ=255)
                                                                                                                                                 ccGCGCTTCGCGGTATCTACCGCTGTTCTTAAAAAAGGCCGCACCGAGGTTGCGGTGGTATACGACCCTAGGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGCTcgc  <  2:2242538/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCACACCATCATCACCGAAGAAAGCGGTGAACTTGAAGGTACTGATCAGGATGTTCAATGGGTTATCGATCCACTGGATGGCACTACCAACTTTATCAAACGTCTGCCGCACTTCGCGGTATCTATCGCTGTTCGTATCAAAGGCCGCACCGAAGTTGCTGTGGTATACGATCCTATGCGTAACGAACTGTTCACCGCCACTCGCGGTCAGGGCGCACAGCTGAACGGCTACCGACTGCGCGGCAGCACCGCTCGC  >  AP012306/2209213‑2209506

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