Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,704,719 G→T P15T (CCG→ACG)  dppC ← dipeptide/heme ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,704,7190GT96.4% 290.2 / ‑0.8 84P15T (CCG→ACG) dppCdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (34/47);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGATGAACAGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAG  >  NC_000913/3704623‑3704815
                                                                                                |                                                                                                
tCGCGATGAACAGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCAc                                                                                               <  1:473112/100‑1 (MQ=255)
      tGAACAGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGCTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTgc                                                                                         <  2:16520/100‑1 (MQ=255)
          cAGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGGTGCGTTTAAAATAGTGCCAGAACTCCGGTAACGGGGTCATCGTCACCGGTGCGCTa                                                                                     <  2:49244/100‑1 (MQ=255)
           aGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTaa                                                                                    >  2:798876/1‑100 (MQ=255)
           aGCACGATGACGACGTAAACCAGCCCGACGACCGCGCATTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTaa                                                                                    >  2:1101260/1‑100 (MQ=255)
            gCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAAt                                                                                   >  2:786608/1‑100 (MQ=255)
            gCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAAt                                                                                   <  2:344211/100‑1 (MQ=255)
             cACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATc                                                                                  <  2:1996048/100‑1 (MQ=255)
             cACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATc                                                                                  <  1:1506642/100‑1 (MQ=255)
                 aTGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACt                                                                               >  2:1795186/1‑99 (MQ=255)
                      gacgTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTAttt                                                                         <  1:269697/100‑1 (MQ=255)
                      gacgTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTAttt                                                                         >  1:1469065/1‑100 (MQ=255)
                       acgTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTAtttt                                                                        <  2:1290227/100‑1 (MQ=255)
                        cgTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTc                                                                       >  1:355212/1‑100 (MQ=255)
                          tAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAg                                                                     <  2:1987305/100‑1 (MQ=255)
                          tAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAg                                                                     <  2:930596/100‑1 (MQ=255)
                           aaaCCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGt                                                                    >  1:915810/1‑100 (MQ=255)
                           aaaCCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGt                                                                    <  1:222363/100‑1 (MQ=255)
                           aaaCCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGt                                                                    >  1:1451408/1‑100 (MQ=255)
                            aaCCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTa                                                                   <  2:1607256/100‑1 (MQ=255)
                             aCCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTaa                                                                  >  2:1610051/1‑100 (MQ=255)
                              ccAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAAc                                                                 >  1:827850/1‑100 (MQ=255)
                              ccAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAAc                                                                 >  1:243276/1‑100 (MQ=255)
                               cAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAAcc                                                                <  2:1427278/100‑1 (MQ=255)
                                 gCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCtg                                                              >  1:1680585/1‑100 (MQ=255)
                                   ccGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCtgtg                                                            >  2:1041669/1‑100 (MQ=255)
                                         cgatctccTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatga                                                      <  2:921065/94‑1 (MQ=255)
                                         acCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatga                                                      >  2:4482/1‑100 (MQ=255)
                                         acCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatga                                                      >  2:84543/1‑100 (MQ=255)
                                          ccGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatga                                                      <  2:1413327/99‑1 (MQ=255)
                                           cgcgCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatgatg                                                    <  1:547640/100‑1 (MQ=255)
                                           cgcgCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatgatg                                                    <  2:465982/100‑1 (MQ=255)
                                              gcCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTcc                                                 <  2:1620138/100‑1 (MQ=255)
                                               ccTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACGTGTGACATGATGTccc                                                >  2:1199953/1‑100 (MQ=255)
                                               ccTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACatgaagatcg                                                >  1:921069/1‑94 (MQ=255)
                                                cTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTcccc                                               <  2:181136/100‑1 (MQ=255)
                                                cTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTcccc                                               <  1:1394940/100‑1 (MQ=255)
                                                 tttGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCt                                              <  1:990890/100‑1 (MQ=255)
                                                 tttGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCt                                              <  2:1469056/100‑1 (MQ=255)
                                                  ttgttgCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCtt                                             <  2:915806/100‑1 (MQ=255)
                                                   tgttgCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTa                                            >  1:1259498/1‑100 (MQ=255)
                                                    gttgCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTAc                                           >  2:938612/1‑100 (MQ=255)
                                                    gttgCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTAc                                           <  2:98306/100‑1 (MQ=255)
                                                        cGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTActtct                                       <  1:1022382/100‑1 (MQ=255)
                                                          tttAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTa                                     <  2:1494417/100‑1 (MQ=255)
                                                          tttAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTa                                     <  2:644631/100‑1 (MQ=255)
                                                          atctaaaTAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTa                                     <  2:1371288/96‑1 (MQ=255)
                                                             aaaaTAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGa                                  <  1:798879/100‑1 (MQ=255)
                                                              ataTAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTaagat                                 >  1:1371297/3‑96 (MQ=255)
                                                                aTAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGa                               >  2:975287/1‑100 (MQ=255)
                                                                 tAGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGaa                              <  1:1041674/100‑1 (MQ=255)
                                                                  attGCCAGAACGCCTGTAACGGTGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACg                                <  2:1959736/95‑1 (MQ=255)
                                                                  aGTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAAt                             <  2:53007/100‑1 (MQ=255)
                                                                   gTGCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGaaga                            >  2:1241412/1‑98 (MQ=255)
                                                                     gCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcg                          >  2:395774/1‑100 (MQ=255)
                                                                     gCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcg                          >  1:36573/1‑100 (MQ=255)
                                                                     gCCAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcg                          >  1:580113/1‑100 (MQ=255)
                                                                      ccAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGCCCCCTTACTTCTTATGACGAATAcgc                         <  1:1241418/100‑1 (MQ=255)
                                                                       cAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcgcg                        >  2:4712/1‑100 (MQ=255)
                                                                       cAGAACTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcgcg                        <  2:355211/100‑1 (MQ=255)
                                                                       cAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcgcg                        <  2:728142/100‑1 (MQ=255)
                                                                       cAGAACTCCTGTAACGGGCTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcgcg                        >  2:143564/1‑100 (MQ=255)
                                                                            cTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTc                   >  1:1158487/1‑100 (MQ=255)
                                                                            cTCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTc                   <  1:866838/100‑1 (MQ=255)
                                                                             tCCTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTca                  <  2:382850/100‑1 (MQ=255)
                                                                              ccTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcac                 >  2:1538065/1‑100 (MQ=255)
                                                                              ccTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcac                 >  2:473982/1‑100 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACGGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  1:169251/100‑1 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  2:404630/100‑1 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  1:1698146/100‑1 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  1:395775/100‑1 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  2:1896254/100‑1 (MQ=255)
                                                                               cTGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacc                <  1:84544/100‑1 (MQ=255)
                                                                                tGTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcacca               <  1:1795202/100‑1 (MQ=255)
                                                                                 gTAACGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTcaccac              <  2:1324031/100‑1 (MQ=255)
                                                                                    aCGGGGTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGcc           >  1:134180/1‑100 (MQ=255)
                                                                                      ggggTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCAcgtcgt         >  2:856853/1‑100 (MQ=255)
                                                                                      ggggTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGt         <  2:457377/100‑1 (MQ=255)
                                                                                        ggTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTAc       <  2:1680572/100‑1 (MQ=255)
                                                                                        ggTCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTAc       >  2:1224552/1‑100 (MQ=255)
                                                                                          tCATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTAcag     <  1:1546077/100‑1 (MQ=255)
                                                                                           cATCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTAcagc    >  2:774301/1‑100 (MQ=255)
                                                                                             tCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTAcagcag  >  1:101214/1‑100 (MQ=255)
                                                                                             tCGTCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTGATGACGAATACGCGGGGTCACCACGCCGTAcagcgg  >  1:1630910/1‑98 (MQ=255)
                                                                                                |                                                                                                
TCGCGATGAACAGCACGATGACGACGTAAACCAGCCCGACGACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAG  >  NC_000913/3704623‑3704815

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