Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,515 G→C G61G (GGC→GGG yigI ← 4HBT thioesterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,5150GC96.7% 328.7 / 1.9 91G61G (GGC→GGGyigI4HBT thioesterase family protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (40/48);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCA  >  NC_000913/4004421‑4004609
                                                                                              |                                                                                              
gTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTg                                                                                           <  1:267074/100‑1 (MQ=255)
gTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTg                                                                                           <  2:723998/100‑1 (MQ=255)
 tAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGc                                                                                          <  2:1745773/100‑1 (MQ=255)
  aGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCa                                                                                         >  1:2043108/1‑100 (MQ=255)
   gTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCaa                                                                                        >  2:1857205/1‑100 (MQ=255)
    ttCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCaaa                                                                                       <  1:1148135/100‑1 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAAt                                                                                     >  2:247134/1‑100 (MQ=255)
        tCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGc                                                                                   >  1:593741/1‑100 (MQ=255)
          ttCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCtt                                                                                 >  2:44025/1‑100 (MQ=255)
            cACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCCGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTg                                                                               >  1:1113985/1‑100 (MQ=255)
            cACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTg                                                                               >  2:65011/1‑100 (MQ=255)
             aCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgc                                                                              <  1:1232113/100‑1 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcg                                                                             <  1:826333/100‑1 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcg                                                                             <  1:1809701/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                  <  2:423284/100‑1 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                <  2:989976/100‑1 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                <  2:303546/100‑1 (MQ=255)
                            gatctttAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccacc                                                               <  1:283982/95‑1 (MQ=255)
                              gggTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcat                                                             <  1:691930/100‑1 (MQ=255)
                               ggTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatc                                                            >  1:1787405/1‑100 (MQ=255)
                               ggTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatc                                                            >  1:530444/1‑100 (MQ=255)
                                 ttAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccaccagatc                                                          >  2:283981/1‑96 (MQ=255)
                                  tAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCAtt                                                         >  1:716582/1‑100 (MQ=255)
                                  tAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCAtt                                                         >  2:1848983/1‑100 (MQ=255)
                                    aGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTgg                                                       <  2:1066586/100‑1 (MQ=255)
                                    aGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTgg                                                       >  1:542483/1‑100 (MQ=255)
                                     ggTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGc                                                      >  2:507594/1‑100 (MQ=255)
                                      gTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCt                                                     <  2:221843/100‑1 (MQ=255)
                                        aCTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGa                                                   >  2:1940769/1‑100 (MQ=255)
                                         cTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                  >  1:1506937/1‑100 (MQ=255)
                                         cTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                  >  1:1932323/1‑100 (MQ=255)
                                          ttCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAtt                                                 <  2:1519526/100‑1 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                >  2:190537/1‑100 (MQ=255)
                                             ccACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttttt                                              <  1:767057/100‑1 (MQ=255)
                                             ccACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttttt                                              >  1:1073274/1‑100 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                             <  1:972095/100‑1 (MQ=255)
                                                cGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaaa                                           >  1:1215874/1‑100 (MQ=255)
                                                cGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaaa                                           <  2:1982474/100‑1 (MQ=255)
                                                cGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaaa                                           <  1:1765777/100‑1 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                          <  2:593739/100‑1 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                          >  1:65452/1‑100 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                          <  1:572369/100‑1 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                          <  2:644997/100‑1 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                         <  1:507595/100‑1 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                         >  1:1047352/1‑100 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                         <  2:736204/100‑1 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                         >  2:545830/1‑100 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                       >  2:681138/1‑100 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                       <  2:647697/100‑1 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                       <  1:1857222/100‑1 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                       >  2:132765/1‑100 (MQ=255)
                                                     acCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCa                                      >  1:1013153/1‑100 (MQ=255)
                                                      ccAGACCTGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                     >  2:1072779/1‑100 (MQ=255)
                                                      ccAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                     >  2:887330/1‑100 (MQ=255)
                                                       cAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                    >  1:1222244/1‑100 (MQ=255)
                                                       cAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                    <  1:748165/100‑1 (MQ=255)
                                                       cAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                    >  1:1027472/1‑100 (MQ=255)
                                                         gACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtg                                  >  1:1684908/1‑100 (MQ=255)
                                                         gACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtg                                  <  2:408016/100‑1 (MQ=255)
                                                         gACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtg                                  >  1:1949242/1‑100 (MQ=255)
                                                            cggcggCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGc                               <  2:236722/100‑1 (MQ=255)
                                                                ggcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAc                           <  1:65011/100‑1 (MQ=255)
                                                                ggcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAc                           >  1:1542143/1‑100 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                          <  2:1787390/100‑1 (MQ=255)
                                                                    aCATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtt                       <  1:1105146/100‑1 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                     <  1:247135/100‑1 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                     >  2:1936047/1‑100 (MQ=255)
                                                                       tCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttttt                    <  2:530825/100‑1 (MQ=255)
                                                                       tCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttttt                    <  2:1542133/100‑1 (MQ=255)
                                                                       tCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttttt                    >  1:984787/1‑100 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                   <  1:1848999/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                   >  2:1666803/1‑100 (MQ=255)
                                                                         cAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCg                  <  2:1717575/100‑1 (MQ=255)
                                                                         cAGCGCCGACGCAATGACCCCCCCGTGCAAAACGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCg                  <  1:479784/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                 <  2:2043090/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                 <  2:173353/100‑1 (MQ=255)
                                                                           gcgcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTa                <  1:920809/100‑1 (MQ=255)
                                                                            cgcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTan               <  2:530443/100‑2 (MQ=255)
                                                                             gcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAAc              <  1:772312/100‑1 (MQ=255)
                                                                              ccGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACg             >  2:1321901/1‑100 (MQ=255)
                                                                              ccGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACg             >  1:1490447/1‑100 (MQ=255)
                                                                               cGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGc            <  2:1684895/100‑1 (MQ=255)
                                                                                gACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCt           <  1:550439/100‑1 (MQ=255)
                                                                                  cGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc         <  1:100002/100‑1 (MQ=255)
                                                                                  cGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc         <  1:190538/100‑1 (MQ=255)
                                                                                   gCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCa        <  1:1940787/100‑1 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg       >  2:51912/1‑100 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg       >  2:834793/1‑100 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg       >  1:366131/1‑100 (MQ=255)
                                                                                      aTGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGtt     <  1:545831/100‑1 (MQ=255)
                                                                                         aCCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa  >  2:841303/1‑100 (MQ=255)
                                                                                              |                                                                                              
GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCA  >  NC_000913/4004421‑4004609

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