Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,571 T→C 6.8% intergenic (+4/+155) ykgQ → / ← rclC pseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5710TC 6.8% 234.5 / 6.3 73intergenic (+4/+155)ykgQ/rclCpseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/39);  new base C (2/3);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTG  >  NC_000913/317427‑317716
                                                                                                                                                |                                                                                                                                                 
ccATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                                                                               >  1:301240/1‑149 (MQ=255)
 cATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                              >  3:468655/1‑149 (MQ=255)
  aTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCAt                                                                                                                                             <  4:69175/149‑1 (MQ=255)
  aTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCAt                                                                                                                                             <  1:601016/149‑1 (MQ=255)
   tCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                                                                            >  4:129975/1‑149 (MQ=255)
     aaaaaTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                                                                          <  2:488428/149‑1 (MQ=255)
         aTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                                                                      >  3:417871/1‑149 (MQ=255)
          tCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                                                                     <  1:184576/149‑1 (MQ=255)
            gATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                                                                   <  2:301240/149‑1 (MQ=255)
              tttCTCACTGGCAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                 >  1:651935/1‑149 (MQ=255)
                 ctcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAga                                                                                                                              <  4:75766/149‑1 (MQ=255)
                  tcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                                                                          <  2:261671/136‑1 (MQ=255)
                  tcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                                                                          >  1:261671/1‑136 (MQ=255)
                  tcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCAt                                                                                                                                             >  2:227993/1‑133 (MQ=255)
                  tcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCAt                                                                                                                                             <  1:227993/133‑1 (MQ=255)
                    aCTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                                                                           <  4:468655/149‑1 (MQ=255)
                    aCTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                                                                           <  4:629978/149‑1 (MQ=255)
                    aCTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                                                                           <  2:321044/149‑1 (MQ=255)
                        aCAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaa                                                                                                                       <  1:154128/149‑1 (MQ=255)
                          aGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaaaa                                                                                                                     <  3:402251/149‑1 (MQ=255)
                           gATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                                                                    <  4:555163/149‑1 (MQ=255)
                                     tGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                                                                          <  2:579574/149‑1 (MQ=255)
                                        cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                                                                       <  3:129975/149‑1 (MQ=255)
                                                aGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                                                               <  4:263729/149‑1 (MQ=255)
                                                   tCGCTTGACTTGGGCTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGTTATGAAAATACATGACTGATAGTAATCATGTATTATACTTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGAt                                                                                            <  3:427522/149‑1 (MQ=255)
                                                     gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                                                                          >  2:624580/1‑149 (MQ=255)
                                                           tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                                                                    >  3:482000/1‑149 (MQ=255)
                                                            ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                                                                   <  1:78346/149‑1 (MQ=255)
                                                               ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                >  4:87495/1‑149 (MQ=255)
                                                                gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                                                                               <  3:87495/149‑1 (MQ=255)
                                                                 ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAt                                                                                                                                                <  1:10820/83‑1 (MQ=255)
                                                                 ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAt                                                                                                                                                >  2:10820/1‑83 (MQ=255)
                                                                     gTTGGGGAGTCGGTAAAACATTTTAATGCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                                                                   >  3:461539/1‑140 (MQ=255)
                                                                     gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                                                                   <  4:461539/140‑1 (MQ=255)
                                                                      ttGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                                                                         >  4:340113/1‑149 (MQ=255)
                                                                      ttGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                                                                         >  3:153207/1‑149 (MQ=255)
                                                                             gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat                                                                  <  4:153207/149‑1 (MQ=255)
                                                                                 gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                                                                   <  1:51943/128‑1 (MQ=255)
                                                                                 gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                                                                   >  2:51943/1‑128 (MQ=255)
                                                                                 gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                              >  1:365828/1‑149 (MQ=255)
                                                                                  tAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGc                                                             <  3:340113/149‑1 (MQ=255)
                                                                                     aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc                                                          >  3:461673/1‑149 (MQ=255)
                                                                                      aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                                                               <  1:158736/111‑1 (MQ=255)
                                                                                      aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                                                               >  2:158736/1‑111 (MQ=255)
                                                                                         ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc                                                      >  2:472577/1‑149 (MQ=255)
                                                                                              aTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcat                                                 <  1:624580/149‑1 (MQ=255)
                                                                                                 cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcat                                              >  2:469156/1‑149 (MQ=255)
                                                                                                  tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACCGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcatc                                             <  4:417871/149‑1 (MQ=255)
                                                                                                         gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                      <  4:308393/149‑1 (MQ=255)
                                                                                                          ttCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAggg                                     <  2:365828/149‑1 (MQ=255)
                                                                                                             gTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGc                                  <  4:482000/149‑1 (MQ=255)
                                                                                                                gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAg                                                        <  1:610385/124‑1 (MQ=255)
                                                                                                                gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAg                                                        >  2:610385/1‑124 (MQ=255)
                                                                                                                     cTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                                                                   <  2:337735/60‑1 (MQ=255)
                                                                                                                     cTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                                                                   >  1:337735/1‑60 (MQ=255)
                                                                                                                     cTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGcc                          >  2:282629/1‑149 (MQ=255)
                                                                                                                      ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTg                                                                                              <  1:610626/80‑1 (MQ=255)
                                                                                                                      ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTg                                                                                              >  2:610626/1‑80 (MQ=255)
                                                                                                                          tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGaaa                     >  2:544259/1‑149 (MQ=255)
                                                                                                                            gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCGTCACTCTCTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGaaaaa                   >  2:494340/3‑149 (MQ=255)
                                                                                                                              aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGcc                                 <  2:588060/133‑1 (MQ=255)
                                                                                                                              aaaaGACATGATTGATATTAATCACGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGcc                                 >  1:588060/1‑133 (MQ=255)
                                                                                                                                 aTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaa              <  1:327146/149‑1 (MQ=255)
                                                                                                                                 aTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaa              <  1:469156/149‑1 (MQ=255)
                                                                                                                                   aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                                                                    <  1:321037/93‑1 (MQ=255)
                                                                                                                                   aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTCTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                                                                    >  2:321037/1‑93 (MQ=255)
                                                                                                                                      tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAacaaca         >  1:344409/1‑149 (MQ=255)
                                                                                                                                        attgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAAGAacaga       <  4:461673/149‑1 (MQ=255)
                                                                                                                                         ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac      >  1:477221/1‑149 (MQ=255)
                                                                                                                                          tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagaca     <  2:344409/149‑1 (MQ=255)
                                                                                                                                          tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagaca     <  2:651935/149‑1 (MQ=255)
                                                                                                                                          taatatCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGGTCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagaca     <  1:494340/147‑1 (MQ=255)
                                                                                                                                             tatTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAgtg  >  1:404533/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTG  >  NC_000913/317427‑317716

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