Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,952 T→A 5.6% intergenic (+56/‑39) potF → / → potG putrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,9520TA 5.6% 262.4 / 5.6 72intergenic (+56/‑39)potF/potGputrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

TCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATG  >  NC_000913/894806‑895099
                                                                                                                                                  |                                                                                                                                                    
tCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTa                                                                                                                                                     >  4:507720/1‑148 (MQ=255)
  gTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACtattat                                                                                                                                                    <  4:928188/147‑1 (MQ=255)
  gTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                                                                                                    >  3:928188/1‑147 (MQ=255)
          cTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCtt                                                                                                                                          <  1:573112/148‑1 (MQ=255)
                aCTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacac                                                                                                                                    <  3:507720/148‑1 (MQ=255)
                 ctctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGc                                                                                                                                        >  4:275442/1‑143 (MQ=255)
                 ctctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGc                                                                                                                                        <  3:275442/143‑1 (MQ=255)
                    tGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc                                                                                                                                >  2:716206/1‑148 (MQ=255)
                           gCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                                                                                                    >  2:508035/1‑122 (MQ=255)
                           gCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                                                                                                    <  1:508035/122‑1 (MQ=255)
                             aGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCAtt                                                                                                                                                      >  4:69758/1‑118 (MQ=255)
                             aGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCAtt                                                                                                                                                      <  3:69758/118‑1 (MQ=255)
                              ggATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgc                                                                                                                                               >  2:941109/1‑124 (MQ=255)
                              ggATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCATAGATGCCAGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgc                                                                                                                                               <  1:941109/124‑1 (MQ=255)
                                aTCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                                                                                                    <  3:244449/117‑1 (MQ=255)
                                aTCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                                                                                                    >  4:244449/1‑117 (MQ=255)
                                 tCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagc                                                                                                                   <  1:716206/148‑1 (MQ=255)
                                  ccGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGcc                                                                                                                          >  2:94057/1‑140 (MQ=255)
                                  ccGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGcc                                                                                                                          <  1:94057/140‑1 (MQ=255)
                                  ccGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGgag                                                                                                                      >  4:676413/1‑144 (MQ=255)
                                  ccGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGgag                                                                                                                      <  3:676413/144‑1 (MQ=255)
                                     aaaaTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcg                                                                                                                                              <  4:747644/118‑1 (MQ=255)
                                     aaaaTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcg                                                                                                                                              >  3:747644/1‑118 (MQ=255)
                                     aaaaTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                               >  2:696042/1‑148 (MQ=255)
                                        aTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTg                                                                                                                                         <  3:957638/119‑1 (MQ=255)
                                        aTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTg                                                                                                                                         >  4:957638/1‑119 (MQ=255)
                                         tCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTg                                                                                                                                                >  2:121276/1‑112 (MQ=255)
                                         tCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTg                                                                                                                                                <  1:121276/112‑1 (MQ=255)
                                            aCCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATgg                                                                                                                                                  <  1:464306/107‑1 (MQ=255)
                                            aCCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATgg                                                                                                                                                  >  2:464306/1‑107 (MQ=255)
                                                    cgcACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATc                                                                                                <  1:696042/148‑1 (MQ=255)
                                                        cccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATg                                                                                                                            <  2:285397/116‑1 (MQ=255)
                                                        cccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATg                                                                                                                            >  1:285397/1‑116 (MQ=255)
                                                         ccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCg                                                                                           >  4:886013/1‑148 (MQ=255)
                                                              cgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGAAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTg                                                                                                            >  3:394341/1‑126 (MQ=255)
                                                              cgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGAAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTg                                                                                                            <  4:394341/126‑1 (MQ=255)
                                                                    ccAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                   >  1:1058003/1‑97 (MQ=255)
                                                                    ccAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                   <  2:1058003/97‑1 (MQ=255)
                                                                     cAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagcagc                                                                                                                <  1:429615/115‑1 (MQ=255)
                                                                     cAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagcagc                                                                                                                >  2:429615/1‑115 (MQ=255)
                                                                        aGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTg                                                                                                                                         <  1:899576/87‑1 (MQ=255)
                                                                        aGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTg                                                                                                                                         >  2:899576/1‑87 (MQ=255)
                                                                        aGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGa                                                                                                           >  4:850450/1‑117 (MQ=255)
                                                                        aGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGa                                                                                                           <  3:850450/117‑1 (MQ=255)
                                                                                                     tAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGc                                                                                      <  4:329227/109‑1 (MQ=255)
                                                                                                     tAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGc                                                                                      >  3:329227/1‑109 (MQ=255)
                                                                                                          gCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGc                                          >  2:350761/1‑148 (MQ=255)
                                                                                                             ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCg                                                                         <  3:1003217/114‑1 (MQ=255)
                                                                                                             ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCg                                                                         >  4:1003217/1‑114 (MQ=255)
                                                                                                             ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcg                                                            <  3:316654/127‑1 (MQ=255)
                                                                                                             ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcg                                                            >  4:316654/1‑127 (MQ=255)
                                                                                                             ggAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAAc                                       <  3:886013/148‑1 (MQ=255)
                                                                                                                    cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                          <  3:528639/74‑1 (MQ=255)
                                                                                                                    cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                          >  4:528639/1‑74 (MQ=255)
                                                                                                                    cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCa                                >  1:460778/1‑148 (MQ=255)
                                                                                                                       accacACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAAt                             >  3:32055/1‑148 (MQ=255)
                                                                                                                         cacaCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGc                                                                                                    <  3:1102240/75‑1 (MQ=255)
                                                                                                                         cacaCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGc                                                                                                    >  4:1102240/1‑75 (MQ=255)
                                                                                                                              ccgccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGa                      <  1:350761/148‑1 (MQ=255)
                                                                                                                               cgccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGaaa                                                                              >  1:842246/1‑91 (MQ=255)
                                                                                                                               cgccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGaaa                                                                              <  2:842246/91‑1 (MQ=255)
                                                                                                                                    gCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACg                                                                                                     <  1:1122858/63‑1 (MQ=255)
                                                                                                                                    gCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACg                                                                                                     >  2:1122858/1‑63 (MQ=255)
                                                                                                                                    gCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                               >  2:223291/1‑133 (MQ=255)
                                                                                                                                    gCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                               <  1:223291/133‑1 (MQ=255)
                                                                                                                                       aTTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcg                                                                  <  4:313790/95‑1 (MQ=255)
                                                                                                                                       aTTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcg                                                                  >  3:313790/1‑95 (MQ=255)
                                                                                                                                         gcGCACCATAATGGTGCGAATTGCACACATTCAATGCCGgag                                                                                                                      >  2:490500/2‑42 (MQ=9)
                                                                                                                                         gcGCACCATAATGGTGCGAATTGCACACATTCAATGCCGgag                                                                                                                      <  1:490500/41‑1 (MQ=9)
                                                                                                                                            cACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCgg        >  1:147907/1‑148 (MQ=255)
                                                                                                                                              ccATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcgc                                                                 <  4:668305/89‑1 (MQ=255)
                                                                                                                                              ccATTATGGTCCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcgc                                                                 >  3:668305/1‑89 (MQ=255)
                                                                                                                                               cATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTgg     <  4:32055/148‑1 (MQ=255)
                                                                                                                                                aTTATGGTGCGCTTGCACACATTCAATGCCGAAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgc                                                           >  1:1006037/1‑93 (MQ=255)
                                                                                                                                                aTTATGGTGCGCTTGCACACATTCAATGCCGAAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgc                                                           <  2:1006037/93‑1 (MQ=255)
                                                                                                                                                 ttATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCgg        <  4:846162/143‑1 (MQ=255)
                                                                                                                                                 ttATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCgg        >  3:846162/1‑143 (MQ=255)
                                                                                                                                                  tATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgatg  <  2:460778/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATG  >  NC_000913/894806‑895099

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