Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,510,368 T→G 10.1% G203G (GGA→GGC yhfX ← putative PLP‑binding protein YhfX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,3680TG10.1% 165.4 / 7.4 58G203G (GGA→GGCyhfXputative PLP‑binding protein YhfX
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base G (4/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATC  >  NC_000913/3510227‑3510508
                                                                                                                                             |                                                                                                                                            
cacaCCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGAAAAGTTTTGCCGACacc                                                                                                                                      <  1:977974/150‑3 (MQ=255)
  cacCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCt                                                                                                                                    <  2:1296969/150‑1 (MQ=255)
  cacCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCAGAAAGcca                                                                                                                                    >  2:174213/1‑149 (MQ=255)
   acCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGAAAAGTTTTGCCGACaccgaa                                                                                                                                   <  1:455188/150‑6 (MQ=255)
     cGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCAt                                                                                                                                 <  1:1224503/150‑1 (MQ=255)
      gTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATc                                                                                                                                >  1:480490/1‑150 (MQ=255)
        aTTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATccc                                                                                                                              <  1:508107/150‑1 (MQ=255)
           gcgcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCaaa                                                                                                                           >  2:1342081/1‑150 (MQ=255)
           gcgcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCaaa                                                                                                                           <  1:967452/150‑1 (MQ=255)
             gcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGc                                                                                                                         <  1:242638/150‑1 (MQ=255)
               cAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAg                                                                                                                       >  1:596887/1‑150 (MQ=255)
                aGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCATGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAgg                                                                                                                      <  1:439115/150‑1 (MQ=255)
                  tAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCGGGCa                                                                                                                    <  1:240978/150‑1 (MQ=255)
                       gCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGa                                                                                                               >  1:1298170/1‑150 (MQ=255)
                         aGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACCGCCTCATCCCAAAGCAGGCAAGGGaaa                                                                                                             >  1:1121184/1‑150 (MQ=255)
                                      tGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCgg                                                                                                >  1:228620/1‑150 (MQ=255)
                                        gTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCt                                                                                              <  2:1298170/150‑1 (MQ=255)
                                         tCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCCGTGTCGCCAAAACTATTAAAACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTa                                                                                             >  2:1218848/1‑150 (MQ=255)
                                                aGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAg                                                                                      <  1:861636/150‑1 (MQ=255)
                                                          aGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGa                                                                            >  2:1148383/1‑150 (MQ=255)
                                                                tCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTgg                                                                      <  1:954937/150‑1 (MQ=255)
                                                                 cAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGa                                                                     <  1:900300/150‑1 (MQ=255)
                                                                 cAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTACGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGa                                                                     >  2:642040/1‑150 (MQ=255)
                                                                   aGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGAtt                                                                   <  1:1342081/150‑1 (MQ=255)
                                                                     tGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTc                                                                 <  1:897160/150‑1 (MQ=255)
                                                                        aaTACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGc                                                              <  2:272218/150‑1 (MQ=255)
                                                                         aTACCAGATTTCGCCAGTTGATCCCGTGTCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCg                                                             <  2:315903/150‑1 (MQ=255)
                                                                         aTACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTAGGTTTCTGCAAAACTATTAAAACATACTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCg                                                             >  2:966317/1‑150 (MQ=255)
                                                                           aCCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGAc                                                           <  1:323648/150‑1 (MQ=255)
                                                                            ccAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATAACGTTTCTCCACAAATAATAAAAAAAAATAATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACa                                                          >  2:408625/1‑150 (MQ=255)
                                                                                    cGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCgg                                                  >  1:1330294/1‑150 (MQ=255)
                                                                                      ccAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCa                                                <  1:994337/150‑1 (MQ=255)
                                                                                          ttGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAacac                                            <  1:1216408/150‑1 (MQ=255)
                                                                                           tGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACAcc                                           >  1:25234/1‑150 (MQ=255)
                                                                                            gATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  1:743707/1‑150 (MQ=255)
                                                                                            gATCCCGTGCCTGTATCAGCGTGTGAAGATTAGGTGTCGGCAAAACTTTTCAGAAAGCCTCATCCCAAAGCAGGTAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  2:550153/1‑150 (MQ=255)
                                                                                            gATCCCGTGCCTGTATCAGCGTGTGAAGATTACGTGTCTGCAAAACTATTCCGAAAGCATCATCCCAAAGCAGGTAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  2:554934/1‑150 (MQ=255)
                                                                                              tCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTg                                        <  2:645991/150‑1 (MQ=255)
                                                                                                ccGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAg                                      >  1:672758/1‑150 (MQ=255)
                                                                                                  gTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAggg                                    <  1:66456/150‑1 (MQ=255)
                                                                                                  gTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAggg                                    <  1:1275162/150‑1 (MQ=255)
                                                                                                          aTCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAACAGGCAAGGGAAATGGGTAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGc                           >  1:888658/1‑150 (MQ=255)
                                                                                                             aGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGctc                         <  2:596887/150‑1 (MQ=255)
                                                                                                             aGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGctc                         <  1:1148383/150‑1 (MQ=255)
                                                                                                                 tgtgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTGGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCt                     >  1:1290671/1‑150 (MQ=255)
                                                                                                                 tgtgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCt                     <  1:800483/150‑1 (MQ=255)
                                                                                                                   tgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTgg                   >  1:266585/1‑150 (MQ=255)
                                                                                                                    gAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGc                  >  1:330923/1‑150 (MQ=255)
                                                                                                                        aTTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCgg              <  2:672758/150‑1 (MQ=255)
                                                                                                                         gaggctgtCGGAAAAGTTTTGCCGACAGTCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGa             >  2:937716/6‑150 (MQ=255)
                                                                                                                           cGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATa           >  2:285142/1‑150 (MQ=255)
                                                                                                                            gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaa          >  1:186165/3‑150 (MQ=255)
                                                                                                                            gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaa          >  2:829643/3‑150 (MQ=255)
                                                                                                                            gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAAACTCCGGCAACACCTTCAGGGCAAAAACGCTCGCCGGGCCCGGAGaa          >  2:1091796/3‑150 (MQ=255)
                                                                                                                             gtgtCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaaa         <  1:642040/150‑1 (MQ=255)
                                                                                                                             gtgtCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaaa         <  1:608542/150‑1 (MQ=255)
                                                                                                                                   gCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAat   >  1:1149877/1‑150 (MQ=255)
                                                                                                                                    cAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAatc  >  2:981308/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATC  >  NC_000913/3510227‑3510508

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