Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 325,740 G→T 94.4% P189P (CCG→CCT betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273325,7400GT94.4% 696.5 / 33.2 249P189P (CCG→CCTbetTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base G (2/12);  new base T (149/86);  total (151/98)
Fisher's exact test for biased strand distribution p-value = 3.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 250 total aligned reads displayed.

GGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTG  >  CP009273/325595‑325884
                                                                                                                                                 |                                                                                                                                                
ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTAt                                                                                                                                                <  1:2776370/148‑1 (MQ=255)
         ttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTc                                                                                                                                       >  2:1073792/1‑148 (MQ=255)
           cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAg                                                                                                                                     <  1:2999612/148‑1 (MQ=255)
            aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGt                                                                                                                                    <  2:2300795/148‑1 (MQ=255)
             cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTg                                                                                                                                   >  2:2514704/1‑148 (MQ=255)
              tACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTgg                                                                                                                                  >  1:78102/1‑148 (MQ=255)
               aCGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGa                                                                                                                                 <  1:1014732/148‑1 (MQ=255)
               aCGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGa                                                                                                                                 >  1:2770180/1‑148 (MQ=255)
                 ggCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGata                                                                                                                               >  2:2890140/1‑148 (MQ=255)
                      aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCa                                                                                                                          <  1:1297092/148‑1 (MQ=255)
                       aCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAg                                                                                                                         <  2:3253883/148‑1 (MQ=255)
                        ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGc                                                                                                                        >  1:969008/1‑148 (MQ=255)
                        ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGc                                                                                                                        >  1:969137/1‑148 (MQ=255)
                         cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCg                                                                                                                       <  1:2425584/148‑1 (MQ=255)
                         cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCg                                                                                                                       <  1:2126094/148‑1 (MQ=255)
                         cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCg                                                                                                                       <  2:2007461/148‑1 (MQ=255)
                            cTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTg                                                                                                                    <  1:1830722/148‑1 (MQ=255)
                             tGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGa                                                                                                                   <  1:823342/148‑1 (MQ=255)
                              ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGAt                                                                                                                  >  2:3203303/1‑148 (MQ=255)
                                tCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCg                                                                                                                <  1:1821032/148‑1 (MQ=255)
                                  gATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGc                                                                                                              >  2:1170669/1‑148 (MQ=255)
                                   aTGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCa                                                                                                             >  1:251370/1‑148 (MQ=255)
                                        tGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATc                                                                                                        <  2:3047584/148‑1 (MQ=255)
                                          cgcTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCtt                                                                                                      <  1:3197196/148‑1 (MQ=255)
                                           gcTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTc                                                                                                     <  1:3212780/148‑1 (MQ=255)
                                             tGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCgg                                                                                                   >  1:135289/1‑148 (MQ=255)
                                             tGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCgg                                                                                                   >  1:55779/1‑148 (MQ=255)
                                              gATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGt                                                                                                  <  1:2248691/148‑1 (MQ=255)
                                               aTGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTa                                                                                                 >  1:1999513/1‑148 (MQ=255)
                                                      tGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCAc                                                                                          >  1:817262/1‑148 (MQ=255)
                                                             cGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTc                                                                                   >  1:2230577/1‑148 (MQ=255)
                                                             cGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTc                                                                                   >  2:700000/1‑148 (MQ=255)
                                                               gATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCgg                                                                                 >  1:2174720/1‑148 (MQ=255)
                                                                atatTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGt                                                                                >  1:1895350/1‑148 (MQ=255)
                                                                atatTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGt                                                                                >  1:1895729/1‑148 (MQ=255)
                                                                aTACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGt                                                                                >  2:654200/1‑148 (MQ=255)
                                                                 tACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTa                                                                               >  1:2694761/1‑148 (MQ=255)
                                                                  aCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTAt                                                                              <  2:800544/148‑1 (MQ=255)
                                                                  aCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTAt                                                                              >  2:1527961/1‑148 (MQ=255)
                                                                     ttAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCgg                                                                           >  1:2905014/1‑148 (MQ=255)
                                                                     ttAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCgg                                                                           <  1:3100262/148‑1 (MQ=255)
                                                                     ttAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCgg                                                                           >  1:1070757/1‑148 (MQ=255)
                                                                     ttAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCgg                                                                           >  1:193299/1‑148 (MQ=255)
                                                                        gCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGgtgt                                                                        >  1:1319368/1‑148 (MQ=255)
                                                                          tATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTgtgg                                                                      <  1:1183616/148‑1 (MQ=255)
                                                                           aTCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTgtggt                                                                     >  1:2328809/1‑148 (MQ=255)
                                                                           aTCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTgtggt                                                                     <  2:1899085/148‑1 (MQ=255)
                                                                             cGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGc                                                                   >  2:2636359/1‑148 (MQ=255)
                                                                             cGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGc                                                                   >  1:1027764/1‑148 (MQ=255)
                                                                                tataATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGc                                                                >  1:3329588/1‑148 (MQ=255)
                                                                                 ataATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCt                                                               >  1:656634/1‑148 (MQ=255)
                                                                                    aTTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTaa                                                            >  2:2692583/1‑148 (MQ=255)
                                                                                    aTTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTaa                                                            >  2:2481440/1‑148 (MQ=255)
                                                                                     tttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAAc                                                           >  1:1518529/1‑148 (MQ=255)
                                                                                     tttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAAc                                                           >  1:1518610/1‑148 (MQ=255)
                                                                                      ttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACt                                                          <  1:2345600/148‑1 (MQ=255)
                                                                                      ttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACt                                                          >  2:2471183/1‑148 (MQ=255)
                                                                                      ttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACt                                                          >  2:2467091/1‑148 (MQ=255)
                                                                                         ccGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATg                                                       >  1:2002593/1‑148 (MQ=255)
                                                                                            cTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCt                                                    >  1:3284023/1‑148 (MQ=255)
                                                                                            cTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCt                                                    <  2:2230577/148‑1 (MQ=255)
                                                                                            cTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCt                                                    <  2:1609746/148‑1 (MQ=255)
                                                                                              cACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTg                                                  <  1:2810645/148‑1 (MQ=255)
                                                                                               aCCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGa                                                 <  1:529026/148‑1 (MQ=255)
                                                                                                   tCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGt                                             >  2:857919/1‑148 (MQ=255)
                                                                                                   tCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGt                                             >  1:1171868/1‑148 (MQ=255)
                                                                                                   tCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGt                                             >  1:194758/1‑148 (MQ=255)
                                                                                                   tCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGt                                             >  1:3269821/1‑148 (MQ=255)
                                                                                                    ccGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGTCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTa                                            >  2:917431/1‑148 (MQ=255)
                                                                                                    ccGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTa                                            >  2:2993837/1‑148 (MQ=255)
                                                                                                        tCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGt                                        <  2:2382651/148‑1 (MQ=255)
                                                                                                        tCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGt                                        <  2:49876/148‑1 (MQ=255)
                                                                                                        tCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGt                                        <  2:1284276/148‑1 (MQ=255)
                                                                                                         cGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGtt                                       >  2:633854/1‑148 (MQ=255)
                                                                                                         cGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGtt                                       >  2:552550/1‑148 (MQ=255)
                                                                                                          ggCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGttt                                      >  1:3356558/1‑148 (MQ=255)
                                                                                                           gcgcTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTg                                     <  2:1679725/148‑1 (MQ=255)
                                                                                                            cgcTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGa                                    <  1:476481/148‑1 (MQ=255)
                                                                                                              cTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGata                                  <  1:2584443/148‑1 (MQ=255)
                                                                                                              cTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGata                                  <  2:393912/148‑1 (MQ=255)
                                                                                                                gTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATAtt                                >  2:3403467/1‑148 (MQ=255)
                                                                                                                  aCCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTcc                              <  1:1073792/148‑1 (MQ=255)
                                                                                                                  aCCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTcc                              >  2:2853582/1‑148 (MQ=255)
                                                                                                                    ccGATCTTCGGTAAACGGATTTACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTTCGCTCGGTATCGGTGTGGTGCTGCTTAACTATGGCTTGAGCGTACTGTTTGATAttacat                            <  2:1305972/148‑5 (MQ=255)
                                                                                                                      gATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCTATACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTTCTGTTTGATATTCCCgtt                          <  2:1731799/148‑3 (MQ=255)
                                                                                                                         cTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCg                       >  2:3326131/1‑148 (MQ=255)
                                                                                                                         cTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCg                       >  1:2849071/1‑148 (MQ=255)
                                                                                                                          ttCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGa                      <  1:1283864/148‑1 (MQ=255)
                                                                                                                            cGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATg                    <  2:1999513/148‑1 (MQ=255)
                                                                                                                              gTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATggc                  >  2:1603104/1‑148 (MQ=255)
                                                                                                                               tacaCGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTCCGCTCGGTNTCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATNTTCCCGATTCGATggcg                 <  2:2586166/145‑1 (MQ=255)
                                                                                                                                  aCGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGgcggcg              >  1:387477/1‑148 (MQ=255)
                                                                                                                                     gATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGaaa           <  2:78102/148‑1 (MQ=255)
                                                                                                                                       ttAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGc         <  2:1399243/148‑1 (MQ=255)
                                                                                                                                       ttAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGc         <  2:1894955/148‑1 (MQ=255)
                                                                                                                                        tAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCg        >  2:1538921/1‑148 (MQ=255)
                                                                                                                                         aaCGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCgg       <  1:3042248/148‑1 (MQ=255)
                                                                                                                                          aCGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGc      >  2:1250844/1‑148 (MQ=255)
                                                                                                                                              gCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTg  >  2:523637/1‑148 (MQ=255)
                                                                                                                                              gCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTg  >  2:524140/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTG  >  CP009273/325595‑325884

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: