Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,065 G→C 6.7% D35E (GAC→GAG) ‡ ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0650GC 6.7% 417.2 / 15.3 132D35E (GAC→GAG) ‡ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (64/58);  new base C (7/2);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 135 total aligned reads displayed.

TCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGC  >  NC_000913/751916‑752214
                                                                                                                                                     |                                                                                                                                                     
tCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCgg                                                                                                                                                       >  2:696077/1‑150 (MQ=255)
         cATGGACTGGTTTTGTCACCCTTTTTTTGCCAGTTAATATTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGATAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAtt                                                                                                                                              <  1:126386/150‑1 (MQ=255)
          aTGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTg                                                                                                                                             <  1:1193610/150‑1 (MQ=255)
           tGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTTTTTTTTTTGTTATTCACCGCAACGCTAACAAAATATTCACATTTAGCAATGATTCCTTTTTCCTTTAATATTGAAATATCAATGCGGTCGCGCATTGa                                                                                                                                            <  2:165109/150‑1 (MQ=255)
           tGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGa                                                                                                                                            <  1:1657346/150‑1 (MQ=255)
           tGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCATTTAATATTAAAATAAAAATGAGTTCGCGAATTGa                                                                                                                                            >  2:62583/1‑150 (MQ=255)
            ggAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGAt                                                                                                                                           <  2:487623/150‑1 (MQ=255)
             gAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCACCGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                          <  1:1077944/150‑2 (MQ=255)
             gAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                          <  1:1114719/150‑2 (MQ=255)
                  gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                                     >  2:1604664/1‑150 (MQ=255)
                  gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                                     <  2:866413/150‑1 (MQ=255)
                   ttttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCga                                                                                                                                    >  2:369788/1‑150 (MQ=255)
                       gTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAAc                                                                                                                                >  1:496810/1‑150 (MQ=255)
                        tCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACa                                                                                                                               <  2:1397566/150‑1 (MQ=255)
                        tCACCCTTTTTTTGCCAGTTAAATTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGCCTCCTTTTTCCTTTAATAATGAAATATCAAGGCGGGCGCGCTTTGATTTATCGAGAACa                                                                                                                               >  2:74466/1‑150 (MQ=255)
                          aCCCTTATTTTGCCAGTTAATTTATTGCCCATTACTGATTTTGTTGTTATTCACCGAAACGCAAACAAAATATTCACCGGGAGCAATGACTCCTTTTACCTTTAATAATGAAATATAAATGCAGTCGAGCAATGATTTATCTAGAACAtt                                                                                                                             >  2:427156/1‑150 (MQ=255)
                            ccTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCa                                                                                                                           <  2:1424067/150‑1 (MQ=255)
                             cTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAg                                                                                                                          >  2:1513022/1‑150 (MQ=255)
                             cTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAg                                                                                                                          >  1:648790/1‑150 (MQ=255)
                                 ttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTg                                                                                                                      <  1:947762/150‑1 (MQ=255)
                                        aGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTGACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGa                                                                                                               <  2:923255/150‑1 (MQ=255)
                                         gTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGAc                                                                                                              <  1:1513022/150‑1 (MQ=255)
                                          ttAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGAcc                                                                                                             >  2:16982/1‑150 (MQ=255)
                                            aaTTTATTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGc                                                                                                           <  2:1584597/150‑1 (MQ=255)
                                                ttttGCCCATTACTGATTTTGTTGTTATTCACCGCAACTCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAtt                                                                                                       <  1:304088/150‑1 (MQ=255)
                                                ttttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAtt                                                                                                       <  1:1514143/150‑1 (MQ=255)
                                                ttttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAtt                                                                                                       <  1:463215/150‑1 (MQ=255)
                                                 tttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTg                                                                                                      >  2:228848/1‑150 (MQ=255)
                                                    gCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTgcgc                                                                                                   >  1:1366145/1‑150 (MQ=255)
                                                    gCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTgcgc                                                                                                   <  2:321779/150‑1 (MQ=255)
                                                       cATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAg                                                                                                <  1:530517/150‑1 (MQ=255)
                                                        aTTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTGGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                               <  2:1364191/150‑1 (MQ=255)
                                                        aTTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                               >  2:1410327/1‑150 (MQ=255)
                                                            cTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGAACGCCATTGCGCAAGTCAta                                                                                           >  2:71670/1‑149 (MQ=255)
                                                               aTTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                        <  2:1067074/150‑1 (MQ=255)
                                                               aTTTTGTTGTTATTCACCGCAACGCTAACAAAATAATCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                        <  1:234228/150‑1 (MQ=255)
                                                                   tgttgtTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                    >  1:348140/1‑150 (MQ=255)
                                                                   tgttgtTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                    >  2:1454837/1‑150 (MQ=255)
                                                                   tgttgtTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATAGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                    >  2:237545/1‑150 (MQ=255)
                                                                    gttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTATTTATGAAATATCAATGCGGTCGCGAATTGATTTATCGAGAACATTAAGATTGAATTCGACCGACATTGAGCAAGCCATCGCCATGAc                                                                                   >  2:830830/1‑150 (MQ=255)
                                                                    gttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAc                                                                                   >  1:1294039/1‑150 (MQ=255)
                                                                     ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAcc                                                                                  <  2:1069255/150‑1 (MQ=255)
                                                                       gttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCTCCATTGCGCAAGGCATCGCCATGACcag                                                                                >  2:1113318/1‑150 (MQ=255)
                                                                       gttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcag                                                                                <  2:64275/150‑1 (MQ=255)
                                                                        ttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcagg                                                                               <  1:228848/150‑1 (MQ=255)
                                                                         tATTCAAAGCAACTATAACAAAATATTCACCGGGAGCATTGACTCATTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggc                                                                              <  1:96101/150‑1 (MQ=255)
                                                                           ttCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcag                                                                            >  1:1179418/1‑150 (MQ=255)
                                                                            tCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcagg                                                                           <  2:1179418/150‑1 (MQ=255)
                                                                            tCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcagg                                                                           <  1:664910/150‑1 (MQ=255)
                                                                                  cAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaa                                                                     >  1:731663/1‑150 (MQ=255)
                                                                                     cGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAag                                                                  <  1:696077/150‑1 (MQ=255)
                                                                                      gCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAaga                                                                 >  1:548090/1‑150 (MQ=255)
                                                                                         aaCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagagag                                                              >  1:1434426/1‑150 (MQ=255)
                                                                                         aaCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagagag                                                              >  2:795676/1‑150 (MQ=255)
                                                                                           cAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGTAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTc                                                            <  1:422241/150‑1 (MQ=255)
                                                                                           cAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                               >  2:191670/1‑115 (MQ=255)
                                                                                            aaaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCg                                                           >  1:827271/1‑150 (MQ=255)
                                                                                            aaaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATAAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCg                                                           >  1:526042/1‑150 (MQ=255)
                                                                                               atatTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATa                                                        <  2:1398505/150‑1 (MQ=255)
                                                                                                tatTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                                >  2:550787/1‑141 (MQ=255)
                                                                                                tatTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                                <  1:550787/141‑1 (MQ=255)
                                                                                                tatTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATaa                                                       <  2:496810/150‑1 (MQ=255)
                                                                                                 atTCACCGGGAGCAATGCCTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCCCCATGACCAGGCAGGATACAAAAGAGAGTCGATaaa                                                      >  2:807305/1‑150 (MQ=255)
                                                                                                 atTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGGTACAAAAGAGAGTCGATaaa                                                      <  1:948056/150‑1 (MQ=255)
                                                                                                    caaccggACCAATGACTACTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAatat                                                   <  1:62583/145‑1 (MQ=255)
                                                                                                     aCCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATAtt                                                  >  1:770084/1‑150 (MQ=255)
                                                                                                      ccGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACACAAGAGAGTCGATAAATATTc                                                 <  1:71670/150‑1 (MQ=255)
                                                                                                      ccGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTc                                                 >  1:11574/1‑150 (MQ=255)
                                                                                                       cGGGAGCAATGACTCATTTTTCCTTTAATAATCAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCATATTGAATTCGACCGACATTGCGCAAGGCATAGCAATGACCAGGCAGGATAACAAAGAGAGTAGATAAATAttaa                                                >  1:1062487/1‑148 (MQ=255)
                                                                                                         ggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                              <  1:449440/150‑1 (MQ=255)
                                                                                                              aaTGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAAAATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTc                                         >  2:1253505/1‑150 (MQ=255)
                                                                                                                 gACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTTAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGTCATGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCAt                                      >  1:1212751/1‑150 (MQ=255)
                                                                                                                 gACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCAt                                      <  1:369788/150‑1 (MQ=255)
                                                                                                                     ccTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCt                                  <  2:264277/150‑1 (MQ=255)
                                                                                                                      cTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCAAGGTGTCCATACCTg                                 >  1:148862/1‑150 (MQ=255)
                                                                                                                       tttttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTTGATAAATATTCACGGTGTCCATACCTGa                                >  1:1231334/1‑150 (MQ=255)
                                                                                                                          ttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATaa                             <  1:856295/150‑1 (MQ=255)
                                                                                                                          ttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATaa                             <  1:1135941/150‑1 (MQ=255)
                                                                                                                           tCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATaaa                            >  1:1023048/1‑150 (MQ=255)
                                                                                                                             cTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATAc                                                                       >  1:319977/1‑105 (MQ=255)
                                                                                                                             cTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAata                          >  1:5176/1‑150 (MQ=255)
                                                                                                                               ttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGAGACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAtt                        >  2:1318387/1‑150 (MQ=255)
                                                                                                                                taataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAttt                       <  1:1410327/150‑1 (MQ=255)
                                                                                                                                    aatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATg                   <  2:446777/150‑1 (MQ=255)
                                                                                                                                    aatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATg                   >  2:598848/1‑150 (MQ=255)
                                                                                                                                       gataaatCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTAGACCGCCATTGCGCAAGGCATCGCAATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTAACGGTTTCCAAACCTGATAAATATTATATTaaa                >  2:1154779/6‑150 (MQ=255)
                                                                                                                                       gAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGaaa                >  2:634048/1‑150 (MQ=255)
                                                                                                                                        ataaatCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAg               >  2:602900/5‑150 (MQ=255)
                                                                                                                                        aaaTATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAg               >  2:1142741/1‑150 (MQ=255)
                                                                                                                                          atatCAATGCGGTCGCGCATTGATTTATCGAGAACAATAAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGATTCTATAAATATTCACGGTGTCAATAACTTATAAATATTTTATGAAAggc             >  2:172856/1‑150 (MQ=255)
                                                                                                                                           tatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcg            >  2:1014591/1‑150 (MQ=255)
                                                                                                                                           tatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcg            >  2:120363/1‑150 (MQ=255)
                                                                                                                                           aatCAATGCGCGACCGCATTGATTTATCGAGAACATTCTGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcg            >  1:213743/2‑150 (MQ=255)
                                                                                                                                           aatCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGTCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcg            >  1:1638406/2‑150 (MQ=255)
                                                                                                                                            atCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCGTTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcgg           <  1:551201/150‑1 (MQ=255)
                                                                                                                                            atCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcgg           >  1:565294/1‑150 (MQ=255)
                                                                                                                                            atCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcgg           >  2:1687509/1‑150 (MQ=255)
                                                                                                                                             tCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcggc          <  1:598848/150‑1 (MQ=255)
                                                                                                                                                 tGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCgatg      <  2:1294039/150‑1 (MQ=255)
                                                                                                                                                    ggTCGCGCATTGATTTATCTAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCATTATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATTTTTTATGAAAGGCGGCGatgatg   <  1:427156/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGC  >  NC_000913/751916‑752214

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: