Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,369,596 G→A 6.0% D197D (GAC→GAT ramA ← (R)‑stereoselective amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,369,5960GA 6.0% 200.7 / 2.1 66D197D (GAC→GATramA(R)‑stereoselective amidase
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base A (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGTAGTCATATTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGTTGGCCGAGCGCCCGAGCA  >  NC_002947/4369448‑4369743
                                                                                                                                                    |                                                                                                                                                    
ggTAGTCATATTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTcgtcgccg                                                                                                                                                      >  5:262718/1‑149 (MQ=255)
  tAGTCATATTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTcgtcgccgtc                                                                                                                                                    <  5:96876/149‑1 (MQ=255)
     tCATATTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCgtcgccgtcgcc                                                                                                                                                 >  6:24032/1‑149 (MQ=255)
       atatTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCa                                                                                                                                               >  3:6473/1‑149 (MQ=255)
        tatTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCAc                                                                                                                                              >  3:36785/1‑149 (MQ=255)
        tatTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCAc                                                                                                                                              <  5:370185/149‑1 (MQ=255)
             gTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACCGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCCCCCGag                                                                                                                                         >  2:105094/1‑148 (MQ=255)
             gTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGAt                                                                                                                                         >  8:162657/1‑149 (MQ=255)
             gTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGAt                                                                                                                                         >  2:372345/1‑149 (MQ=255)
                gcgggcggTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCgtcgccgtcgcc                                                                                                                                                 >  5:368420/1‑138 (MQ=255)
                gcgggcggTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCgtcgccgtcgcc                                                                                                                                                 <  6:368420/138‑1 (MQ=255)
                gcgggcggTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGa                                                                                                                                      <  4:6473/149‑1 (MQ=255)
                    gcggTGCGCAGCAGGTCGAGATCTAATTCCACCACCTACCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCAt                                                                                                                                  >  1:257973/1‑149 (MQ=255)
                        tGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCAcc                                                                                                                              <  1:372345/149‑1 (MQ=255)
                                  gTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCt                                                                                                                   <  4:36785/149‑1 (MQ=255)
                                  gTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCt                                                                                                                   >  3:337345/1‑149 (MQ=255)
                                   tCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTg                                                                                                                  <  1:234372/149‑1 (MQ=255)
                                    cGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTgg                                                                                                                 >  2:128614/1‑149 (MQ=255)
                                      agaTCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGtt                                                                                                               >  4:19201/1‑149 (MQ=255)
                                       gaTCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGttt                                                                                                              >  8:326616/1‑149 (MQ=255)
                                         gcTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTc                                                                                                            <  2:99643/148‑1 (MQ=255)
                                          cTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTcc                                                                                                           >  3:332394/1‑149 (MQ=255)
                                             aTTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTcc                                                                                                           <  8:301393/146‑1 (MQ=255)
                                             aTTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTc                                                                                                            <  2:257973/145‑1 (MQ=255)
                                             aTTCCACCACCTGCCGACACTCCGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGTGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            <  6:163423/130‑1 (MQ=255)
                                             aTTCCACCACCTGCCGACACTCAGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGTGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            >  5:163423/1‑130 (MQ=255)
                                              ttCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATc                                                                                                                                 <  7:79796/124‑1 (MQ=255)
                                              ttCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATc                                                                                                                                 >  8:79796/1‑124 (MQ=255)
                                                    caccTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAg                                                                                                 >  6:82703/1‑149 (MQ=255)
                                                             acacTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGtt                                                                                                               <  8:204061/126‑1 (MQ=255)
                                                             acacTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGtt                                                                                                               >  7:204061/1‑126 (MQ=255)
                                                             acacTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTg                                                                                        >  8:354910/1‑149 (MQ=255)
                                                             acacTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTAGTTTTCCTGAGCGCGAGCCATTATTg                                                                                        <  5:24032/149‑1 (MQ=255)
                                                             acacTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTAGCCGTCGCCCACCCGATTGACCATAACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCAGTATTg                                                                                        >  4:224501/1‑149 (MQ=255)
                                                                 tctGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGt                                                                                    >  2:131306/1‑149 (MQ=255)
                                                                     cgcgGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCATCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGa                                                                                                        <  2:225848/125‑1 (MQ=255)
                                                                     cgcgGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCATCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGa                                                                                                        >  1:225848/1‑125 (MQ=255)
                                                                      gcgGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            >  4:84062/1‑105 (MQ=255)
                                                                      gcgGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            <  3:84062/105‑1 (MQ=255)
                                                                             gCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAgg                                                                        <  6:262718/149‑1 (MQ=255)
                                                                              cTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCAcc                                                                                                                                             <  8:212519/80‑1 (MQ=255)
                                                                              cTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCAcc                                                                                                                                             >  7:212519/1‑80 (MQ=255)
                                                                              cTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAggg                                                                       <  3:353870/149‑1 (MQ=255)
                                                                              cTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAggg                                                                       <  4:332394/149‑1 (MQ=255)
                                                                              cTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAggg                                                                       <  3:131146/149‑1 (MQ=255)
                                                                                  gAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCg                                                                    >  8:51926/1‑148 (MQ=255)
                                                                                   aaTAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGt                                                                   >  8:34036/1‑148 (MQ=255)
                                                                                        cGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGAt                                                             >  3:67096/1‑149 (MQ=255)
                                                                                          tGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGagc                                                            >  7:157776/1‑146 (MQ=255)
                                                                                            cgcCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAATCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCAt                                                         <  7:162657/149‑1 (MQ=255)
                                                                                               ccGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGt                                                       >  5:87675/1‑148 (MQ=255)
                                                                                               ccGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCATCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGt                                                       >  8:233302/1‑148 (MQ=255)
                                                                                                cGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTa                                                     >  7:247117/1‑149 (MQ=255)
                                                                                                cGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGACGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGGCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTa                                                     >  5:256584/1‑149 (MQ=255)
                                                                                                     ggATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGtt                                                <  5:82703/149‑1 (MQ=255)
                                                                                                          gACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGa                                           <  7:326616/149‑1 (MQ=255)
                                                                                                           aCAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATcc                                                           <  4:172582/132‑1 (MQ=255)
                                                                                                           aCAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATcc                                                           >  3:172582/1‑132 (MQ=255)
                                                                                                           aCAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGAc                                          <  7:51926/149‑1 (MQ=255)
                                                                                                             aGCGGCGCTGCCGCCGGCAAAAACGACCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCTTTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCa                                        <  1:84118/149‑1 (MQ=255)
                                                                                                                 gcgcTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAatat                                    <  7:354910/149‑1 (MQ=255)
                                                                                                                 gcgcTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCATCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAatat                                    <  7:233302/149‑1 (MQ=255)
                                                                                                                            ggCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGt                                                 <  6:302758/125‑1 (MQ=255)
                                                                                                                            ggCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGt                                                 >  5:302758/1‑125 (MQ=255)
                                                                                                                               aaaaaCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            <  6:31662/48‑1 (MQ=255)
                                                                                                                               aaaaaCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACcgc                                                                                                                            >  5:31662/1‑48 (MQ=255)
                                                                                                                                 aaaCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGt                    >  2:12980/1‑149 (MQ=255)
                                                                                                                                 aaaCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGt                    >  5:318486/1‑149 (MQ=255)
                                                                                                                                   aCGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAatat                                    >  1:81600/1‑131 (MQ=255)
                                                                                                                                      aGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGTTGGcc               <  1:105094/149‑1 (MQ=255)
                                                                                                                                          cTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCAAGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGTTGGCCGAg            >  1:269083/1‑149 (MQ=255)
                                                                                                                                             cgtcgccgtcgcCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGTTGGCCGAgcgc         <  4:67096/148‑1 (MQ=255)
                                                                                                                                                    gtcgccCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTCGACTGCGGTGCACAGGGCCGTAACGATCCATGTTACCCTTAGTGACAAATATGAGCTCGCCGCCGAGTTGGCCCAGCGCCCGAGCa  >  6:275585/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTAGTCATATTCGTGGCGGGCGGTGCGCAGCAGGTCGAGATCTAATTCCACCACCTGCCGACACTCTGCGCGGCCCGCTTCGAATAGCGTGCGCCCGAAGGGATCGACAGCGGCGCTGCCGCCGGCAAAAACGAGCCCTTCGTCGCCGTCGCCCACCCGATTGACCATCACCGCGAATGCCTGGTTTTCCTGAGCGCGAGCCATTATTGCAGTGCGGTGGACAGGGCCGTAAGGATCCATGTTACCGTTAGTGACCAATATGAGCTCGGCGCCGAGTTGGCCGAGCGCCCGAGCA  >  NC_002947/4369448‑4369743

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