Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,208,927 A→T 5.2% S65T (TCT→ACT)  tfaE ← e14 prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,9270AT 5.2% 345.7 / 4.8 114S65T (TCT→ACT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base A (57/51);  new base T (2/4);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 114 total aligned reads displayed.

ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCT  >  NC_000913/1208781‑1209075
                                                                                                                                                  |                                                                                                                                                     
aTCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACg                                                                                                                                                    >  2:536723/1‑150 (MQ=255)
   gTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCa                                                                                                                                                 >  1:1199528/1‑150 (MQ=255)
    tATGGAGATAAAGGGGCGTTAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTGTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                <  1:190815/150‑1 (MQ=255)
    tATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                >  1:882357/1‑150 (MQ=255)
    tATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                                                                                                                                <  1:1153040/150‑1 (MQ=255)
     aTGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGa                                                                                                                                               >  1:1307429/1‑150 (MQ=255)
      tGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAt                                                                                                                                              >  2:1128936/1‑150 (MQ=255)
        gagaTATAGGGGCGATAGTGGTTGTATTTTGAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                            <  1:62706/150‑1 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCCGTTCTGCAGAttt                                                                                                                                            <  2:1237089/150‑2 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                            <  2:961702/150‑1 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                            >  2:754841/1‑150 (MQ=255)
        gagaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCGTTGATTTCATTTGATTCTCCAGTTTTGGTGCTATAGACCTTTTCACCGCTATTGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTg                                                                                                                                            <  2:1061969/150‑1 (MQ=255)
         agaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCCGTTCTGCAGAtttt                                                                                                                                           <  1:507036/150‑3 (MQ=255)
         agaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGc                                                                                                                                           <  1:768943/150‑1 (MQ=255)
          gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                                                                                                                          >  2:649697/1‑150 (MQ=255)
          gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                                                                                                                          <  2:1517376/150‑1 (MQ=255)
          gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                                                                                                                          <  2:854406/150‑1 (MQ=255)
            tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATTGTCTGGCACATATTCCCATGAGTTATAATCTGCCGTTCTGCAGAttttaac                                                                                                                                        <  2:212927/150‑6 (MQ=255)
            tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATa                                                                                                                                        >  1:794286/1‑150 (MQ=255)
             aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                                                                                                                       >  2:912686/1‑150 (MQ=255)
                ggggCGATAGTTGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGATTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCa                                                                                                                                    <  1:1395949/150‑1 (MQ=255)
                 gggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                                                                                                                                   >  1:207223/1‑150 (MQ=255)
                  ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                                                                                                                                  >  1:1168577/1‑150 (MQ=255)
                     gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                                                                                                                               <  1:1225131/150‑1 (MQ=255)
                        tgTGGTTGTTTTTTCAGGGTAATCACCCGGTTCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTTTAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                                                                                                                            <  1:748176/149‑1 (MQ=255)
                         gTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                                                                                                                           >  2:920610/1‑150 (MQ=255)
                             ttGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAc                                                                                                                       <  2:1436547/150‑1 (MQ=255)
                              tGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                                                                                                                      <  2:599299/150‑1 (MQ=255)
                                  ttttCTGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTTCTATAGACCGTTTCACCGAGATTGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg                                                                                                                 <  1:1210231/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg                                                                                                                 <  2:371022/150‑1 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg                                                                                                                 >  1:315664/1‑150 (MQ=255)
                                  ttttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg                                                                                                                 <  2:315664/150‑1 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGTTTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAgggg                                                                                                                <  1:148371/150‑1 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATTTGTACCAgggg                                                                                                                >  1:177135/1‑150 (MQ=255)
                                   tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAgggg                                                                                                                <  1:1028461/150‑1 (MQ=255)
                                    ttgaGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGc                                                                                                               <  1:1149145/147‑1 (MQ=255)
                                     tCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTAACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCa                                                                                                              <  2:741669/150‑1 (MQ=255)
                                             aTCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACa                                                                                                      <  1:495877/150‑1 (MQ=255)
                                              tCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAg                                                                                                     <  1:1447814/150‑1 (MQ=255)
                                               cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg                                                                                                    <  2:1307429/150‑1 (MQ=255)
                                               cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg                                                                                                    <  2:1115071/150‑1 (MQ=255)
                                                    ggAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACAt                                                                                               >  1:1694169/1‑150 (MQ=255)
                                                       gCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGcc                                                                                            >  2:1042334/1‑150 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCg                                                                                           <  1:1366245/150‑1 (MQ=255)
                                                        cTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCg                                                                                           >  1:808418/1‑150 (MQ=255)
                                                         tttGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCgg                                                                                          <  1:54848/150‑1 (MQ=255)
                                                            gATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAAt                                                                                       >  1:1560109/1‑150 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                   <  2:1364089/150‑1 (MQ=255)
                                                                tCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCg                                                                                   >  1:1116511/1‑150 (MQ=255)
                                                                  ttttGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCACCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGaa                                                                                 >  1:1034191/1‑149 (MQ=255)
                                                                     tgtttCTCCATTTTTGGTTCTATTGACCGTTTCCCCGCGTTGGTCTGGCCCATATTCCCATTAGTTAAAATTTGCAGAAAGGCAGATTGCATTACCAGGCTTATGTGTAACAGGGGCATATAAACAGGAACATGCCGGAATGCCGACAcc                                                                              <  1:1102019/147‑1 (MQ=255)
                                                                         tctcCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACg                                                                          >  1:1287952/1‑150 (MQ=255)
                                                                          ttctaGTTTTGGTGCTATATACCGTTTCACCGAGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGACTTATGTGTACCAGGTGAATCTAAACAGGAACATGCCGGAATGCCGACACCAACgg                                                                         <  1:84423/146‑1 (MQ=255)
                                                                          ctcCAGTTTTGGAGCTATAGACCTTTTCAACGCGATTGTCTGGAACATATTCCCATTAGTTAAAATCTTCATAACGCCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACgg                                                                         <  1:62885/150‑1 (MQ=255)
                                                                            cctgtattGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCAGATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAATGCATAACCACCCTTATGTGTACCAGGGGCAGCTAAACAGGCACATGCCGGAATGCCGACACCACCGGCa                                                                       >  2:828096/7‑150 (MQ=255)
                                                                            ccAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCa                                                                       >  2:27921/1‑150 (MQ=255)
                                                                            ccAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCa                                                                       >  2:1620477/1‑150 (MQ=255)
                                                                                 tttGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGata                                                                  <  2:1425510/150‑1 (MQ=255)
                                                                                      tGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACAACAACGGCAAGATATTCAt                                                             >  2:1030513/1‑150 (MQ=255)
                                                                                       gCTATAGACCGTTTCACCTCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGTCAGATTGCATAACCATCCTTATGTGTACCAGGGGCATCTAAACAGTAACATGCCGGAATGCCGACACCAACGGCAAGATATTCAtt                                                            <  2:819829/150‑1 (MQ=255)
                                                                                         tataGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGTGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCAtttt                                                          >  1:56975/1‑149 (MQ=255)
                                                                                         tataGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTg                                                          >  2:704164/1‑150 (MQ=255)
                                                                                            aGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                       >  2:430175/1‑150 (MQ=255)
                                                                                            aGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                       <  2:1554833/150‑1 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                                       >  2:1333855/1‑150 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                      >  2:785933/1‑150 (MQ=255)
                                                                                             gACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGt                                                      <  2:287940/150‑1 (MQ=255)
                                                                                              aCCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTg                                                     >  2:1084375/1‑150 (MQ=255)
                                                                                                 gTTTCACCGCGATGGTCTGGCACATATTCCCAAGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGaa                                                  >  2:180212/1‑150 (MQ=255)
                                                                                                   ttCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAat                                                <  2:750478/150‑1 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCACCGGCAAGATATTCATTTGAAGTGGAAata                                               >  1:179505/1‑150 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                               <  2:1352112/150‑1 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                               >  1:1622072/1‑150 (MQ=255)
                                                                                                    tCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAata                                               <  1:1435779/150‑1 (MQ=255)
                                                                                                     cACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAatat                                              >  2:578777/1‑150 (MQ=255)
                                                                                                       ccGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAatatat                                            >  2:1178790/1‑150 (MQ=255)
                                                                                                          cgATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTcc                                         <  1:362138/150‑1 (MQ=255)
                                                                                                            aTGGTCTGGCACATATACCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCg                                       >  2:1368896/1‑150 (MQ=255)
                                                                                                             tGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCgt                                      <  2:757673/150‑1 (MQ=255)
                                                                                                               gTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCgtgt                                    >  1:1256781/1‑150 (MQ=255)
                                                                                                                tCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACTGCAGATTGCATAACCATCATTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGtt                                   <  1:1521697/150‑1 (MQ=255)
                                                                                                                   ggCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATACCGACACCAACGGCAAGATATTCATTTGAATTGGAAATATATTCCCGTGTTTCa                                >  1:878079/1‑150 (MQ=255)
                                                                                                                     cacaTATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATACCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCAcc                              >  1:994284/1‑150 (MQ=255)
                                                                                                                      acaTATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCa                             >  2:1114365/1‑150 (MQ=255)
                                                                                                                       caTATTACCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCAt                            >  2:549204/1‑150 (MQ=255)
                                                                                                                          atTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGt                         >  1:452202/1‑150 (MQ=255)
                                                                                                                            tCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAg                       <  2:1694169/150‑1 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                    <  1:1128936/150‑1 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                    >  1:134647/1‑150 (MQ=255)
                                                                                                                               cATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTta                    >  2:860806/1‑150 (MQ=255)
                                                                                                                                    gTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCACGTGTTTCACCATCCTAGTTATAAAc               >  1:446265/1‑150 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTTCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACg              >  2:347948/5‑150 (MQ=255)
                                                                                                                                     atgcaaTCTGCCGTTCTGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATTTATTCCCGTGTTTCACCATCGTAGTTATAAACg              >  2:1323720/5‑150 (MQ=255)
                                                                                                                                         aaTCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGACCATGCCGGAATGCCGACACGAACGGCAAGATATTCATTTGAAGTGGAAAGATATTCCCGTGTTTCACCCTCGTAGTTTTAAACGGTaa          >  2:572734/1‑150 (MQ=255)
                                                                                                                                           tCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCTCGTGTTTCACCATCGTAGTTATAAACGGTAACa        <  1:1042334/150‑1 (MQ=255)
                                                                                                                                             tGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCCCCATCGTAGTTATAAACGGTAACATc      >  1:1442661/1‑150 (MQ=255)
                                                                                                                                               cAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATccc    >  2:510236/1‑150 (MQ=255)
                                                                                                                                                aGAACGGCAGATTTCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCATAGTTATAAACGGTAACATcccc   <  2:556115/150‑1 (MQ=255)
                                                                                                                                                aGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATcccc   >  1:348721/1‑150 (MQ=255)
                                                                                                                                                 gAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGCCACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCt  >  1:1617533/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTATAAACGGTAACATCCCCT  >  NC_000913/1208781‑1209075

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: