Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,501,908 C→T 9.1% intergenic (+319/+195) yjgZ → / ← insG KpLE2 phage‑like element; uncharacterized protein YjgZ/KpLE2 phage‑like element; IS4 putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,501,9080CT 9.1% 289.0 / 14.0 100intergenic (+319/+195)yjgZ/insGKpLE2 phage‑like element; uncharacterized protein YjgZ/KpLE2 phage‑like element; IS4 putative transposase
Reads supporting (aligned to +/- strand):  ref base C (50/39);  new base T (5/4);  total (56/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAG  >  NC_000913/4501760‑4502055
                                                                                                                                                    |                                                                                                                                                   
gTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCa                                                                                                                                                    <  2:1226922/150‑1 (MQ=255)
 tACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCaa                                                                                                                                                   <  1:890325/150‑1 (MQ=255)
  aCCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAAc                                                                                                                                                  >  1:1184562/1‑150 (MQ=255)
     aGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTgg                                                                                                                                               >  2:457481/1‑150 (MQ=255)
        ctagctagGCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAATTGCGTTAATGTTGAGTTCCTCGAAGCCTTAGCAATTAAGGGACATAGCCAGTTGCGTATCTgggta                                                                                                                                            <  2:778213/150‑4 (MQ=255)
         tagctagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGa                                                                                                                                           <  2:462292/150‑1 (MQ=255)
           gctagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGGAACTGGCTGATg                                                                                                                                         <  1:920599/150‑1 (MQ=255)
           gctagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATg                                                                                                                                         >  2:890951/1‑150 (MQ=255)
           gctagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATg                                                                                                                                         >  1:1504842/1‑150 (MQ=255)
             tagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGTTGCGTATCTGGCtatgtcc                                                                                                                                       <  1:1332631/150‑7 (MQ=255)
             tagcCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGTTGCGTATCTGGCtatgtcc                                                                                                                                       <  1:785632/150‑7 (MQ=255)
              acccccTACTCAACCAAATAAATCAACAATCGTCAAAATGAAATAGTTCAAATACACTTTCAATAAAACGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCaa                                                                                                                                      >  2:245422/3‑150 (MQ=255)
                        cAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACACTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAAt                                                                                                                            >  2:260094/1‑150 (MQ=255)
                          aCCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCATACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGt                                                                                                                          >  2:409071/1‑150 (MQ=255)
                           ccAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTc                                                                                                                         >  2:157333/1‑150 (MQ=255)
                            cAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCt                                                                                                                        >  1:1314336/1‑150 (MQ=255)
                              aataaatCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCttt                                                                                                                      <  1:798224/150‑1 (MQ=255)
                              aataaatCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCttt                                                                                                                      <  2:407884/150‑1 (MQ=255)
                               ataaatCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTa                                                                                                                     <  1:860194/150‑1 (MQ=255)
                               ataaatCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTa                                                                                                                     <  2:354388/150‑1 (MQ=255)
                                   atCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAg                                                                                                                    <  2:1255252/147‑1 (MQ=255)
                                   atCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAg                                                                                                                    >  1:1255252/1‑147 (MQ=255)
                                   atCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATc                                                                                                                 >  1:722155/1‑150 (MQ=255)
                                    tCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATcc                                                                                                                <  2:369451/150‑1 (MQ=255)
                                      aaCCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGc                                                                                                              >  2:1525208/1‑150 (MQ=255)
                                        ccATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTg                                                                                                            <  1:812419/150‑1 (MQ=255)
                                             gTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCt                                                                                                                               <  2:765860/126‑1 (MQ=255)
                                             gTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCt                                                                                                                               >  1:765860/1‑126 (MQ=255)
                                             gTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTg                                                                                                       <  1:457481/150‑1 (MQ=255)
                                              tCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGc                                                                                                      >  2:1303995/1‑150 (MQ=255)
                                               cAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCg                                                                                                     >  2:415721/1‑150 (MQ=255)
                                                 gagTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCggg                                                                                                   >  1:746072/1‑150 (MQ=255)
                                                 gagTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCggg                                                                                                   >  1:751474/1‑150 (MQ=255)
                                                     gaaaTCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATg                                                                                               <  2:1446394/148‑1 (MQ=255)
                                                      cAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGa                                                                                              >  1:382234/1‑150 (MQ=255)
                                                      cAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGa                                                                                              >  1:639292/1‑150 (MQ=255)
                                                       aaTCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGAc                                                                                             <  1:1017261/150‑1 (MQ=255)
                                                       aaTCGTTCATATACACTGTCGCTGATGCGCTAAAAAATGCGTTAATGTTCCGTTCCTCGAACCGTGACCAATTAAGGTCCATAGCAAGATCCGAACCTGGCTGATGAAGGAAGTATAATGTGTTTAGCTCCGGGGGCCAGCGGGCATGAc                                                                                             >  2:697006/1‑150 (MQ=255)
                                                              cATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCa                                                                                      >  1:1144073/1‑150 (MQ=255)
                                                                tataCACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCaaa                                                                                    >  2:270660/1‑150 (MQ=255)
                                                                 ataCACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATACGCTGGCCTGCGGGCATGACAAAGTCaaaa                                                                                   >  2:1112043/1‑149 (MQ=255)
                                                                  taCACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGCCAAACa                                                                                  >  1:416878/1‑150 (MQ=255)
                                                                    cacTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGt                                                                                >  1:27158/1‑150 (MQ=255)
                                                                      cTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCa                                                                              >  2:844355/1‑150 (MQ=255)
                                                                        gTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCaaa                                                                            <  2:561673/150‑1 (MQ=255)
                                                                         tCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCATACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAg                                                                           <  1:409071/150‑1 (MQ=255)
                                                                           ggTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGAGAATGTCAAACAGTCAAAGAt                                                                         <  2:1210098/150‑1 (MQ=255)
                                                                                   gcTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGAAGTCTAATGTCTTTAGATCCTCTGGCATGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGt                                                                 >  1:414584/1‑150 (MQ=255)
                                                                                     tAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTg                                                               <  1:245422/150‑1 (MQ=255)
                                                                                      aaaaaaTGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGa                                                              >  1:425329/1‑150 (MQ=255)
                                                                                      aaaaaaTGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGa                                                              >  2:1534085/1‑150 (MQ=255)
                                                                                           aTGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTc                                                         <  1:191170/150‑1 (MQ=255)
                                                                                              cGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGCCAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGt                                                      >  1:1150972/1‑150 (MQ=255)
                                                                                               gTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTc                                                     >  1:1041122/1‑150 (MQ=255)
                                                                                                  aaTGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTg                                                  <  2:947956/150‑1 (MQ=255)
                                                                                                  aaTGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTg                                                  <  1:879702/150‑1 (MQ=255)
                                                                                                   aTGTTGAGTTCCTCGAAGCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCGGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTgg                                                 <  1:1526927/150‑1 (MQ=255)
                                                                                                    tGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTggg                                                <  2:877588/150‑1 (MQ=255)
                                                                                                     gTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTgggg                                               <  1:1534085/150‑1 (MQ=255)
                                                                                                     gTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTgggg                                               <  1:1303995/150‑1 (MQ=255)
                                                                                                       tGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGAc                                             <  1:844355/150‑1 (MQ=255)
                                                                                                       tGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGAc                                             <  1:1139284/150‑1 (MQ=255)
                                                                                                        gAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACa                                            >  1:517876/1‑150 (MQ=255)
                                                                                                          gTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACaaa                                          >  1:1146758/1‑150 (MQ=255)
                                                                                                             ccTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACaaaaaa                                       >  2:279171/1‑150 (MQ=255)
                                                                                                               tCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGTCGTTGACGGTCAGTCTTGGGGACAAAAAAGc                                     <  2:1314336/150‑1 (MQ=255)
                                                                                                               tCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGc                                     <  1:1058468/150‑1 (MQ=255)
                                                                                                               tCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGc                                     >  1:352428/1‑150 (MQ=255)
                                                                                                               tCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGc                                     <  2:722155/150‑1 (MQ=255)
                                                                                                                cGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAACGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCt                                    >  1:509581/1‑150 (MQ=255)
                                                                                                                     ccTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTa                               <  1:1112043/150‑1 (MQ=255)
                                                                                                                      cTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTAt                              >  1:645018/1‑150 (MQ=255)
                                                                                                                         aCCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTc                                   <  1:648386/142‑1 (MQ=255)
                                                                                                                         aCCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTc                                   >  2:648386/1‑142 (MQ=255)
                                                                                                                          ccAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATtctc                          >  2:792057/1‑150 (MQ=255)
                                                                                                                            aaTTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGt                        >  1:1040940/1‑150 (MQ=255)
                                                                                                                             aTTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGtt                       >  2:42815/1‑150 (MQ=255)
                                                                                                                             aTTAAGGCACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTTGATCCGATGCACTGAGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGCCGGTCAGTCTTGGGGACAAAAAACCTAAGTATTATCGtt                       >  2:1526147/1‑150 (MQ=255)
                                                                                                                                   ggACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAAc                 >  1:1531333/1‑150 (MQ=255)
                                                                                                                                    gcatcaGCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAAcc                >  2:392234/6‑150 (MQ=255)
                                                                                                                                    gACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGAGTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAAcc                >  2:975606/1‑150 (MQ=255)
                                                                                                                                     catcaGCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCa               <  1:392234/146‑1 (MQ=255)
                                                                                                                                     catcaGCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCa               >  1:942409/5‑150 (MQ=255)
                                                                                                                                     catcaGCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCa               >  1:936806/5‑150 (MQ=255)
                                                                                                                                     aCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCa               <  2:817508/150‑1 (MQ=255)
                                                                                                                                        tAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGAt            >  2:191241/1‑150 (MQ=255)
                                                                                                                                        caGCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGAt            >  2:800133/2‑150 (MQ=255)
                                                                                                                                         aGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGATGGCAACCAGATa           >  1:661084/1‑150 (MQ=255)
                                                                                                                                          gCCAGTTGCGTATCTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATaa          >  2:263113/1‑150 (MQ=255)
                                                                                                                                          gCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATaa          >  2:787554/1‑150 (MQ=255)
                                                                                                                                             aGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAAttt       <  2:517876/150‑1 (MQ=255)
                                                                                                                                             aGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAAttt       <  1:270660/150‑1 (MQ=255)
                                                                                                                                             aGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAAttt       >  1:1210447/1‑150 (MQ=255)
                                                                                                                                              gATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTa      >  1:397480/1‑150 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGTACTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTATCGTTGGCAACCAGAAAATTTAACAg  >  2:613388/1‑150 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTg                      >  1:36173/1‑130 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTg                      <  2:36173/130‑1 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAg  <  2:1504842/150‑1 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAg  <  2:525984/150‑1 (MQ=255)
                                                                                                                                                  cGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTAAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAg  <  1:415721/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAG  >  NC_000913/4501760‑4502055

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