Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,073,225 T→C W563R (TGG→CGG)  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,073,2250TC98.5% 492.0 / ‑1.3 135coding (1687/3120 nt)fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (63/70);  minor base . (1/1);  total (64/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 135 total aligned reads displayed.

ACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGA  >  NC_000913/2073090‑2073360
                                                                                                                                       |                                                                                                                                        
acgctgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGATCGTCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  2:605651/139‑1 (MQ=255)
 cgctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGaa                                                                                                                                      >  2:149192/1‑139 (MQ=255)
 cgctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  2:1180211/138‑1 (MQ=255)
            aaCGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATa                                                                                                                           >  2:995688/1‑139 (MQ=255)
                   gTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCAc                                                                                                                    <  1:1094721/139‑1 (MQ=255)
                       cGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTg                                                                                                                <  1:500307/139‑1 (MQ=255)
                         tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGTACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATcc                                                                                                                                 >  1:717788/1‑120 (MQ=255)
                         tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGTACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATcc                                                                                                                                 <  2:717788/120‑1 (MQ=255)
                         tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGAt                                                                                                                            >  1:103531/1‑125 (MQ=255)
                         tCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGAt                                                                                                                            <  2:103531/125‑1 (MQ=255)
                           accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGa                                                                                                                             <  2:551725/122‑1 (MQ=255)
                           accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGACACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGa                                                                                                                             >  1:551725/1‑122 (MQ=255)
                                ggATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCggtg                                                                                                       >  1:402735/1‑139 (MQ=255)
                                    gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTg                                                                                                                <  2:787515/126‑1 (MQ=255)
                                    gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTg                                                                                                                >  1:787515/1‑126 (MQ=255)
                                    gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGtggtgg                                                                                                   <  1:149192/139‑1 (MQ=255)
                                     tCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  <  2:402735/139‑1 (MQ=255)
                                     tCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  <  1:27465/139‑1 (MQ=255)
                                       aTTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATg                                                                                                >  1:1078381/1‑139 (MQ=255)
                                       aTTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATg                                                                                                >  2:1385819/1‑139 (MQ=255)
                                           cTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGaa                                                                                                                                      >  1:458376/1‑97 (MQ=255)
                                           cTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGaa                                                                                                                                      <  2:458376/97‑1 (MQ=255)
                                           cTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCt                                                                                            >  2:516418/1‑139 (MQ=255)
                                            tCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTc                                                                                           >  2:1369809/1‑139 (MQ=255)
                                            tCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTc                                                                                           >  1:1187811/1‑139 (MQ=255)
                                              aGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                                         <  2:52125/139‑1 (MQ=255)
                                                cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTg                                                                                                                <  1:466125/114‑1 (MQ=255)
                                                cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTg                                                                                                                >  2:466125/1‑114 (MQ=255)
                                                cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  >  1:782951/1‑128 (MQ=255)
                                                cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  <  2:782951/128‑1 (MQ=255)
                                                cgcgGTACAGCCATGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGGCCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCAGAATATCCGCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGcg                                                                                       <  2:939307/139‑2 (MQ=255)
                                                          cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       >  2:203144/1‑81 (MQ=255)
                                                          cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  1:203144/81‑1 (MQ=255)
                                                          cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATcc                                                                                                                                 <  1:1497191/87‑1 (MQ=255)
                                                          cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATcc                                                                                                                                 >  2:1497191/1‑87 (MQ=255)
                                                          cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  >  2:151670/1‑118 (MQ=255)
                                                            cTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATc                                                                                                                                  >  2:1109657/1‑84 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  2:96811/76‑1 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       >  1:96811/1‑76 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATc                                                                                                                                  <  2:1474381/81‑1 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATc                                                                                                                                  >  1:1474381/1‑81 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCa                                                                                                                     >  2:1012122/1‑94 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCa                                                                                                                     <  1:1012122/94‑1 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                                        >  2:338833/1‑139 (MQ=255)
                                                               aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                                        >  2:521084/1‑139 (MQ=255)
                                                                       cGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCACGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                <  2:920519/139‑1 (MQ=255)
                                                                            gcaCTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                           <  2:1471481/139‑1 (MQ=255)
                                                                                 gtgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACg                                                                                                                        <  1:679995/73‑1 (MQ=255)
                                                                                 gtgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACg                                                                                                                        >  2:679995/1‑73 (MQ=255)
                                                                                  tgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  1:416611/57‑1 (MQ=255)
                                                                                  tgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       >  2:416611/1‑57 (MQ=255)
                                                                                   gCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGcaca                                                    <  1:338833/139‑1 (MQ=255)
                                                                                      gAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAggg                                                 <  1:999430/139‑1 (MQ=255)
                                                                                      gAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAggg                                                 <  2:29671/139‑1 (MQ=255)
                                                                                      gAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAggg                                                 <  1:1482633/139‑1 (MQ=255)
                                                                                            tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                >  1:34607/1‑118 (MQ=255)
                                                                                            tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                <  2:34607/118‑1 (MQ=255)
                                                                                            tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTc                                          <  1:521084/139‑1 (MQ=255)
                                                                                            tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACACGCACAGGGAAGTTc                                          >  1:107853/1‑139 (MQ=255)
                                                                                              ccATTGACCCTACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTTCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCa                                                            <  1:1020277/120‑1 (MQ=255)
                                                                                                 ttGACCCCACGAATGTCACTCTTGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGAc                                                                                              <  1:930281/83‑1 (MQ=255)
                                                                                                 ttGACCCCACGAATGTCACTCTTGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGAc                                                                                              >  2:930281/1‑83 (MQ=255)
                                                                                                    aCCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGc                                  <  1:996925/139‑1 (MQ=255)
                                                                                                      cccACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                <  2:1200998/108‑1 (MQ=255)
                                                                                                      cccACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                >  1:1200998/1‑108 (MQ=255)
                                                                                                        cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       <  1:799125/35‑1 (MQ=255)
                                                                                                        cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                       >  2:799125/1‑35 (MQ=255)
                                                                                                        cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                           <  2:244864/111‑1 (MQ=255)
                                                                                                        cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                           >  1:244864/1‑111 (MQ=255)
                                                                                                             aTGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGaa                         <  2:505485/139‑1 (MQ=255)
                                                                                                                tCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGt                      >  2:1456437/1‑139 (MQ=255)
                                                                                                                  aCTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGa                    <  1:1456437/139‑1 (MQ=255)
                                                                                                                   ctctCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAtt                                                                     >  1:1154701/1‑90 (MQ=255)
                                                                                                                   ctctCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAtt                                                                     <  2:1154701/90‑1 (MQ=255)
                                                                                                                     ctcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGaaaa                 >  1:533370/1‑139 (MQ=255)
                                                                                                                      tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGaaaaa                <  1:1385819/139‑1 (MQ=255)
                                                                                                                       cGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAAc               >  1:925421/1‑139 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                                                 <  2:229176/56‑1 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa                                                                               >  1:247223/1‑74 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa                                                                               <  2:247223/74‑1 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                <  1:1167788/89‑1 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                                                >  2:1167788/1‑89 (MQ=255)
                                                                                                                         ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCCGTACCGGCAACCCTGAAAGTGAAAAAcc              >  2:1047189/1‑139 (MQ=255)
                                                                                                                          cTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                                                 >  1:1246446/1‑55 (MQ=255)
                                                                                                                          cTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                                                 <  2:1246446/55‑1 (MQ=255)
                                                                                                                          cTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATg                                                                                    >  1:172973/1‑68 (MQ=255)
                                                                                                                          cTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATg                                                                                    <  2:172973/68‑1 (MQ=255)
                                                                                                                           tCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGa           <  1:1119043/139‑1 (MQ=255)
                                                                                                                           tCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGa           <  2:655110/139‑1 (MQ=255)
                                                                                                                           tCCGGTGCCACCCGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTAc                                      <  2:551504/111‑1 (MQ=255)
                                                                                                                           tCCGGTGCCACCCGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTAc                                      >  1:551504/1‑111 (MQ=255)
                                                                                                                             cGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAccc                                                         <  2:1484899/92‑1 (MQ=255)
                                                                                                                                tGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCg                                                                                                          <  2:1418920/40‑1 (MQ=38)
                                                                                                                                tGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCg                                                                                                          >  1:1418920/1‑40 (MQ=38)
                                                                                                                                tGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGa      >  1:40994/1‑139 (MQ=255)
                                                                                                                                  ccACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                  >  1:372447/1‑46 (MQ=255)
                                                                                                                                  ccACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACa    >  1:668458/1‑139 (MQ=255)
                                                                                                                                  ccACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACa    <  1:516418/139‑1 (MQ=255)
                                                                                                                                    aCCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGAc                                                                              <  1:182849/64‑1 (MQ=255)
                                                                                                                                    aCCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAggg                                                 <  1:1421442/93‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
ACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGA  >  NC_000913/2073090‑2073360

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