Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,378,146 G→C 9.1% L372V (CTG→GTG) ‡ menD ← 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,378,1460GC 9.1% 230.7 / 10.7 74L372V (CTG→GTG) ‡menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base C (5/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CGATACCGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTCCGC  >  NC_000913/2378005‑2378295
                                                                                                                                             |                                                                                                                                                     
cGATACCGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCtt                                                                                                                                               >  1:840541/1‑150 (MQ=255)
   tACCGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGc                                                                                                                                            <  1:415633/150‑1 (MQ=255)
    aCCGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGAAGCTCAACTGGCGCATCGca                                                                                                                                           <  2:1312882/150‑2 (MQ=255)
     ccGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCGTCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCa                                                                                                                                          >  1:870539/1‑150 (MQ=255)
      cGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGAAGCGCAACTGGCGCATCGcatc                                                                                                                                         <  1:1311724/150‑4 (MQ=255)
        cTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAAc                                                                                                                                       <  1:465941/150‑1 (MQ=255)
         tGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACg                                                                                                                                      <  2:840541/150‑1 (MQ=255)
            ccccACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCAt                                                                                                                                   >  2:1032010/1‑150 (MQ=255)
            ccccACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCAt                                                                                                                                   <  1:134730/150‑1 (MQ=255)
                     tGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTATCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCggg                                                                                                                          <  2:1352817/150‑1 (MQ=255)
                     tGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGTCAGATAGTCGCAGATGCGATGCGCCAGTTGCTCTTCGCCAAACGCATCAGGGCggg                                                                                                                          <  2:867119/150‑1 (MQ=255)
                      gCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACTAACAAAAAATTGCCCCTTTTCAGTCAGATAGTCTCAGATGCGATGCGCCAGTTGGGCTTCGCCAAACGCATCACGGCGTGc                                                                                                                         <  2:708162/150‑1 (MQ=255)
                       cTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCa                                                                                                                        >  1:1227186/1‑150 (MQ=255)
                       cTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCa                                                                                                                        >  1:601104/1‑150 (MQ=255)
                            cacCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAAc                                                                                                                   >  1:1074851/1‑150 (MQ=255)
                               ctggTAACTTGCGGTATGTTGCTAAAGTGCATAATTCCCTCTTACCACAAGGTTTTTACCAACAAACAATTGACCTGTTCAGGCAGATAGTCGGAGATGCGATGCGCCAGTTGCGCTTCGCCAAACTCATCACGGCGTGCAATAACCGcc                                                                                                               <  2:1478520/148‑1 (MQ=255)
                                gggTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGTAGATATTCGTAGATGCTATGCGCCAGTTTCGCTTAGCCAAACGCATCACGGCGGGCAATAAcctgc                                                                                                               <  2:708760/150‑4 (MQ=255)
                                         gCCGTAAGTTGCCAAAGCGCATCAATCACACGTACCACCCGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTGGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGcc                                                                                                      <  1:1168610/150‑1 (MQ=255)
                                         gCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTTCCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGcc                                                                                                      <  1:767976/150‑1 (MQ=255)
                                         gCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGcc                                                                                                      >  1:29422/1‑150 (MQ=255)
                                         gCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGcc                                                                                                      <  2:601104/150‑1 (MQ=255)
                                         gCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGcc                                                                                                      <  2:659524/150‑1 (MQ=255)
                                            ggAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGt                                                                                                   >  2:363764/1‑150 (MQ=255)
                                                          cccATCAATCAGAAGTACAACCACGCTGTTACCAACAAACAATAGCCCCTGTTCAGGAACATAGTCGCAGATGCGATGCGCCACTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCgg                                                                                     >  2:808118/3‑150 (MQ=255)
                                                           gcATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGGATCACGGCGGGCAATAACCGCCTGCATTGCGTGTTCCGCCATGCGCggg                                                                                    <  2:386324/150‑1 (MQ=255)
                                                            cATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGa                                                                                   >  2:482507/1‑150 (MQ=255)
                                                                aatcaGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAATCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTc                                                                               <  2:428995/150‑1 (MQ=255)
                                                                aatcaGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTc                                                                               <  1:707530/150‑1 (MQ=255)
                                                                aatcaGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTc                                                                               <  1:1328428/150‑1 (MQ=255)
                                                                   caGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAAc                                                                            >  2:1242673/1‑150 (MQ=255)
                                                                       cGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCAc                                                                        >  1:727379/1‑150 (MQ=255)
                                                                          accaccAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAg                                                                     >  1:950851/1‑150 (MQ=255)
                                                                          accaccAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAg                                                                     >  1:623447/1‑150 (MQ=255)
                                                                           ccaccaGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAgg                                                                    >  2:68859/1‑150 (MQ=255)
                                                                            caccaGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAggg                                                                   <  1:293823/150‑1 (MQ=255)
                                                                             accaGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGc                                                                  >  2:1277134/1‑150 (MQ=255)
                                                                              ccaGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCt                                                                 >  2:1293971/1‑150 (MQ=255)
                                                                                  gCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCg                                                             >  1:1480732/1‑150 (MQ=255)
                                                                                       tACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGttttt                                                        >  2:69817/1‑150 (MQ=255)
                                                                                        aCCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTc                                                       >  1:496083/1‑150 (MQ=255)
                                                                                        aCCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTc                                                       >  2:1120687/1‑150 (MQ=255)
                                                                                        aCCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTc                                                       <  2:644742/150‑1 (MQ=255)
                                                                                         ccAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTc                                                                                                 <  2:1603365/107‑1 (MQ=255)
                                                                                         ccAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTc                                                                                                 >  1:1603365/1‑107 (MQ=255)
                                                                                             caaacaaTTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCg                                                  <  2:950851/150‑1 (MQ=255)
                                                                                              aaacaaTTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCgg                                                 <  1:125682/150‑1 (MQ=255)
                                                                                              aaacaaTTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCgg                                                 <  2:427763/150‑1 (MQ=255)
                                                                                                  aaTTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGc                                             <  2:1480732/150‑1 (MQ=255)
                                                                                                  aaTTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGc                                             <  2:448273/150‑1 (MQ=255)
                                                                                                      gCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCt                                         >  1:1281636/1‑150 (MQ=255)
                                                                                                      gCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCt                                         >  2:1426262/1‑150 (MQ=255)
                                                                                                       ccccTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTc                                        <  2:1525571/150‑1 (MQ=255)
                                                                                                         ccTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTccc                                      >  2:203721/1‑149 (MQ=255)
                                                                                                           tGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTccagc                                    <  1:1242673/150‑1 (MQ=255)
                                                                                                             ttCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATg                                              >  1:828599/1‑138 (MQ=255)
                                                                                                             ttCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATg                                              <  2:828599/138‑1 (MQ=255)
                                                                                                              tCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTccagccag                                 >  1:643439/1‑150 (MQ=255)
                                                                                                                aGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAACCGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCACCGCACCAGGGCTGGCGGTTTTCTGCCGGATGCCGCGCCAGCCACTc                               >  2:721402/1‑150 (MQ=255)
                                                                                                                 ggCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATAACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCg                              >  2:715349/1‑150 (MQ=255)
                                                                                                                        aGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAatat                       >  1:833074/1‑150 (MQ=255)
                                                                                                                           cGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAAACGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTgg                    >  2:719126/1‑150 (MQ=255)
                                                                                                                            gCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGc                   <  1:1512735/150‑1 (MQ=255)
                                                                                                                             gtttggcgaagcgCAACTGGCGCATCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCa                  >  1:576/6‑150 (MQ=255)
                                                                                                                             gtttggcgaagcgCAACTGGCGCATCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGGTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCa                  >  2:1207823/6‑150 (MQ=255)
                                                                                                                                tggcgaagcgCAACTGGCGCATCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCCGCCAGTCGGCAATATTGGCAAtt               >  2:161232/3‑150 (MQ=255)
                                                                                                                                atgcgatgcgCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAAtt               <  2:598709/150‑1 (MQ=255)
                                                                                                                                  gcgaagcgCAACTGGCGCATCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTaa             >  1:849743/1‑150 (MQ=255)
                                                                                                                                      tgcgCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGa         >  1:247071/1‑150 (MQ=255)
                                                                                                                                      tgcgCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGa         <  2:870539/150‑1 (MQ=255)
                                                                                                                                         gcCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGAATTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGt      >  2:714437/1‑150 (MQ=255)
                                                                                                                                           cAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTcc    >  2:1326620/1‑150 (MQ=255)
                                                                                                                                           cAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTcc    >  2:442306/1‑150 (MQ=255)
                                                                                                                                            aGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTCcg   >  1:258291/1‑150 (MQ=255)
                                                                                                                                            aGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTCcg   >  1:264110/1‑150 (MQ=255)
                                                                                                                                             gTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTCcgc  >  2:1134086/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CGATACCGCTGGCCCCACGGTTGCTGTACACCGGGTAACCTGCCGGAAGTTGCGAAAGCGCATCAATCAGACGTACCACCAGGCTGTTACCAACAAACAATTGCCCCTGTTCAGGCAGATAGTCGCAGATGCGATGCGCCAGTTGCGCTTCGCCAAACGCATCACGGCGGGCAATAACCGCCTGCATTGCCTGTTCCGCCAGGCGCGGGATTTCAACGCACCAGGGCTGGCGTTTTTCTGCCGGATGCAGCTCCAGCCAGTCGGCAATATTGGCAATTAAGCGACGTCCGC  >  NC_000913/2378005‑2378295

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: