Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,777,373 A→C 5.3% K545T (AAG→ACG)  betA‑II → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,777,3730AC 5.3% 275.7 / 4.0 93K545T (AAG→ACG) betA‑IIcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (43/45);  new base C (3/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATC  >  NC_002947/5777225‑5777518
                                                                                                                                                    |                                                                                                                                                 
gTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGACCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAg                                                                                                                                                  <  1:189128/150‑1 (MQ=255)
    aCGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGtact                                                                                                                                              >  1:1717865/1‑150 (MQ=255)
              gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACTACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGAAAGCCACTGCCGCGAAGCAGCGCCAATTACTGCGTGGTggg                                                                                                                                    >  2:884921/1‑150 (MQ=255)
              gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCggg                                                                                                                                    <  1:1567566/150‑1 (MQ=255)
                   tGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATg                                                                                                                               >  2:1394910/1‑150 (MQ=255)
                      gcgTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCAc                                                                                                                            <  2:884722/150‑1 (MQ=255)
                          ggTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGt                                                                                                                        <  2:1061694/150‑1 (MQ=255)
                              taTGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGaaa                                                                                                                    >  2:747171/2‑150 (MQ=255)
                              gATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGaaa                                                                                                                    <  1:1394910/150‑1 (MQ=255)
                               aTGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGTACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAg                                                                                                                   >  2:82915/1‑150 (MQ=255)
                                   cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCCCCGCCACGTAGTACTTGGCggt                                                                                                                                    >  1:1415215/1‑128 (MQ=255)
                                   cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCCCCGCCACGTAGTACTTGGCggt                                                                                                                                    <  2:1415215/129‑2 (MQ=255)
                                   cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCgg                                                                                                                                     <  1:1435229/128‑1 (MQ=255)
                                   cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCgg                                                                                                                                     >  2:1435229/1‑128 (MQ=255)
                                    tCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGt                                                                                                                                                 >  1:1186353/1‑115 (MQ=255)
                                    tCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGt                                                                                                                                                 <  2:1186353/115‑1 (MQ=255)
                                     cGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGc                                                                                                             <  2:1134923/150‑1 (MQ=255)
                                      gATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGcc                                                                                                            >  1:783799/1‑150 (MQ=255)
                                        tCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGa                                                                                                          <  2:1024457/150‑1 (MQ=255)
                                        tCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGa                                                                                                          >  2:1285615/1‑150 (MQ=255)
                                              cGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtg                                                                                                    >  1:684098/1‑150 (MQ=255)
                                               gCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtga                                                                                                   >  1:1595235/1‑150 (MQ=255)
                                                cTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGGGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACGGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                  >  2:43171/1‑150 (MQ=255)
                                                  gATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGTACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaagt                                                                                                <  1:82915/150‑1 (MQ=255)
                                                  gATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCggg                                                                                                                                    >  2:1062734/1‑114 (MQ=255)
                                                  gATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCggg                                                                                                                                    <  1:1062734/114‑1 (MQ=255)
                                                     catcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaagtgaa                                                                                             >  2:745537/1‑150 (MQ=255)
                                                       tcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAAc                                                                                           <  2:1621725/150‑1 (MQ=255)
                                                       tcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACAGCCGCGTACTACGTAGCGGGGGCTGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAAc                                                                                           >  1:961022/1‑150 (MQ=255)
                                                          ccGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAgg                                                                                        >  2:1902378/1‑150 (MQ=255)
                                                           cGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAggg                                                                                       <  2:1689802/150‑1 (MQ=255)
                                                             gCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCt                                                                                     <  2:437770/150‑1 (MQ=255)
                                                               aaCCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGa                                                                                   >  2:149761/1‑150 (MQ=255)
                                                               aaCCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGa                                                                                   <  1:64929/150‑1 (MQ=255)
                                                                aCCTCAACGCCACCACGATCATGATAGACGAGAAGAGTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGGGGAGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAg                                                                                  <  2:778326/150‑1 (MQ=255)
                                                                 ccTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAgg                                                                                 <  1:1285615/150‑1 (MQ=255)
                                                                  cTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGta                                                                                                                                                >  1:1790675/1‑86 (MQ=255)
                                                                  cTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGta                                                                                                                                                <  2:1790675/86‑1 (MQ=255)
                                                                    cAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCt                                                                              >  1:712972/1‑150 (MQ=255)
                                                                     aaCGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCtt                                                                             <  2:1889029/150‑1 (MQ=255)
                                                                      aCGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTg                                                                            >  1:59441/1‑150 (MQ=255)
                                                                      aCGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTg                                                                            <  1:149761/150‑1 (MQ=255)
                                                                        gCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAg                                                                                                                                                  <  2:555085/78‑1 (MQ=255)
                                                                        gCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAg                                                                                                                                                  >  1:555085/1‑78 (MQ=255)
                                                                          caccacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGGAGTGAAACAGGGCTGCGGCCTTGAtaa                                                                        >  2:339867/1‑148 (MQ=255)
                                                                           accacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAGCAGGGCTGAGGCCTTGATCAg                                                                       <  2:1595235/150‑1 (MQ=255)
                                                                           accacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAg                                                                       >  2:936809/1‑150 (MQ=255)
                                                                            ccacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAgg                                                                      <  1:154553/150‑1 (MQ=255)
                                                                            ccacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAgg                                                                      >  1:794076/1‑150 (MQ=255)
                                                                             cacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAggg                                                                     <  2:1843914/150‑1 (MQ=255)
                                                                             cacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAggg                                                                     >  2:395746/1‑150 (MQ=255)
                                                                                  tCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAg                                                                                                                   >  1:1297812/1‑99 (MQ=255)
                                                                                  tCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAg                                                                                                                   <  2:1297812/99‑1 (MQ=255)
                                                                                     tGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGACGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                  >  1:1124716/1‑113 (MQ=255)
                                                                                     tGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGACGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                  <  2:1124716/113‑1 (MQ=255)
                                                                                       aTCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGc                                                           <  2:783799/150‑1 (MQ=255)
                                                                                        tCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCg                                                          >  2:606817/1‑150 (MQ=255)
                                                                                            cGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGccc                                                      <  1:747171/150‑1 (MQ=255)
                                                                                             gagaAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGcccc                                                     >  1:82883/1‑150 (MQ=255)
                                                                                               gaagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCtt                                                   >  1:1041720/1‑150 (MQ=255)
                                                                                                aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCttt                                                  <  2:1717865/150‑1 (MQ=255)
                                                                                                aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCttt                                                  >  2:1498472/1‑150 (MQ=255)
                                                                                                aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCttt                                                  <  2:951557/150‑1 (MQ=255)
                                                                                                   aTTTCGGACAAGATCCTGGGTCGCACGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCAc                                               <  2:712972/150‑1 (MQ=255)
                                                                                                   aTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCAc                                               <  1:1498472/150‑1 (MQ=255)
                                                                                                    tttCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACg                                              >  1:1143696/1‑150 (MQ=255)
                                                                                                     ttCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGt                                             <  1:1902378/150‑1 (MQ=255)
                                                                                                        ggACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTg                                          >  2:913154/1‑150 (MQ=255)
                                                                                                              gATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAggg                                                                                       >  1:902317/1‑99 (MQ=255)
                                                                                                              gATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAggg                                                                                       <  2:902317/99‑1 (MQ=255)
                                                                                                               aTCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGc                                   <  1:1446457/150‑1 (MQ=255)
                                                                                                                   gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGc                               >  2:845444/1‑150 (MQ=255)
                                                                                                                   gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGc                               >  1:756109/1‑150 (MQ=255)
                                                                                                                   gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGc                               >  2:646972/1‑150 (MQ=255)
                                                                                                                       tCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCCTTCACGTTTGCTCAGCCGGGCAgcc                           >  2:6498/1‑149 (MQ=255)
                                                                                                                         gCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCa                         >  2:553168/1‑150 (MQ=255)
                                                                                                                             gCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCaaaaa                     >  2:984483/1‑150 (MQ=255)
                                                                                                                              ccACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCaaaaaa                    >  1:1669324/1‑150 (MQ=255)
                                                                                                                               cACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCCGGGGTCGCGCAGCGGCCCCTTTCAAGTTTGCTTAGTCGGGCAGCGCAAAAAAt                   <  1:43171/150‑1 (MQ=255)
                                                                                                                                aCTGCCGCGCAGCACCGACAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtg                  <  2:59441/150‑1 (MQ=255)
                                                                                                                                aCTGCCGCGCACCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGCGCTGAGGCCTTGATCAGGGCCCGCGCAGCGGCGCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtg                  <  1:884921/150‑1 (MQ=255)
                                                                                                                                   gCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc               <  1:913154/150‑1 (MQ=255)
                                                                                                                                        gcagcaCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATc          <  1:553168/150‑1 (MQ=255)
                                                                                                                                         atcccccGCCACGTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCt                                                    <  1:824710/102‑1 (MQ=255)
                                                                                                                                         atcccccGCCACGTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCt                                                    >  2:824710/6‑107 (MQ=255)
                                                                                                                                            caCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACg      >  2:1701251/1‑150 (MQ=255)
                                                                                                                                              ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                  <  2:659374/56‑1 (MQ=255)
                                                                                                                                              ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaa                                                                                                  >  1:659374/1‑56 (MQ=255)
                                                                                                                                              ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa    <  1:606817/150‑1 (MQ=255)
                                                                                                                                               cGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCgg                                 >  2:731404/1‑120 (MQ=255)
                                                                                                                                               cGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCgg                                 <  1:731404/120‑1 (MQ=255)
                                                                                                                                                gCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAatc  <  2:1143696/150‑1 (MQ=255)
                                                                                                                                                gCCAAGTACTACGTAGCGGGGACTTCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGACCCTTTCACGTTGGCTCATGCGGGCAGGGCAAAAAATGTGCGCATCAACGAatc  <  2:961022/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATC  >  NC_002947/5777225‑5777518

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: