Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 860,669 G→T 9.9% I384I (ATC→ATA moeA ← molybdopterin molybdenumtransferase; molybdopterin biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273860,6690GT 9.9% 1108.0 / 56.9 298I384I (ATC→ATAmoeAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (145/123);  new base T (22/8);  total (167/131)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

TCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGT  >  CP009273/860521‑860818
                                                                                                                                                    |                                                                                                                                                     
 cAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGAt                                                                                                                                                     >  1:9631/1‑150 (MQ=255)
   aaGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATaa                                                                                                                                                   >  1:4303192/1‑150 (MQ=255)
       cgcgCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCa                                                                                                                                               <  2:1671176/150‑1 (MQ=255)
        gagaaGAATGATTAGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAAAGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTACCACAAAGCCGCGATCGCGTTCCAGCACGATAAAGCAg                                                                                                                                              >  2:3726158/5‑150 (MQ=255)
          gcAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGtt                                                                                                                                            <  1:3853160/150‑1 (MQ=255)
          gcAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGtt                                                                                                                                            >  2:2002642/1‑150 (MQ=255)
           cAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTg                                                                                                                                           <  1:3722747/150‑1 (MQ=255)
           cAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTg                                                                                                                                           >  2:616358/1‑150 (MQ=255)
            aGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGc                                                                                                                                          >  2:1202312/1‑150 (MQ=255)
             gAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGcc                                                                                                                                         <  1:895911/150‑1 (MQ=255)
                tGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAg                                                                                                                                      <  2:3498022/150‑1 (MQ=255)
                tGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAg                                                                                                                                      <  1:1248209/150‑1 (MQ=255)
                 gATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAgg                                                                                                                                     >  2:509816/1‑150 (MQ=255)
                 gATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAgg                                                                                                                                     >  1:4320203/1‑150 (MQ=255)
                  aTTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGc                                                                                                                                    >  1:4386578/1‑150 (MQ=255)
                  aTTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGc                                                                                                                                    >  2:1243087/1‑150 (MQ=255)
                  aTTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTAAAGGCCTCCGAACAACGCTTTAAAAGGTTCTACTTCAAACCATTCGCCAACTTCCACATTGAAGCTATCGCGTTCCAGCACGATAAATCAGTTGCCTAgga                                                                                                                                    >  1:4524401/1‑149 (MQ=255)
                    ttGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTa                                                                                                                                  <  2:1695777/150‑1 (MQ=255)
                      gCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTaaa                                                                                                                                <  2:4374476/150‑1 (MQ=255)
                         ggTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGa                                                                                                                             <  1:2852846/150‑1 (MQ=255)
                          gTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAg                                                                                                                            <  2:3640438/150‑1 (MQ=255)
                          gTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAg                                                                                                                            <  2:3935100/150‑1 (MQ=255)
                              tAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTaa                                                                                                                        >  1:2820627/1‑150 (MQ=255)
                                 cgcAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                                                                                     <  1:2228663/150‑1 (MQ=255)
                                 cgcAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAata                                                                                                                     <  1:2610675/150‑1 (MQ=255)
                                  gcagcaTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAatat                                                                                                                    <  1:3603816/150‑1 (MQ=255)
                                      caTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATAtgt                                                                                                                <  2:4075827/150‑1 (MQ=255)
                                       aTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCGCCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATAtgtg                                                                                                               >  2:1442705/1‑150 (MQ=255)
                                       aTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATAtgtg                                                                                                               >  1:995409/1‑150 (MQ=255)
                                       aTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATAtgtg                                                                                                               >  2:695200/1‑150 (MQ=255)
                                                     gAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCg                                                                                                 >  1:2582608/1‑150 (MQ=255)
                                                          ccGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCggtggt                                                                                            >  1:3736880/1‑150 (MQ=255)
                                                          ccGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCggtggt                                                                                            >  1:3951393/1‑150 (MQ=255)
                                                           cGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtc                                                                                           <  2:774748/150‑1 (MQ=255)
                                                           cGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtc                                                                                           <  2:3776689/150‑1 (MQ=255)
                                                            gCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtcg                                                                                          <  2:1083659/150‑1 (MQ=255)
                                                             ccATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtcgt                                                                                         >  1:646328/1‑150 (MQ=255)
                                                             ccATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtcgt                                                                                         >  1:3710016/1‑150 (MQ=255)
                                                             caaTTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtcgt                                                                                         >  1:4015606/3‑150 (MQ=255)
                                                              cATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGgtcgtc                                                                                        <  1:1888641/150‑1 (MQ=255)
                                                               aTTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCa                                                                                       >  2:186562/1‑150 (MQ=255)
                                                                    aGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTcc                                                                                  <  1:4094673/150‑1 (MQ=255)
                                                                        cTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGtt                                                                              <  2:1608735/150‑1 (MQ=255)
                                                                         tCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTc                                                                             <  1:3545003/150‑1 (MQ=255)
                                                                         tCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTc                                                                             <  2:4493615/150‑1 (MQ=255)
                                                                          ccGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCg                                                                            >  2:3644338/1‑150 (MQ=255)
                                                                           cGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGc                                                                           >  1:2076704/1‑150 (MQ=255)
                                                                            gAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGcc                                                                          <  2:1100155/150‑1 (MQ=255)
                                                                               caacGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATc                                                                       <  2:157899/150‑1 (MQ=255)
                                                                                   gcgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcg                                                                   <  2:2500847/150‑1 (MQ=255)
                                                                                    cgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgt                                                                  >  1:599015/1‑150 (MQ=255)
                                                                                    cgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgt                                                                  >  1:1156819/1‑150 (MQ=255)
                                                                                     gTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgtt                                                                 <  1:1002149/150‑1 (MQ=255)
                                                                                     gTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgtt                                                                 >  1:3774576/1‑150 (MQ=255)
                                                                                        aaaCGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgttgcg                                                              >  2:1648089/1‑150 (MQ=255)
                                                                                          aCGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGgcgttgcgtt                                                            <  2:2634423/150‑1 (MQ=255)
                                                                                            ggTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCCAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCgttgcgttgc                                                          <  1:1499591/150‑1 (MQ=255)
                                                                                             gTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTgca                                                         <  2:1003334/150‑1 (MQ=255)
                                                                                               tCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTgcagc                                                       >  1:3941506/1‑150 (MQ=255)
                                                                                               tCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTgcagc                                                       >  2:4548112/1‑150 (MQ=255)
                                                                                                cTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTACATATATGTGAACCCTGATGTCCGGTGGTCGACACTTCCAGTTCGCCATCGGCGATGCGTTgcagca                                                      >  2:4154133/1‑150 (MQ=255)
                                                                                                cTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTgcagca                                                      <  1:4135125/150‑1 (MQ=255)
                                                                                                cTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTgcagca                                                      >  2:1163629/1‑150 (MQ=255)
                                                                                                  aCTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCAcg                                                    >  1:3509073/1‑150 (MQ=255)
                                                                                                  aCTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATAGATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCAcg                                                    <  2:3103565/150‑1 (MQ=255)
                                                                                                   cTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCAcgc                                                   >  2:4215746/1‑150 (MQ=255)
                                                                                                       cACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCcgcg                                               <  2:218315/150‑1 (MQ=255)
                                                                                                       cACCCATTCGCCCACGTCCACATTGCCGCGATCGCGTTCCAGCACGATAACGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCcgcg                                               <  1:793321/150‑1 (MQ=255)
                                                                                                         cccATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCt                                             >  1:529180/1‑150 (MQ=255)
                                                                                                          ccATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTg                                            >  1:4166782/1‑150 (MQ=255)
                                                                                                            aTTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGa                                          >  2:3492713/1‑150 (MQ=255)
                                                                                                             ttCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGaa                                         <  1:834858/150‑1 (MQ=255)
                                                                                                                gCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATc                                      >  1:3683523/1‑150 (MQ=255)
                                                                                                                 cccACTTCCACATTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCa                                     >  2:2279966/1‑150 (MQ=255)
                                                                                                                 cccACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCa                                     >  2:4399031/1‑150 (MQ=255)
                                                                                                                   cACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAg                                   >  1:1723598/1‑150 (MQ=255)
                                                                                                                    aCTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGa                                  >  2:2789201/1‑150 (MQ=255)
                                                                                                                       tCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGt                               >  2:947353/1‑150 (MQ=255)
                                                                                                                          acaTTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCt                            >  2:3584930/1‑150 (MQ=255)
                                                                                                                          acaTTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCt                            >  2:3197751/1‑150 (MQ=255)
                                                                                                                           caTTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTg                           <  1:1908125/150‑1 (MQ=255)
                                                                                                                            aTTGCCGCGATCGCGTTCCAGCACTATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTgg                          >  2:3881075/1‑150 (MQ=255)
                                                                                                                                 cgcgATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGttt                     <  2:1367752/150‑1 (MQ=255)
                                                                                                                                  gcgATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGtttt                    <  1:484173/150‑1 (MQ=255)
                                                                                                                                  gcgATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGtttt                    <  1:186562/150‑1 (MQ=255)
                                                                                                                                    gATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGtttttt                  >  2:3272890/1‑150 (MQ=255)
                                                                                                                                    gATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGtttttt                  >  1:199903/1‑150 (MQ=255)
                                                                                                                                      tCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTc                >  2:856112/1‑150 (MQ=255)
                                                                                                                                      tCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTc                <  1:1418579/150‑1 (MQ=255)
                                                                                                                                       cgcgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCa               <  2:216660/150‑1 (MQ=255)
                                                                                                                                         cgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAgg             >  2:1850834/1‑150 (MQ=255)
                                                                                                                                         cgTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAgg             >  1:2184955/1‑150 (MQ=255)
                                                                                                                                          gTTCCAGCTCGTTAAAGCAGTTGCCGAGGATAAAGGAGCTATATATATGTGAACCCTGATGTCCGTTGGTCGTCACTTCCTGTTCGCCTTCGGCGTTGCGTTGCAGCACGCCGGGGTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGc            <  2:3426749/150‑1 (MQ=255)
                                                                                                                                          gTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGc            >  2:4239673/1‑150 (MQ=255)
                                                                                                                                          gTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGc            >  1:3295054/1‑150 (MQ=255)
                                                                                                                                           ttCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCg           >  2:2959103/1‑150 (MQ=255)
                                                                                                                                            tCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCgg          >  1:2535689/1‑150 (MQ=255)
                                                                                                                                             ccAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCggg         >  1:71564/1‑150 (MQ=255)
                                                                                                                                                 cACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGc     >  1:2658766/1‑150 (MQ=255)
                                                                                                                                                   cGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCtg   >  1:1111997/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGT  >  CP009273/860521‑860818

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: