Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,291,262 T→C F817L (TTT→CTT)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,2620TC100.0% 255.1 / NA 64F817L (TTT→CTT) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/42);  total (22/42)

CCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCC  >  NC_000913/3291118‑3291403
                                                                                                                                                |                                                                                                                                             
cccTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCt                                                                                                                                          <  1:582291/150‑1 (MQ=255)
  cTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTc                                                                                                                                       >  2:358276/1‑151 (MQ=255)
     tGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAg                                                                                                                                     <  1:677298/150‑1 (MQ=255)
      ggCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGt                                                                                                                                    <  1:668405/150‑1 (MQ=255)
              tGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGAt                                                                                                                           >  1:245197/1‑151 (MQ=255)
               gCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGAtt                                                                                                                          >  1:187021/1‑151 (MQ=255)
                   gCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGa                                                                                                                      <  1:613354/151‑1 (MQ=255)
                     aaCGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAAc                                                                                                                    <  2:406302/151‑1 (MQ=255)
                      aCGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAAcc                                                                                                                   <  1:251457/151‑1 (MQ=255)
                                      acaaAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTAtt                                                                                                    >  1:5194/1‑150 (MQ=255)
                                           gTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAgg                                                                                              <  1:613553/151‑1 (MQ=255)
                                                gtcgATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAg                                                                                         <  2:450595/151‑1 (MQ=255)
                                                     tGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAgcgc                                                                                      >  2:204318/1‑149 (MQ=255)
                                                      gAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAg                                                                                    <  2:226617/150‑1 (MQ=255)
                                                           cggcATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCAt                                                                              <  1:576021/151‑1 (MQ=255)
                                                                 aCTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAgg                                                                         <  1:8144/150‑1 (MQ=255)
                                                                 aCTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGa                                                                        <  2:110801/151‑1 (MQ=255)
                                                                 aCTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGa                                                                        <  1:367868/151‑1 (MQ=255)
                                                                    atatCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGAt                                                                       <  2:66502/149‑1 (MQ=255)
                                                                         tctctGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATaaaaaa                                                                 >  2:259418/1‑150 (MQ=255)
                                                                          ctctGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATaaaaaa                                                                 <  2:566386/149‑1 (MQ=255)
                                                                          ctctGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATaaaaaa                                                                 <  1:516064/149‑1 (MQ=255)
                                                                           tctGGGCTATCGAGAAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGCCACCTTGCTTTCAATACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAAt                                                               >  2:723815/1‑150 (MQ=255)
                                                                              gggATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGaa                                                            <  1:198887/150‑1 (MQ=255)
                                                                              gggATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGaa                                                            <  1:174333/150‑1 (MQ=255)
                                                                                gATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAag                                                         <  2:140120/151‑1 (MQ=255)
                                                                                gATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAag                                                         <  2:483097/151‑1 (MQ=255)
                                                                                 atatCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGaaa                                                           >  1:198711/1‑148 (MQ=255)
                                                                                 atatCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAaga                                                        <  1:535206/151‑1 (MQ=255)
                                                                                 atatCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAaga                                                        <  1:520073/151‑1 (MQ=255)
                                                                                    tCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAgcg                                                     >  1:237151/1‑151 (MQ=255)
                                                                                       aGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCct                                                  <  2:525231/151‑1 (MQ=255)
                                                                                        gTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCct                                                  <  2:110917/150‑1 (MQ=255)
                                                                                          aaaTCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCtt                                               <  2:440395/151‑1 (MQ=255)
                                                                                           aaTCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTta                                              <  1:407502/151‑1 (MQ=255)
                                                                                            aTCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATaaa                                                                    >  1:110850/1‑128 (MQ=255)
                                                                                            aTCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATaaa                                                                    <  2:110848/128‑1 (MQ=255)
                                                                                            aTCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTtat                                             >  2:594826/1‑151 (MQ=255)
                                                                                              cacaAAGCTGGACTGTACGCTGGGGAAATCAGGCTGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTtata                                            >  2:458811/1‑150 (MQ=255)
                                                                                                caAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACa                                         <  1:169656/151‑1 (MQ=255)
                                                                                                 aaaGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAg                                        <  2:5194/151‑1 (MQ=255)
                                                                                                    gCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTc                                     <  1:248111/151‑1 (MQ=255)
                                                                                                       ggACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCttt                                  <  2:112883/151‑1 (MQ=255)
                                                                                                       ggACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCtt                                   <  2:348221/150‑1 (MQ=255)
                                                                                                        gACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCtttt                                 >  1:232607/1‑151 (MQ=255)
                                                                                                            gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttgttg                             >  1:263192/1‑151 (MQ=255)
                                                                                                               cGCTGGGGAAATCAGTCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGat                           >  2:42540/1‑150 (MQ=255)
                                                                                                                     ggAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCcaca                    <  2:510785/151‑1 (MQ=255)
                                                                                                                     ggAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCcaca                    <  1:128126/151‑1 (MQ=255)
                                                                                                                     ggAAATCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCcaca                    <  1:49065/151‑1 (MQ=255)
                                                                                                                       aaaTCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATc                  <  1:256344/151‑1 (MQ=255)
                                                                                                                        aaTCAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATcc                 >  2:734196/1‑151 (MQ=255)
                                                                                                                           cAGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTgc              <  1:579190/151‑1 (MQ=255)
                                                                                                                            aGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTgcg             >  1:607551/1‑151 (MQ=255)
                                                                                                                            aGGCAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTgcg             >  2:624037/1‑151 (MQ=255)
                                                                                                                               cAGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCtt          <  2:474722/151‑1 (MQ=255)
                                                                                                                                aGATCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTa         <  1:173607/151‑1 (MQ=255)
                                                                                                                                  aTCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATg       >  2:443478/1‑151 (MQ=255)
                                                                                                                                  aTCAACAATGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATg       <  1:662652/151‑1 (MQ=255)
                                                                                                                                    caacaaTGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGc      <  2:76511/150‑1 (MQ=255)
                                                                                                                                     aacaaTGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGcct    >  1:350972/1‑151 (MQ=255)
                                                                                                                                      acaaTGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGcctc   <  2:528239/151‑1 (MQ=255)
                                                                                                                                      acaaTGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGcctc   >  1:226128/1‑151 (MQ=255)
                                                                                                                                       caaTGTCAGCTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGcctcc  >  1:306720/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCC  >  NC_000913/3291118‑3291403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: