Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,237,247 +C intergenic (+60/‑337) ygjT → / → ygjU putative transport protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,237,2451.C99.0% 274.9 / ‑4.2 103intergenic (+58/‑339)ygjT/ygjUputative transport protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (54/48);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAC  >  NC_000913/3237173‑3237315
                                                                          |                                                                       
cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                       <  1:2793553/76‑1 (MQ=255)
 agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                      <  2:2412133/76‑2 (MQ=255)
 agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                      <  1:537691/76‑2 (MQ=255)
 agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                      >  1:3219633/1‑75 (MQ=255)
   cGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccctt                                                                    <  2:2470364/76‑4 (MQ=255)
      ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCTCCGTTTTCTCCCTCTCCCTTCCa                                                                  <  1:3407960/76‑1 (MQ=255)
      ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                  >  1:238219/1‑76 (MQ=255)
      ggTGGATAATTTTTAATCCGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                  <  1:3165344/76‑1 (MQ=255)
       gTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                <  2:1580608/76‑1 (MQ=255)
        tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                               <  2:3592774/76‑1 (MQ=255)
        tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                               >  2:1615877/1‑76 (MQ=255)
         ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                              >  1:4584162/1‑76 (MQ=255)
         ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                              >  1:3369466/1‑76 (MQ=255)
         ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                              >  2:4482060/1‑76 (MQ=255)
         ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                              >  2:3299158/1‑76 (MQ=255)
          gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                             >  2:406015/1‑76 (MQ=255)
           aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                       <  2:2569904/65‑1 (MQ=255)
            tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                           <  1:4557008/76‑1 (MQ=255)
             aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                          <  1:4627068/76‑1 (MQ=255)
             aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                          <  1:4920755/76‑1 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagagt                                                         <  1:1598257/76‑2 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                         <  2:3221219/76‑1 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                         >  1:319399/1‑76 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                         <  1:1470788/76‑1 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                         <  1:2124652/76‑1 (MQ=255)
              aTTTTTAATCTGCCTAAGCCGTGTACCCTGGCACTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                         >  2:393851/1‑76 (MQ=255)
               tgtttAATCTGCCTCAGCCGTGTGCCCGGTCATTAACATGAGCGCCTTCTTCTCCCTCTCCCTTCCAGGGAGAggg                                                        <  2:155397/74‑1 (MQ=255)
                ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                       <  1:3436644/76‑1 (MQ=255)
                ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGCTTTATCCCTCTCCCTTCCACAGAGAGGGt                                                       >  1:2653529/1‑76 (MQ=255)
                 tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                      <  2:3571411/76‑1 (MQ=255)
                 tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                      <  1:3373337/76‑1 (MQ=255)
                 tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACACGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                      <  2:5188254/76‑1 (MQ=255)
                  ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGCCg                                                     >  2:134405/1‑76 (MQ=38)
                   tAATCTGCCTAAGCCGTGTCCCCTGTCCTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGTAGCGGGTCgg                                                    >  1:3310159/1‑76 (MQ=255)
                   tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                    <  1:1824726/76‑1 (MQ=255)
                    aaTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                   >  2:1931416/1‑76 (MQ=255)
                    aaTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                   <  2:3754070/76‑1 (MQ=255)
                    aaTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                   >  2:4411764/1‑76 (MQ=255)
                     aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                  >  2:560003/1‑76 (MQ=255)
                     aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                  <  1:3072099/76‑1 (MQ=255)
                      tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                 <  2:2965943/76‑1 (MQ=255)
                      tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                 <  2:3879911/76‑1 (MQ=255)
                      tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                 <  1:5652134/76‑1 (MQ=255)
                      tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                                 >  1:2822798/1‑76 (MQ=255)
                       cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCCGGCGAGGGCCcggggg                                                >  2:244278/1‑74 (MQ=255)
                       cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGCCGGGGTg                                                >  2:720258/1‑76 (MQ=25)
                         gCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                              <  1:519968/76‑1 (MQ=38)
                         gCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGGGAg                                              >  2:4420495/1‑76 (MQ=38)
                          ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                             >  2:3161527/1‑76 (MQ=38)
                          ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                             >  2:2168850/1‑76 (MQ=38)
                           cTAAGCCTTGTCCCCGGTCATTAACATGAGCCCCGTTTCCTCCCTCCCCCTTCGGGGGAGGGGGTCGGGGTGAggg                                            <  2:3420005/76‑1 (MQ=255)
                            tAAGCCGTGTACCCTGTCATTGACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                           <  2:1863358/76‑1 (MQ=25)
                            tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgggg                                           >  2:4949445/1‑75 (MQ=37)
                            tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                           >  2:380059/1‑76 (MQ=37)
                            tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                           >  1:96984/1‑76 (MQ=38)
                             aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                          >  1:4586931/1‑76 (MQ=37)
                              aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTaa                                         >  2:285630/1‑76 (MQ=37)
                                ccGTGTACCCTGTCATTAACATGAGCACGGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                       <  1:1332737/76‑1 (MQ=25)
                                ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                       >  1:1476050/1‑76 (MQ=37)
                                ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                       >  2:3748032/1‑76 (MQ=37)
                                 cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                      >  1:2598187/1‑76 (MQ=37)
                                 cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                      >  2:5780698/1‑76 (MQ=25)
                                 cGTGTAACCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                      >  1:405174/1‑76 (MQ=25)
                                 cGTCTACGCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                      <  1:93946/76‑1 (MQ=21)
                                  gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                     >  1:655807/1‑76 (MQ=37)
                                  gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                     >  1:5279145/1‑76 (MQ=37)
                                   tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttgt                                    >  2:3820931/1‑74 (MQ=25)
                                   tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGCGGTGAGGGTAAttttt                                    >  1:94538/1‑76 (MQ=37)
                                    gtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGGGGGGGCGGGGGGGGGGGTATTTTTTc                                   >  1:1005182/1‑76 (MQ=16)
                                    gtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                   >  1:5050606/1‑76 (MQ=37)
                                    gtACCCTGTCAGTAACATGAGCGCCGTTTTCTCCCTCTCCCTTCCAGGGAGCGGGGCGGGGTGAGGGTAATTTTTc                                   <  1:5393306/76‑1 (MQ=21)
                                     gaCCCTTTCCTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCg                                  <  2:148381/75‑1 (MQ=21)
                                      aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGc                                 <  2:4239152/76‑1 (MQ=37)
                                       cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                                >  2:627684/1‑76 (MQ=25)
                                       cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGGGAGGGCAATTTTTCGcc                                >  2:649309/1‑75 (MQ=37)
                                        ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAc                               <  2:4586930/76‑1 (MQ=37)
                                        ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAc                               <  1:2679689/76‑1 (MQ=37)
                                          tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                             >  2:1563020/1‑76 (MQ=25)
                                           gTCATTAGCAGGAGCGCGTTTTTCTCCCTCTCCCGGCGGGGGAGGGGGTGGGGGTGAGGGTAATTTTTCGCACCGa                            <  2:1979482/76‑1 (MQ=255)
                                           gTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGa                            >  1:4070455/1‑76 (MQ=37)
                                            tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                           <  1:891445/76‑1 (MQ=37)
                                               ttAACATGAGCGCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                        <  2:4638529/76‑1 (MQ=25)
                                               ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                        <  2:4591298/76‑1 (MQ=37)
                                                tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                       >  1:4433836/1‑76 (MQ=37)
                                                tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                       >  2:4774780/1‑76 (MQ=25)
                                                 aaCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGCTGCTgg                      >  2:5655137/1‑76 (MQ=25)
                                                   cATGAGCACCTTTTCCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                    <  2:2585605/76‑1 (MQ=25)
                                                   cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                    <  1:1683232/76‑1 (MQ=37)
                                                   cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAACGATGCTGGcc                    >  1:107211/1‑76 (MQ=25)
                                                      gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                 <  2:291588/76‑1 (MQ=37)
                                                       aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                <  2:4288690/76‑1 (MQ=25)
                                                          aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc             >  2:887574/1‑76 (MQ=25)
                                                             gTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTc          >  2:4759917/1‑76 (MQ=25)
                                                              ttttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCa         <  2:5508773/76‑1 (MQ=21)
                                                              ttttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCa         >  1:3871058/1‑76 (MQ=21)
                                                                  ctccctctcccTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCt     >  2:3570786/4‑76 (MQ=25)
                                                                   tccctctcccTTCCAGGGAGAGGGCCGGGGTGAGGGTAATTTTTCGCACCGCTGCTGGCCTGTTCCCCTCCCCCTa    >  2:580536/3‑76 (MQ=17)
                                                                    ccctctcccTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa   >  1:3123956/2‑76 (MQ=255)
                                                                     cctctcccTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGCTCCCCTCACCCTAAc  >  2:5074992/1‑76 (MQ=25)
                                                                     cctctcccTTCCAGGGAGAGGGTCGGGGTCAGGGTAATTTTTCGCACCGATGCTGGCCTGCTCCCCCCACCCCAAc  >  2:729690/1‑76 (MQ=21)
                                                                          |                                                                       
CAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAC  >  NC_000913/3237173‑3237315

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