Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,072,013 C→T 5.4% H159Y (CAT→TAT)  flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,072,0130CT 5.4% 239.5 / 2.4 74H159Y (CAT→TAT) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base T (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

AGCCAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTGTTC  >  NC_000913/2071865‑2072160
                                                                                                                                                    |                                                                                                                                                   
aGCCAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGc                                                                                                                                                     <  1:292687/149‑1 (MQ=255)
 gCCAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCa                                                                                                                                                    <  1:104925/149‑1 (MQ=255)
 gCCAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGATATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCa                                                                                                                                                    >  5:372987/1‑149 (MQ=255)
  ccAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCAt                                                                                                                                                   <  2:41842/149‑1 (MQ=255)
    aacaaAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCAt                                                                                                                                                   >  6:394648/1‑147 (MQ=255)
    aacaaAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGa                                                                                                                                                 <  7:80577/149‑1 (MQ=255)
     acaaAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAg                                                                                                                                                >  3:204578/1‑149 (MQ=255)
     acaaAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAg                                                                                                                                                >  8:266184/1‑149 (MQ=255)
         aaCGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAggggg                                                                                                                                            >  2:209598/1‑149 (MQ=255)
          aCGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGc                                                                                                                                           <  7:181346/149‑1 (MQ=255)
           cGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCg                                                                                                                                          <  7:147363/149‑1 (MQ=255)
            gACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGa                                                                                                                                         >  6:350315/1‑149 (MQ=255)
                   aCCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCcac                                                                                                                                  <  7:266184/149‑1 (MQ=255)
                     cAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAg                                                                                                                                >  8:184126/1‑149 (MQ=255)
                       gtggtgGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAgg                                                                                                                               >  7:265308/1‑148 (MQ=255)
                       gtggtgGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGCTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGa                                                                                                                              <  1:86402/149‑1 (MQ=255)
                          gtggtCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGTGATAGCCACAGGAAcc                                                                                                                           >  8:174298/1‑149 (MQ=255)
                                ttCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCAt                                                                                                                      >  1:213123/1‑148 (MQ=255)
                                  cAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTa                                                                                                                    <  6:193699/148‑1 (MQ=255)
                                   agagagTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGc                                                                                                                                                     >  4:369339/1‑114 (MQ=255)
                                   agagagTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGc                                                                                                                                                     <  3:369339/114‑1 (MQ=255)
                                     agagTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGa                                                                                                                <  1:209598/149‑1 (MQ=255)
                                      gagTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGAt                                                                                                               <  2:249155/149‑1 (MQ=255)
                                         tGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGGCTTCAGGGACGGGCTGTGAACACCAGGCTGATTGTTTGCGAACAGTGGATGCATGAGGGTGCGATAGCCACAGGAACCGTCATTAATGATAAg                                                                                                            <  4:181307/149‑1 (MQ=255)
                                         tGAACCCCGGGGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAg                                                                                                            <  3:332581/149‑1 (MQ=255)
                                           aaCCCCGGTGGAAGTGTCTCAGGCACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCGGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAgg                                                                                                           <  5:85713/148‑1 (MQ=255)
                                            aCCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATg                                                                                                                 <  4:332656/141‑1 (MQ=255)
                                            aCCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATg                                                                                                                 >  3:332656/1‑141 (MQ=255)
                                                cGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGc                                                                                                     <  3:345115/149‑1 (MQ=255)
                                                  gTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAg                                                                                                   <  4:204578/149‑1 (MQ=255)
                                                   tGGAAGTGTCTCAGACACGGTTATGAGTGCCGGAGTCGGACAGTGCCTTCTGGGTCGGGCTGTGAACACCACGCTGATTGGTGGCGAACAGTGGATGCATGAGGGGGCTATAGCCACTGGAACCGTCATTAATGATAATGGCTGGCAgg                                                                                                  <  3:347137/149‑1 (MQ=255)
                                                        gtgtCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGTATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGgtcgtc                                                                                             >  4:219668/1‑149 (MQ=255)
                                                        gtgtCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGTATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGgtcgt                                                                                              >  2:107293/1‑148 (MQ=255)
                                                           tctcAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAg                                                                                          >  7:35244/1‑149 (MQ=255)
                                                           tctcAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAg                                                                                          >  2:248144/1‑149 (MQ=255)
                                                             tcAGACACGGATATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGATGAATGGAGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGcc                                                                                        >  4:15596/1‑149 (MQ=255)
                                                               aGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGTATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGcc                                                                                        <  1:107293/147‑1 (MQ=255)
                                                               aGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGTATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGcc                                                                                        <  3:219668/147‑1 (MQ=255)
                                                               aGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCg                                                                                      <  7:221107/149‑1 (MQ=255)
                                                               aGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCg                                                                                      <  7:290438/149‑1 (MQ=255)
                                                                    cGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACa                                                                                 <  7:84141/149‑1 (MQ=255)
                                                                                ccGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAAt                                                                                                                  >  2:270025/1‑104 (MQ=255)
                                                                                ccGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAAt                                                                                                                  <  1:270025/104‑1 (MQ=255)
                                                                                   gAGGCGGACAGAGCCTTCAGGGACGGCCGGTGAACACCACGAGGAATTGTGGAGAACAGTGGATGCCTGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCGGGCAGGTCGTCAAGCCCGGTACAGGGGCAACGGATAcc                                                                  <  4:235127/149‑1 (MQ=255)
                                                                                       cGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCgttg                                                              <  5:394648/149‑1 (MQ=255)
                                                                                        ggACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCgttgt                                                             <  8:35244/149‑1 (MQ=255)
                                                                                                  ttCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCggg                                                   >  5:101063/1‑149 (MQ=255)
                                                                                                          cGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGGACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAgg                                           >  8:41980/1‑149 (MQ=255)
                                                                                                           gggCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATATGGGCTGGCAGGTCGTCACGCCCGGTACAGGGGCAACGGATACCGGTGTTATTACCGGGGCGGCAggg                                          >  4:317261/1‑149 (MQ=255)
                                                                                                           gggCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAggg                                          <  5:350315/149‑1 (MQ=255)
                                                                                                           gggCTGTGAACACCACGCGGACTGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAggg                                          <  1:259207/149‑1 (MQ=255)
                                                                                                            ggCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAgggg                                         >  5:152433/1‑149 (MQ=255)
                                                                                                             gCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAggggg                                        >  2:90958/1‑149 (MQ=255)
                                                                                                             gCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAgggg                                         >  5:215259/1‑148 (MQ=255)
                                                                                                                 tgAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGAc                                      <  3:68258/147‑1 (MQ=255)
                                                                                                                 tgAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGAc                                      >  4:68258/1‑147 (MQ=255)
                                                                                                                 tgAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCg                                    >  3:383336/1‑149 (MQ=255)
                                                                                                                   tacacCACGCTGAATGTGGGCGAACAGTGGAGGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGa                                  <  6:372987/148‑1 (MQ=255)
                                                                                                                       ccacGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCa                              <  1:248144/149‑1 (MQ=255)
                                                                                                                        cacGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATg                                >  7:37279/1‑146 (MQ=255)
                                                                                                                        cacGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAg                             <  2:213123/149‑1 (MQ=255)
                                                                                                                           gCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGaa                           >  5:169003/1‑148 (MQ=255)
                                                                                                                             tGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCggg                                                   <  3:365477/122‑1 (MQ=255)
                                                                                                                             tGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCggg                                                   >  4:365477/1‑122 (MQ=255)
                                                                                                                                 tggtggCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATTATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTg                    <  4:295130/149‑1 (MQ=255)
                                                                                                                                     ggCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATAc                <  6:101063/149‑1 (MQ=255)
                                                                                                                                      gCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATAcc               >  4:181966/1‑149 (MQ=255)
                                                                                                                                         aaCAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCggg            >  2:250218/1‑149 (MQ=255)
                                                                                                                                         aaCAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCggg            <  8:62461/149‑1 (MQ=255)
                                                                                                                                           cAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCa          <  6:169003/149‑1 (MQ=255)
                                                                                                                                           cAGTGGATGCAGGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCa          <  5:312218/149‑1 (MQ=255)
                                                                                                                                               ggATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGttt      >  5:81560/1‑149 (MQ=255)
                                                                                                                                                gATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGttt      <  4:383336/148‑1 (MQ=255)
                                                                                                                                                gATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTg     >  4:160259/1‑149 (MQ=255)
                                                                                                                                                 aTGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTGt    <  8:265308/149‑1 (MQ=255)
                                                                                                                                                  ggCATGAGGGGGTGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTGtt   <  7:174298/148‑1 (MQ=255)
                                                                                                                                                   gCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTGTTc  >  3:89649/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCCAACAAAACGACTGTCACCAGTGGTGGTCTTCAGAGAGTGAACCCCGGTGGAAGTGTCTCAGACACGGTTATCAGTGCCGGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGATACCGTTGTTAATACCGGGGCGGAAGGGGGACCGGATGCAGAAAACGGTGATACCGGGCAGTTTGTTC  >  NC_000913/2071865‑2072160

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