Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,201 T→C 5.9% pseudogene (810/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,2010TC 5.9% 169.5 / 3.0 68pseudogene (810/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (26/37);  new base C (2/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

ATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATG  >  NC_000913/1466059‑1466346
                                                                                                                                              |                                                                                                                                                 
aTGCCACAGTAAACAACAATGGCAAAACCACCGATGACGGTAAAGATTCCACCGTTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgc                                                                                                                                             <  2:121562/149‑1 (MQ=255)
  gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgcgc                                                                                                                                           <  1:39433/149‑1 (MQ=255)
       aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACgg                                                                                                                                      >  1:737545/1‑149 (MQ=255)
        gTAAACAACAATGGCAATACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                                                                                     >  1:915083/1‑149 (MQ=255)
        gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                                                                                     <  1:733795/149‑1 (MQ=255)
          aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                                                                                   >  2:991184/1‑149 (MQ=255)
           aacaacAATCGCAACACCACCGGTGACGGTAAAGATTCCACCGGAACGGAAATCAACGCGAATCACGGGACAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTACTGATATCACCGGTGACAGCGCGACGGTGGa                                                                                                                                  <  2:782314/149‑1 (MQ=14)
            acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGTTGACAGCGCGACGGTGGAt                                                                                                                                 <  2:367973/149‑1 (MQ=255)
              aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAATTCCACGGTAATAACGGGAAAGTGATTCAGGACGGCGCTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATaa                                                                                                                               >  2:365076/1‑149 (MQ=255)
                   tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAgg                                                                                                                          <  1:1175988/149‑1 (MQ=255)
                     gCAAAACCACCGTTGACGGTAAAGATGCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                                                                                        <  1:843731/149‑1 (MQ=255)
                       aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAc                                                                                                                      >  1:883002/1‑149 (MQ=255)
                        aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                                                                                     >  1:95940/1‑149 (MQ=255)
                        aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                                                                                                             >  2:232227/1‑125 (MQ=38)
                        aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                                                                                                             <  1:232227/125‑1 (MQ=38)
                        aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCAcc                                                                                                                     >  1:416120/1‑149 (MQ=255)
                           ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATg                                                                                                                  >  2:726714/1‑149 (MQ=255)
                            caccGGTGACGTTATAGTTTACATCTGTACGGAAATCTACGTTAATAACGGGAAAGTGCTTCAGGTCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGa                                                                                                                 <  2:697785/149‑1 (MQ=14)
                              ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGAGATTCAGGACGGCGATATGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCTCGACCGTGGATAACAAGGGCACCATGccc                                                                                                               >  2:1096319/1‑149 (MQ=14)
                               cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                              <  1:144753/149‑1 (MQ=255)
                                   gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                          >  1:1152427/1‑149 (MQ=255)
                                    aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                                         <  2:416120/149‑1 (MQ=21)
                                          aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                   <  2:276557/149‑1 (MQ=40)
                                                              aTCAACGGTAATAACGCGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGGCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATGGGTATCCAGAt                                                                               <  2:170264/149‑1 (MQ=255)
                                                                 aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGa                                                                            >  1:734493/1‑149 (MQ=255)
                                                                 aaCGGTAAAAACGGGAAAGTGATTCAGTACGGCGATCTGGATGTCAGCGTCGGCGGTCACGGTATTGATATCACCGGTGACACCGCGACAGTGGATAACAAGGGCACCATGACCGTCACCGATCCGCAGTCCATGGATATCAAAATCGa                                                                            >  2:643744/1‑149 (MQ=14)
                                                                  aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                           >  1:1016531/1‑149 (MQ=255)
                                                                   cGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCAGCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACg                                                                          <  2:220484/149‑1 (MQ=255)
                                                                    ggTAATAATTGGAAAGTGATTCAGGACGGCTATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACAGATCCGGAGTCCATGGGTATCCAGATCGACgg                                                                         <  1:756671/149‑1 (MQ=255)
                                                                    ggTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCAAGATCGACgg                                                                         >  1:1107908/1‑149 (MQ=255)
                                                                     gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt                                                                        <  1:991184/149‑1 (MQ=255)
                                                                         taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCCCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa                                                                    >  2:822860/1‑149 (MQ=255)
                                                                         taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa                                                                    >  2:309284/1‑149 (MQ=255)
                                                                                      aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCaa                                                       >  2:609446/1‑149 (MQ=255)
                                                                                             ccGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATGACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGaa                                                <  1:221617/148‑1 (MQ=255)
                                                                                               ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAgg                                              >  1:1061057/1‑149 (MQ=255)
                                                                                               ggCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAgg                                              <  1:1089719/149‑1 (MQ=255)
                                                                                                  gATCTGTATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCga                                           <  1:1104449/149‑1 (MQ=255)
                                                                                                  gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCga                                           >  1:910301/1‑149 (MQ=255)
                                                                                                  gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCga                                           <  2:734493/149‑1 (MQ=255)
                                                                                                   aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgag                                          <  2:425686/149‑1 (MQ=255)
                                                                                                     cTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAtcg                                                                                         >  1:1094063/1‑98 (MQ=21)
                                                                                                        gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGCCCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                     <  2:95940/149‑1 (MQ=255)
                                                                                                          tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCa                                   <  1:365076/149‑1 (MQ=255)
                                                                                                          tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCa                                   >  1:845717/1‑149 (MQ=255)
                                                                                                              aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                               >  2:923971/1‑149 (MQ=255)
                                                                                                              aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                               >  2:265701/1‑149 (MQ=255)
                                                                                                              aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                               <  2:682746/149‑1 (MQ=255)
                                                                                                               gcggcggcgGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAcc                              <  2:883002/149‑1 (MQ=255)
                                                                                                                     gcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGt                        <  1:728312/149‑1 (MQ=255)
                                                                                                                     gcggTCACGGTATTGATATCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGt                        <  2:915083/149‑1 (MQ=255)
                                                                                                                       ggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTgg                      <  1:1096319/149‑1 (MQ=255)
                                                                                                                           aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcacc                  <  2:737545/149‑1 (MQ=255)
                                                                                                                                aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                              <  1:225151/68‑1 (MQ=40)
                                                                                                                                aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                              >  2:225151/1‑68 (MQ=40)
                                                                                                                                   gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGCTCCGGAGTCCATGGGTATCCAGAACGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCccccc          >  2:933365/1‑148 (MQ=255)
                                                                                                                                   gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa          <  1:726714/149‑1 (MQ=255)
                                                                                                                                     tatCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa        <  2:1152427/149‑1 (MQ=255)
                                                                                                                                      atCACCGGTTACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGTCCATGGGTATCAAGATCGACGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt       <  2:425269/149‑1 (MQ=255)
                                                                                                                                      atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt       <  2:429882/149‑1 (MQ=255)
                                                                                                                                       tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt      <  1:625379/149‑1 (MQ=255)
                                                                                                                                       tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGCCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGAtt      <  1:864326/149‑1 (MQ=255)
                                                                                                                                        cACCGGTGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                          <  2:313623/48‑1 (MQ=255)
                                                                                                                                        cACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                          >  1:313623/1‑48 (MQ=255)
                                                                                                                                         aCCGGTGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                      <  1:865043/51‑1 (MQ=255)
                                                                                                                                         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                      >  2:865043/1‑51 (MQ=255)
                                                                                                                                          ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt   <  2:210156/149‑1 (MQ=255)
                                                                                                                                           cGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg  <  1:976325/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATG  >  NC_000913/1466059‑1466346

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