Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,046,172 A→T 86.6% K256N (AAA→AAT hemN → coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,046,1720AT86.6% 600.5 / 74.0 201K256N (AAA→AAThemNcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base A (11/16);  new base T (99/75);  total (110/91)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 206 total aligned reads displayed.

CATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCT  >  CP009273/4046025‑4046321
                                                                                                                                                   |                                                                                                                                                     
 aTCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGAt                                                                                                                                                    <  2:2429608/150‑1 (MQ=255)
  tCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAatgatg                                                                                                                                                   <  1:3968871/150‑1 (MQ=255)
  tCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAatgatg                                                                                                                                                   <  2:3680886/150‑1 (MQ=255)
   cGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGc                                                                                                                                                  >  1:157346/1‑150 (MQ=255)
   cGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGc                                                                                                                                                  <  2:561749/150‑1 (MQ=255)
    gACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                                 <  2:1248171/150‑1 (MQ=255)
    gACCTGATTAACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCt                                                                                                                                                 >  2:784563/1‑150 (MQ=255)
      ccTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGa                                                                                                                                               <  2:3312195/150‑1 (MQ=255)
           tttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGc                                                                                                                                          >  1:4640276/1‑150 (MQ=255)
           tttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGc                                                                                                                                          <  1:1190923/150‑1 (MQ=255)
            ttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGcc                                                                                                                                         <  2:2773979/150‑1 (MQ=255)
            ttACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGcc                                                                                                                                         <  2:1295789/150‑1 (MQ=255)
             tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCg                                                                                                                                        <  2:4382443/150‑1 (MQ=255)
               cGGCCTGCCGAAACAGACGCCGGAGAGTTTCACCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCCGCGCAAAATCAATGATGCTGACCTGCCGAg                                                                                                                                      >  2:4484380/1‑150 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                     >  1:1782665/1‑150 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                     >  1:1782709/1‑150 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                     <  2:2688125/150‑1 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                     <  2:4325111/150‑1 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGt                                                                                                                                     <  1:2236445/150‑1 (MQ=255)
                ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGt                                                                                                                                     <  1:4536196/150‑1 (MQ=255)
                 gCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTc                                                                                                                                    <  2:2972797/150‑1 (MQ=255)
                   cTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCg                                                                                                                                  >  2:1673388/1‑150 (MQ=255)
                    tGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgc                                                                                                                                 >  1:47258/1‑150 (MQ=255)
                      ccGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCgcag                                                                                                                               <  2:4340001/150‑1 (MQ=255)
                        gAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgcagca                                                                                                                             >  2:3627739/1‑150 (MQ=255)
                            cAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCaaaaa                                                                                                                         >  2:3856241/1‑150 (MQ=255)
                                cGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                                     <  2:4374406/150‑1 (MQ=255)
                                 gCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGa                                                                                                                    >  2:3047733/1‑150 (MQ=255)
                                 gCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGa                                                                                                                    >  2:3047741/1‑150 (MQ=255)
                                   cGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGata                                                                                                                  <  1:734663/150‑1 (MQ=255)
                                     gtgagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATc                                                                                                                <  1:3832611/148‑1 (MQ=255)
                                      agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATcc                                                                                                               >  1:3090973/1‑150 (MQ=255)
                                      agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATcc                                                                                                               >  1:768835/1‑150 (MQ=255)
                                      agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATcc                                                                                                               >  1:3401845/1‑150 (MQ=255)
                                          tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCa                                                                                                           <  2:3101352/150‑1 (MQ=255)
                                           tttgCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAg                                                                                                          <  1:1246290/147‑1 (MQ=255)
                                           ttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAg                                                                                                          <  1:3189307/150‑1 (MQ=255)
                                           ttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAg                                                                                                          <  1:3557016/150‑1 (MQ=255)
                                            tCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAgg                                                                                                         <  2:2367058/150‑1 (MQ=255)
                                              gCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                                                                       >  2:3573630/1‑150 (MQ=255)
                                              gCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                                                                       >  1:2850307/1‑150 (MQ=255)
                                              gCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                                                                       <  2:3983326/150‑1 (MQ=255)
                                              gCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaa                                                                                                       <  2:3983106/150‑1 (MQ=255)
                                               ccTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaaa                                                                                                      >  1:2911755/1‑150 (MQ=255)
                                               ccTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGaaa                                                                                                      >  2:1884668/1‑150 (MQ=255)
                                                cTTTACCCTGAACCGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                                     <  2:1923712/150‑1 (MQ=255)
                                                 tttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAcc                                                                                                    >  1:3047947/1‑150 (MQ=255)
                                                    aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                                                                 <  1:1516477/150‑1 (MQ=255)
                                                    aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                                                                 <  1:2500376/150‑1 (MQ=255)
                                                    aCCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                                                                 <  1:2500755/150‑1 (MQ=255)
                                                       cTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGcc                                                                                              <  1:4172886/150‑1 (MQ=255)
                                                             cGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTg                                                                                        >  2:762490/1‑150 (MQ=255)
                                                             cGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTg                                                                                        >  1:1877753/1‑150 (MQ=255)
                                                                    cGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                                                                 >  1:1142397/1‑150 (MQ=255)
                                                                    cGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                                                                 >  1:2579498/1‑150 (MQ=255)
                                                                     ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATc                                                                                <  2:651349/150‑1 (MQ=255)
                                                                     ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATc                                                                                >  1:39831/1‑150 (MQ=255)
                                                                       aGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                              >  1:1133298/1‑150 (MQ=255)
                                                                           gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                                                          <  1:1645409/150‑1 (MQ=255)
                                                                           gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                                                          <  1:1645891/150‑1 (MQ=255)
                                                                            aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                         >  1:269620/1‑150 (MQ=255)
                                                                               cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                      >  2:4101890/1‑150 (MQ=255)
                                                                               cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                      >  1:2790010/1‑150 (MQ=255)
                                                                               cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                      >  1:2906457/1‑150 (MQ=255)
                                                                               cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                      >  1:3540868/1‑150 (MQ=255)
                                                                                ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGt                                                                     <  2:2230008/150‑1 (MQ=255)
                                                                                  gATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGttt                                                                   >  2:4174900/1‑150 (MQ=255)
                                                                                   aTCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                                  >  1:4264661/1‑150 (MQ=255)
                                                                                      gTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCg                                                               <  2:863190/150‑1 (MQ=255)
                                                                                         tGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTa                                                            >  1:1005092/1‑150 (MQ=255)
                                                                                           aGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATg                                                          >  1:926287/1‑150 (MQ=255)
                                                                                             tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGa                                                        >  1:2269311/1‑150 (MQ=255)
                                                                                             tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCGGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGACACCATCGCCTTCCTGTCGCAATCGGGCTATCAGTTTATCGGTATGGa                                                        <  1:2110351/150‑1 (MQ=255)
                                                                                                   tAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACtt                                                  >  1:1499507/1‑150 (MQ=255)
                                                                                                       tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGcc                                              >  1:4104775/1‑150 (MQ=255)
                                                                                                            gcATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                         >  1:3877538/1‑150 (MQ=255)
                                                                                                             cATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCg                                        <  1:3374808/150‑1 (MQ=255)
                                                                                                               tCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGa                                      <  2:1116791/150‑1 (MQ=255)
                                                                                                                cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                     >  1:3491527/1‑150 (MQ=255)
                                                                                                                cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                     >  1:3491448/1‑150 (MQ=255)
                                                                                                                cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                     <  1:1525133/150‑1 (MQ=255)
                                                                                                                 tGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATg                                    >  1:2012263/1‑150 (MQ=255)
                                                                                                                   ccGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGAc                                  >  1:1406631/1‑150 (MQ=255)
                                                                                                                    cGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                                 <  2:4216789/150‑1 (MQ=255)
                                                                                                                    cGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                                 <  2:3847325/150‑1 (MQ=255)
                                                                                                                     gACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                                >  1:329604/1‑150 (MQ=255)
                                                                                                                         aTTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                            >  2:4033758/1‑150 (MQ=255)
                                                                                                                         aTTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                            >  2:1819715/1‑150 (MQ=255)
                                                                                                                         aTTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                            >  2:281311/1‑150 (MQ=255)
                                                                                                                          tttttGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                           >  1:431884/1‑150 (MQ=255)
                                                                                                                               gctgctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                      <  1:4222130/150‑1 (MQ=255)
                                                                                                                                    tCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa                 <  1:2192236/150‑1 (MQ=255)
                                                                                                                                     cAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAg                >  2:2059239/1‑150 (MQ=255)
                                                                                                                                      aGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGc               >  1:4377637/1‑150 (MQ=255)
                                                                                                                                      aGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGc               >  2:360382/1‑150 (MQ=255)
                                                                                                                                        cgcAAAATCAATGATGCTTACCTTCCGAGTCCGCAGCAAAAACTCGAAATCCTGAAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt             >  2:657711/1‑150 (MQ=255)
                                                                                                                                            aaaaCAATGATGCTGACCTGCCGCGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAg         <  2:254496/150‑1 (MQ=255)
                                                                                                                                              aTCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGc       >  1:1844813/1‑150 (MQ=255)
                                                                                                                                                cAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGt     <  2:2126662/150‑1 (MQ=255)
                                                                                                                                                   aGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgct  >  1:1562213/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCT  >  CP009273/4046025‑4046321

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: