Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,280,489 C→G 6.2% A327P (GCG→CCG)  PP_3751 ← cyanate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,280,4890CG 6.2% 265.8 / 4.9 85A327P (GCG→CCG) PP_3751cyanate transporter
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base G (4/2);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGCGGTCACCAGCAGCATTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGAAGCCCGCC  >  NC_002947/4280343‑4280637
                                                                                                                                                  |                                                                                                                                                     
gCACGCGGTCACCAGCAGCATTGCCATCACGGTCACGGTCAACAGCGTCCACGAATCCTTGCAGTCGCCACGGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACGCTGTGCGAAAGCGGGCAg                                                                                                                                                    >  2:213640/1‑150 (MQ=255)
   cgcgGTCACCAGCAGCATTGCCATCACGGGCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCt                                                                                                                                                 >  2:2085123/1‑150 (MQ=255)
          aCCAGCAGCATTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCt                                                                                                                                         >  2:1052314/1‑150 (MQ=255)
            cagcagCATTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTc                                                                                                                                       <  2:608261/150‑1 (MQ=255)
                agcaTTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCa                                                                                                                                   <  1:725516/150‑1 (MQ=255)
                 gcaTTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGc                                                                                                                                               <  2:783786/138‑1 (MQ=255)
                 gcaTTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGc                                                                                                                                               >  1:783786/1‑138 (MQ=255)
                 gcaTTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAg                                                                                                                                  <  1:182107/150‑1 (MQ=255)
                 gcaTTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAg                                                                                                                                  <  1:1052314/150‑1 (MQ=255)
                       ccATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGttt                                                                                                                            >  1:2309783/1‑150 (MQ=255)
                       ccATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGttt                                                                                                                            >  2:1948290/1‑150 (MQ=255)
                         aTCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTc                                                                                                                          <  2:2089909/150‑1 (MQ=255)
                          tcacggtcacggTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCa                                                                                                                         >  1:1373829/1‑150 (MQ=255)
                          tcacggtcacggTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCa                                                                                                                         >  2:867085/1‑150 (MQ=255)
                           cacggtcacggtCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGTCTGCCGGCTTCGGCAGGCGTTTTCAg                                                                                                                        <  1:1507160/150‑1 (MQ=255)
                                   gggtaaACAGCGTCCACGAGGCTTGAAGGTGGCCACGGACACCCGGAAGCAAACCCGTGATATACGGCAGGAGGCCGGTGAGGATGAGGCCGATACCCGGTACGAACGCCGCCAGGCGGCCGGCTTCGGCAGGCGTTTTCAGGTGTTcag                                                                                                                <  2:2194406/149‑3 (MQ=255)
                                   cggtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCCAGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGACGCCAGGCTGCCGGCTGCGGCAGGCGATTTCAAGTCTTcta                                                                                                                >  2:110051/1‑148 (MQ=255)
                                    ggtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAg                                                                                                               >  2:911291/1‑150 (MQ=255)
                                    ggtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAg                                                                                                               >  2:1320029/1‑150 (MQ=255)
                                     gtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAgg                                                                                                              <  1:840210/150‑1 (MQ=255)
                                     gtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAgg                                                                                                              <  1:525866/150‑1 (MQ=255)
                                     gtcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAgg                                                                                                              <  2:2309783/150‑1 (MQ=255)
                                     ctcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCGCCCGGTGGCAGGCGTTTTCAGGTGTTCCAgg                                                                                                              <  1:1871396/149‑1 (MQ=255)
                                      tcaACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAggg                                                                                                             <  2:578776/150‑1 (MQ=255)
                                        aaCAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATAACATGTGTGCTCAGTGGCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAgggag                                                                                                           >  2:1871396/1‑148 (MQ=255)
                                           aGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAAc                                                                                                        <  1:483057/150‑1 (MQ=255)
                                              gTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTc                                                                                                     >  2:989498/1‑150 (MQ=255)
                                                ccACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAg                                                                                                   <  2:1303171/150‑1 (MQ=255)
                                                 cACGAACACTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAgg                                                                                                  >  2:1198138/1‑150 (MQ=255)
                                                  aCGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGc                                                                                                 >  1:433156/1‑150 (MQ=255)
                                                       gCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTgg                                                                                            <  2:1858889/150‑1 (MQ=255)
                                                                tCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCATGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAg                                                                                   <  1:1948290/150‑1 (MQ=255)
                                                                  gCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGTTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTc                                                                                                                                       >  1:811542/1‑96 (MQ=255)
                                                                  gCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGTTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTc                                                                                                                                       <  2:811542/96‑1 (MQ=255)
                                                                    cACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCAc                                                                               <  1:911291/150‑1 (MQ=255)
                                                                    cACTGACATCCCGCAGCCAACCCGTGATACACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCAc                                                                               <  2:1628592/150‑1 (MQ=255)
                                                                        gACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGTCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCaa                                                                           <  1:205767/150‑1 (MQ=255)
                                                                         aCATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCAGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAg                                                                          >  1:1377731/1‑150 (MQ=255)
                                                                           aTCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGcc                                                                        >  2:460215/1‑150 (MQ=255)
                                                                              ccGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCtt                                                                                                                                        >  1:1544360/1‑83 (MQ=255)
                                                                              ccGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCtt                                                                                                                                        <  2:1544360/83‑1 (MQ=255)
                                                                               cGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGgtg                                                                                                                     >  2:152611/1‑101 (MQ=255)
                                                                               cGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGgtg                                                                                                                     <  1:152611/101‑1 (MQ=255)
                                                                                gCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGc                                                                   <  1:1175943/150‑1 (MQ=255)
                                                                                gCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGc                                                                   <  2:1373829/150‑1 (MQ=255)
                                                                                 cAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGcc                                                                  <  1:92517/150‑1 (MQ=255)
                                                                                   gCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGa                                                                >  1:1136584/1‑150 (MQ=255)
                                                                                    ccAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCa                                                                                                                >  1:135254/1‑101 (MQ=255)
                                                                                    ccAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCa                                                                                                                <  2:135254/101‑1 (MQ=255)
                                                                                      aaCCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCa                                                             >  2:2409888/1‑150 (MQ=255)
                                                                                        cccGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCATGCGTTTTCAg                                                                                                                        >  2:416965/1‑89 (MQ=255)
                                                                                        cccGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCATGCGTTTTCAg                                                                                                                        <  1:416965/89‑1 (MQ=255)
                                                                                         ggggtatatACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATgg                                                          >  1:84409/6‑150 (MQ=255)
                                                                                         ccGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATgg                                                          <  1:989498/150‑1 (MQ=255)
                                                                                          cGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGc                                                         >  1:2039648/1‑150 (MQ=255)
                                                                                              atatACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCTGCAGGCGTGTTCAGTTGTTCGAGGGGCAGCGGCTGGCACTGTCGAAGGCAGGCACCCAGGCCAAAGCCGATCATGGCGaaa                                                     >  2:1706268/1‑148 (MQ=255)
                                                                                               tataCGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGAccc                                                    <  2:1779358/150‑1 (MQ=255)
                                                                                                  aCGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCaaa                                                 <  1:867085/150‑1 (MQ=255)
                                                                                                   cGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAg                                                >  2:1018839/1‑150 (MQ=255)
                                                                                                    ggCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACTTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGc                                               >  1:1836328/1‑150 (MQ=255)
                                                                                                     gCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGcc                                              <  2:2039648/150‑1 (MQ=255)
                                                                                                      cACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGccc                                             <  1:1198138/150‑1 (MQ=255)
                                                                                                      cACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGccc                                             >  2:1285790/1‑150 (MQ=255)
                                                                                                         ggtGCCGGCGATGATGAAGCCGAGACCCTGTACGAAGGCCGCGAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGa                                          >  2:2060123/3‑150 (MQ=255)
                                                                                                         gATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGa                                          >  1:11169/1‑150 (MQ=255)
                                                                                                         gATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGa                                          <  2:1752984/150‑1 (MQ=255)
                                                                                                          aTGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGaa                                         <  2:11169/150‑1 (MQ=255)
                                                                                                                gTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAAGAATTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAgccgcc                                   >  2:330896/1‑150 (MQ=255)
                                                                                                                  gatgatgaAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCcgccgc                                 <  1:1285790/150‑1 (MQ=255)
                                                                                                                         aaGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCCGGTGTTCCAGGGTCGACAACGGTAACGGGCTGGAGCAGGCACCCAAGCCAAATCCGATCATGGCCACCCGAAGCCCCGAAGCCGCCGCAGGCGcc                          >  2:1706136/1‑150 (MQ=255)
                                                                                                                           gCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCCTGGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAg                        <  1:1706268/150‑1 (MQ=255)
                                                                                                                                  cccTGTACGAACGCCGCCAGGCGGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGc                 <  2:1985752/150‑1 (MQ=255)
                                                                                                                                   ccTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGcc                >  2:1388183/1‑150 (MQ=255)
                                                                                                                                   ccTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTACAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGa                                                       <  1:1712108/111‑1 (MQ=255)
                                                                                                                                   ccTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTACAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGa                                                       >  2:1712108/1‑111 (MQ=255)
                                                                                                                                    cTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGccc               <  1:1388183/150‑1 (MQ=255)
                                                                                                                                    cTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGccc               >  2:97539/1‑150 (MQ=255)
                                                                                                                                     ctgcCGAAGCCGGCAGCCTGGCGGCGTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGaa                                         <  1:950298‑M1/119‑1 (MQ=255)
                                                                                                                                     ctgcCGAAGCCGGCAGCCTGGCGGCGTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGaa                                         >  2:950298‑M1/5‑123 (MQ=255)
                                                                                                                                      gTACGAACGCCGCCAGGCTGCCGGATTCGTCAGGGAGTTTCATTTGTTCCAGTGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCCAAGCCCGCAAGCCGCCGCAGGTGCCAGCAAAAGcccct             <  2:1468844/150‑2 (MQ=255)
                                                                                                                                        cCGAAGCCGGCAGCCTGGCGGCGTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGa                                                       <  1:1786754‑M1/105‑1 (MQ=255)
                                                                                                                                        cCGAAGCCGGCAGCCTGGCGGCGTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGa                                                       >  2:1786754‑M1/2‑106 (MQ=255)
                                                                                                                                          GAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGaa         >  2:433156‑M1/1‑150 (MQ=255)
                                                                                                                                          GAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGaa         <  2:2405067‑M1/150‑1 (MQ=255)
                                                                                                                                          GAACGCCGACAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGaa         <  1:2020215‑M1/150‑1 (MQ=255)
                                                                                                                                                 gccAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGAAGCCCGcc  >  1:1246988/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCACGCGGTCACCAGCAGCATTGCCATCACGGTCACGGTCAACAGCGTCCACGAAGCCTGGAAGTCGCCACTGACATCCCGCAGCCAACCCGTGATATACGGCACGATGCCGGTGATGATGAAGCCGATACCCTGTACGAACGCCGCCAGGCTGCCGGCTTCGGCAGGCGTTTTCAGGTGTTCCAGGGTCAACGTCAGGCACTGGCTGAAGCAGGCACCCAAGCCAAAGCCGATCATGGCGACCCAAAGCCCGGAAGCCGCCGCAGGTGCCAGCAAAAGCCCCAGGAAGCCCGCC  >  NC_002947/4280343‑4280637

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